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Synopsis 


In modern times, due to various environmental factors such as industrialization, pollution, 
deforestation, etc., the animal habitats are affected all around the world. The deforestation 
caused by industrial and human activities has a very detrimental effect on biological popula- 
tions in a forest habitats leading to loss of biodiversity, the important reasons being depletion 
of resource, patchiness, of habitats and population migration (diffusion). One example, where 
forest habitat has been destroyed by industrialization, causing patchiness, is the Doon Valley 
in Utter Pradesh, India, where ecological stability is threatened. Therefore, the study of pop- 
ulation dynamics with diffusion and explicit dependence of populations on resources in both 
homogeneous and patchy habitats is an important area of research in the field of mathematical 
ecology dealing with survival of species. 

It may be noted here that the study releted to the dynamics of two or more interacting 
populations in a patchy habitat with diffusion and resource dependence is very limited. Also, the 
effect of supplementary resource on population dynamics in patchy habitats is a new area 
of research in mathematical ecology. In view of the above, in this thesis, the effect of patchiness 
as well as of supplementary resource on a single species, prey-predator and competing species 
systems are investigated by using qualitative theory of differential equations. 

The work embodied in this thesis is divided into seven chapters. Chapter 1, is concerned 
with general introduction and literature survey relevant to population dynamics with diffusion 
in homogeneous and patchy habitats. The literature related to alternative or supplementary 
resource on interacting populations has also been reviewed briefly in this chapter. 

in chapter 2. a general single species harvesting model with diffusion in a linear two- 
patch habitat, under both the reservoir and no-flux boundary conditions and a continuous flux 
matching condition at the interface, is studied. It is shown that there exists a unique positive, 
continuous, monotonic steady state solution under both setsof boundary conditions. The 

linear and nonlinear stability behavoir of the non-uniform and uniform steady state cases, have 
been discussed and the corresponding conditions are determind. 



In chapter 3, a general single species model with diffusion in a linear two-patch "habitat, 
under both the reservoir and no-flux boundary conditions and a continuous flux matching con- 
dition at the interface, is investigated again by considering the dependence of the population 
on a supplementary resource. It is shown that the behavoir of the steady state solution is same 
as in chapter 2, but bit level of steady state distribution of the species population is higher at 
each location in the habitat due to the presence > of the supplementary self-renewable resource. 

In chapter 4, a prey-predator type model with diffusion in a homogeneous and two-patch 
habitats, under both the reservoir and no-flux boundary conditions and a continuous flux match- 
ing condition at the.interface, is studied. It is shown that the steady state solutions of both 
the species are positive, continuous and monotonic in the two-patch habitat. By comparing the 
results with the homogeneous habitat, it has been observed that the diffusion may increase the 
stability, but the patchiness destabilizes the system. 

In chapter 5, the same prey-predator system as in chapter 4, but with a supplementary 
self-renewable resource for the prey population in a two-patch habitat, is investigated. The 
behavoir of the steady state solutions have been shown to be the same as in chapter 4, but the 
level of steady state distributions of prey and predator population are correspondingly greater, 
in the presence of the supplementary resource for the prey. 

In chapter 6, a two species competition model with diffusion in a homogeneous and two- 
patch habitats, under both the reservoir and no-flux boundary conditions and a continuous 
flux matching condition at the interface, is studied. It is found that in this case also, similar 
as in chapter 4, the steady state solutions of both the species are positive, continuous and 
monotonic in the two-patch habitat. Further, by comparing the results with the homogeneous 
case, it is noted that in this case also the diffusion may increase the stability, but the patchiness 
destabilizes the system. 

In chapter 7, the same competing system as in chapter 6, with a common supplementary 
resource for both the species in two-patch habitat, is modelled and analyzed. The behavoir of 
the steady state solutions of both the species are same as in chapter 6, but the level of steady 
state distributions of both the populations are correspdndingly greater, in the presence of a | 
common self-renewable supplementary resource. 


11 



ACKNOWLEDGEMENT 


I express my heartfelt gratitude to my thesis supervisor, Professor J. B. Shukla, for his 
interest, constructive criticism, excellent rapport, indefatigable encouragement throughout the 
research program. I extend my thanks and profound respect to Professor H. I. Freedman, Applied 
Mathematics Institute, University of Alberta, Edmonton, Canada, for valuable suggestion and 
constructive criticism. 


I acknowledge with pleasure the cooperation and help received from all faculty members, 
Department of Mathematics, l.l.T. Kanpur. 


I would like to thank all my friends and well wishers especially Dipak, Pramananda, Alok, 
Arvind, Bhuvanesh, G.K., Pandayji, Johuri, Gautam, Paritosh, Swapan, B.P., Jal, Bishu, Indra, 
Shusmita, Mini and many others for their enthusiastic help throughout my stay at l.l.T. Kanpur 
making it a memorable and pleasant one. 


I take this opportunity to express my respect and thankfulness to my well wishers and 
colleagues Prof. J.K. Sharma, Prof. S.N. Saha, Dr. Sanjay Srivastava, Dr. Gagan Gupta, 
Mr. H.S. Payel, Mr. Narayan Misra, Mr. Surinder Singh and many others at Beant College of 
Engineering and Technology, Gurdaspur. 


Above all, I fell very much obliged to my parents and all others family members for what 
I have received from them by way of inspiration, love, encouragement and moral support. 



August, 1997 


JOYDIP DHAR 



Contents 



Contents 


List of Figures 


List of Tables 

1 General Introducti^j^ 

1.1 Introduction . 

1.2 Migration of Pon,,! 

l^-tion 

1.3 Patchy Habitat 

1.4 Effect of Suppletjj- 

^ Resource . . 

1.5 Summary . . . 


1 

1 

2 

6 

8 

10 


2 A Single Species 

tat ^ 'testing Model with Diffusion in a Two-Patch Habi- 

j . 21 

2.1 Introduction . . 

■■■••.. 21 

2.2 Mathematical 

2.3 Analysis of the Mocjgj 


2.3.1 Analysis , 

Reservoir Boundary Conditions 25 

2.3.2 Analysis 

No-Piux Boundary Conditions 31 



CONTENTS 


11 


2.3.3 The Case of Uniform Equilibrium State 35 

2.4 Numerical Examples 36 

2.5 Summary 44 

3 A Single Species Model with Diffusion and Supplementary Resource in 

Homogeneous and Two-patch Habitats 47 

3.1 Introduction 47 

3.2 Mathematical Model 49 

3.3 Analysis of the Model in a Homogeneous Habitat 52 

3.3.1 Model without Diffusion 52 

3.3.2 Model with Diffusion 53 

3.4 Analysis of the Model with Diffusion in a Two-Patch Habitat 54 

3.4.1 The Non-Uniform Steady State under Reservoir Boundary Conditions 54 

3.4.2 The Non-Uniform Steady State under No-Flux Boundary Conditions 60 

3.4.3 The Uniform Steady State under Both Sets of Boundary Conditions 63 

3.5 A Particular Case 65 

3.5.1 Analysis of the Model in a Homogeneous Habitat 65 

3.5.2 Analysis of the Model with Diffusion in a Two-Patch Habitat ... 67 

3.5.3 The Uniform Steady State under Both Sets of Boundary Conditions 74 

3.6 Numerical Examples 75 

3.7 Summary 79 

4 A Prey-Predator Type Model with Diffusion in Homogeneous and Two- 

Patch Habitats 82 

4.1 Introduction 82 

4.2 Mathematical Model 84 

4.3 Analysis of the Model in a Homogeneous Habitat 86 

4.3.1 Model without Diffusion 86 



CONTENTS 


111 


4.3.2 Model with Diffusion 95 

4.4 Analysis of the Model with Diffusion in a Two-Patch Habitat 98 

4.4.1 The Uniform Steady State under Both Sets of Boundary Conditions 98 

4.4.2 The Non-Uniform Equilibrium State 102 

4.4.3 The Model under Reservoir Boundary Conditions: When 

and y 2 >yl 103 

4.4.4 The Alodel under Reservoir Boundary Conditions: When > xj 

and yl> y 2 Ill 

4.4.5 The Model under No-flux Boundary Conditions: When > x^ and 

y^^Vi 112 

4.4.6 Both the Species have Uniform Equilibrium State in the Second PatchllG 

4.5 Numerical Examples 117 

4.6 Summary 125 

5 A Prey-Predator Type Model with DiflFusion and a Supplementary Re- 
source for the Prey in a Two-Patch Habitat 129 

5.1 Introduction 129 

5.2 Mathematical Model 131 

5.3 Analysis of the Model in Two-Patch Habitat 133 

5.3.1 The Non-uniform Steady State: Under Both Sets of Boundary Con- 
ditions 133 

5.3.2 The Uniform Equilibrium State: Under Both Sets of Boundary Con- 
ditions 140 

5.4 Numerical Example 144 

5.5 Summary 146 

6 A Two Competing Species Model with DiflFusion in Homogeneous and 

Two-Patch Habitats 149 

6.1 Introduction 149 

6.2 Mathematical Model 151 


CONTENTS 


IV 


6.3 Analysis of the Model in a Homogeneous Habitat 153 

6.3.1 Model without Diffusion 153 

6.3.2 Model with Diffusion 160 

6.4 Analysis of the Model with Diffusion in a Two-Patch Habitat 162 

6.4.1 The Uniform Equilibrium State under Both Sets of Boundary Con- 
ditions 162 

6.4.2 The Non-Uniform Equilibrium State 164 

6.4.3 The Model under Reservoir Boundary Conditions: When 

and yl> yl 165 

6.4.4 The Model under Reservoir Boundary Conditions: When > a:* 

andyl> y^ 175 

6.4.5 The Model under No-Flux Boundary Conditions: When > x\ and 

y*2>yl 175 

6.4.6 Both The Species have Uniform Steady State In The Second Patch: 180 

6.5 Numerical Examples 181 

6.6 Summary 185 

7 A Competing Species Model with Diffusion and a Common Supplemen- 
tary Resource for Both the Species in a Two-Patch Habitat 188 

7.1 Introduction - 188 

7.2 Mathematical Model 189 

7.3 Analysis of the Model in a Two Patch Habitat 192 

7.3.1 The Non-uhiform Steady State: Under Both Sets of Boundary Con- 
ditions 192 

7.3.2 The Uniform Steady State: Under Both Sets of Boundary Conditions 199 

7.4 Numerical Example 203 

7.5 Summary 206 



List of Figures 


2.1 The unsteady state and steady state solutions, for Di = 0.7, D 2 = 0 . 6 , 

K\ = 100, K 2 — 200, T\ = 0.2 and T 2 — 0.3 39 

2.2 The unsteady state and steady state solutions, for Di = 2.0, D 2 = 2.0, 

Ki = 100, K 2 = 200, Ti = 0.2 and ra = 0.2 40 

2.3 The steady state solution, when .^1 = 50 < J ?’2 = 100 , for Di = D 2 = 0.6 

and ri =r 2 = 0.08 41 

2.4 The steady state solution when Ki = 100 > K 2 = 50, for Di = D 2 — 0.6 

and ri = r 2 = 0.08 42 

2.5 Effects of diffusion coefficients Dj’s on the steady state distribution, when 

Ki = 50<K2 = 100, for fi = r 2 = 0.08 43 

3.1 The existence of {R*, K^) 51 

3.2 The steady state solutions, when KI > and Pi/ti = / 32 / 7’2 76 

3.3 The steady state solutions, when KI < K^ and PxJti = ^21^2 77 

3.4 The steady state solutions, when KI < K^ and y^i/ri 7 ^ ^ 2/^2 78 

4.1 The function H* with respect to a:* 91 

4.2 Existence of F* 94 

4.3 The steady state solutions, when x\ = 600 < Xj = 1000, y\ = 100 < 2/2 = 300.119 

4.4 The steady state solutions, when xj = 100 > xl = 60, yj = 80 > j/j = 50. . 120 



LIST OF FIGURES 


VI 


4.5 The steady state solutions, when = 50 < = 100, = 80 > ^2 = 40. . 121 

4.6 The steady state solutions, when rcj = 100 > = 40, yl = bO < = 80. . 122 

4.7 Effects of diffusion coefficients on the steady state solutions, for three set of 

values of T>y’s 123 

4.8 The Steady State distributions when the both the species have uniform equi- 
librium state {K 2 = X 2 = 120, M 2 = ?/| = 110) in the second patch 124 

5.1 The steady state solutions for both the species, with and without supple- 
mentary resource for the prey 145 

6.1 Existence of jE*, i.e. (a;*,y*) [for case (i)] 155 

6.2 Existence of i.e. {x*,y*) [for case (ii)] 156 

6.3 The steady state solutions, when X 2 = 162.162 > xj = 55.882, y^ = 

141.892 >y*i= 32.353 182 

6.4 The steady state solutions, when xl = 100.0 < = 149.99, = 74.999 < 

yl = 149.99 183 

6.5 Effects of diffusion coefl&cients in steady state solutions, when ajj = 55.882 < 

x’^ = 127.728, 2 / 1 * = 32.353 < 2/2 = 113.614 184 

7.1 The steady state solutions for both the species, with and without a common 

supplementary resource for both the species 205 


List of Tables 


3.3 Effects of various parameters on R* and K* 53 

4.3 Effects of various parameters on x* and y* 87 

6.3 Effects of various parameters on x* and y* 160 




Chapter 1 


General Introduction 


1.1 Introduction 


In modern times, due to various factors such as industrialization, pollution, deforestation 
etc., the animal habitats are affected all around the world. The deforestation caused by 
industrial and human activities has detrimental effects on biological populations in a forest 
habitat leading to extinction of certain rare species affecting biodiversity, the important 
reasons being depletion of resources, patchiness of the habitat and population migration 
(diffusion). One example, where forest habitat has been destroyed by industrialization, 
causing patchiness and resource depletion, is the Doon Valley in Utter Pradesh, India [118], 
where ecological stability is threatened. Therefore, the study of population dynamics in 
both homogeneous and patchy habitats with diffusion and explicit dependence of population 
on resource is an important area of research in the field of mathematical biology dealing 
with survival of species. 


The evolution and existence of species has been the subject of scientific investigation 
since the days of Darwin and studies were mainly in the form of experimental observations. 


The first major m?^ 
Volterra (1926), wh: 


lematical attempt in this direction was due to Lotka (1924), and 
1 is the main theme of deterministic theory of population dynamics in 
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mathematical biology. Over the last seventy years, many models for two or more interacting 
species have been proposed on the basis of Lotka and Volterra models by taking into account 
the effects of crowding, age structure, time delay, functional response, harvesting, toxicants, 
habitat changes, etc. [ Rescigno [105], May [85], Cushing [22, 23], Harada and Fukao [53], 
Gopalsamy [42, 43, 44], Hsu and co-authors [61, 62, 66, 67], Freedman [28, 29, 31], Cheng 
et al. [15], Freedman and Gopalsamy [33, 34], Yodzis [137], Hallam [49, 50, 51, 52], Kapur 
[74], Hutson and Viskers [70], Takeuchi [126, 127], Mitra et al. [91], Bonan [9], Shukla et 
al. [115, 118, 119, 120, 121], Shukla and Dubey [114, 116, 117] etc.]. 

However, in the original Lotka- Volterra model, spatial variation has not been con- 
sidered. It may be noted that Lotka- Volterra model focuses on population interactions 
at a point in the space ignoring movement (migration) which means a perfect mixing of 
the species in a given region. Mathematically, this is equivalent to assuming that the dis- 
persal rates are sufficiently high, and the population in the habitat are well mixed. It is 
nevertheless true that time and space are inseparable “sister coordinates” , and only when 
population dynamics of organisms are studied in both time and space can the ecological 
situation be understood in it’s reality [ Levin [82, 83, 84], Okubo [100] ]. 


1.2 Migration of Population 


As mentioned above, due to environmental factors and other related effects, the tendency 
of the species living in the habitat is to migrate to better suited regions for their survival 
and existence. In general the movement of the species in the habitat may arise due to 
certain factors such as overcrowding, anticlimate, predator chasing prey, refuge, and fugitive 
strategies and more importantly due to resource ( food ) limitation (i.e. availability ) in 
the habitat, Verma [131]. 

Experimental investigation of the phenomenon of animal dispersion was first con- 
ducted, for insects. Dobzhansky and Wright [24, 25] have done the famous experiments on 
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the release of Drosophila flies, and later Dobzhansky and Powell studied the rates of disper- 
sal of Drosophila pseudoobscura in the habitat. Some Japanese entomologists [ Watanabe 
et al. [133], Kono [77], Ito [71] and Morisita [92, 93] ], also contributed to the concept 
of biological diffusion by including interactive forces between dispersing individuals in the 
form of “density-dependent dispersion” . They emphasized that a model of dispersion must 
consider the forces operating between population individuals, and it cannot be limited to 
the simple random walk. Several other experimental studies have also been conducted on 
dispersion of Predator-Prey systems [ Okubo [100], Huffaker [68, 69], Bergerud et al. [6], 
Bergerud and Page [7] ]. 

A typical species may have dispersive migration in a particular direction due to favor- 
able environmental and ecological conditions in that direction but it may have no migration 
in other direction. The first successful attempt to study the migration of the species math- 
ematically is due to the Skellam [122] using the concept of random dispersal. Since then, 
several investigators [ Rosen [106, 107], Rothe [108, 109, 110], Fife [27], Shigesada [112], 
Okubo [100], Shigesada and Roughgarden [113], Allen [1, 2, 3], Cantrell and Cosner [11, 12], 
Pao [102, 103, 104] ], studied the effects of migration on local and global stability of inter- 
acting species system by considering Lotka-Volterra and other types of prey predator and 
competition models [ see also Levin [82, 83, 84] and Verma [131] ]. 

The effects of diffusion on population dynamics, since individuals interact and move 
about in the habitat, can be described by the continuous space-time interaction-migration 
models. The principal ingredients of these models are equations of the form 

^ = F{u) + DAu ( 1 . 1 ) 

where u is an n- vector (each component represents the density of one species), D is a matrix 
(diffusion coefficients matrix), and A is the Laplace operator in the spatial coordinates. The 
vector F is a catch-all term describing all reactions and interactions. Using (1.1) various 
population models vsrith diffusion for single species and interacting populations in both 
homogeneous and heterogeneous environment, have been proposed and analyzed, [see the 
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review article by Levin [83, 84] and McMurtrie [86]]. Several good examples of this type can 
be found in the monographs by Fife [27] and Okubo [100], which also provide a overview 
of the study of self-organizing phenomena in biology. 

Some workers have also studied the coexistence, persistence and extinction in single 
species and the Lotka^Volterra reaction diffusion models [ Gopalsamy [42], Allen [1, 2, 3], 
Mimmura et al. [88] ] and the global stability in generalized Lotka-Volterra with diffusion, 
Takeuchi [126]. 

Hadeler and Rothe [48] have studied the effects of dispersion on linear stability of 
interacting species in the one dimensional homogeneous finite habitat for prey predator 
system where the dispersion coefficients are equal and constant for both the species. It 
has been shown that the effect of such dispersion is to stabilize the equilibrium state 
under nonhomogeneous boundary conditions. Similar analysis has also been carried out 
by Gopalsamy [42] for competing species in case of homogeneous finite and semi-finite 
habitats. It is evident from his work that dispersive migration has no effect on the otherwise 
unstable equilibrium state in the case of semi-infinite habitat while migration may stabilizes 
the equilibrium state in the case of finite habitat. For unequal but constant dispersion 
coefficients of the two species, it has been shown that for competition model diffusion 
may increase the stability of equilibrium state (at least non decreasing) , but in the case of 
prey-predator model diffiision instability occurs, Mimura and Nishida [89], and Jorne [73]. 

In general, a diffusion process in an ecosystem tends to give rise to a uniform density 
of population in the habitat. As a consequence, it may be expected that diffusion, when 
it occurs, plays the general role of increasing stability in a system of mixed populations 
and resources. However there is an important exception, known as “diffusion induced 
instability” or “diffusive instability” . This exception might not be a rare event especially 
in prey-predator systems [129, 136]. 
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Levin [83] has given two examples of prey-predator model where he studied the role 
of diffusion by suitably choosing function F(u) in equation (1.1). 

Model 1: Lotka-Volterra Model (5i:prey, 52:predator) 

Fi = aiSi — P1S1S2 and F2 = — q;2‘S2 + P2S1S2, (1-2) 

where ai, pi, 0:2, P2 are positive constants. In this model no diffusive instability occurs. 
Model 2: Phytoplankton-Herbivore System (Levin and Segel, 1976) 

jFi = cbiSi -)- eSi — biSiS2 and F2 = — CS2 + ^28182, (f-3) 

this system is seen to have a diffusive instability if 6162 > ce, 62 > e, and furthermore 



Hastings [56] has derived the sufficient conditions for global stability in n-species 
Lotka-Volterra model with variable diffusion coefficients in a linear habitat (i.e. 0 < a: < L) 
under no-flux boundary conditions. He studied the following model, 


m 

dt 


Ni 




■ n 


(1.5) 


i=i 


where the DiS are measures of dispersal rate and may depend on space variable x and time 
t or any of the Ni- The following boundary conditions, 


8 N- 

-7—^ =0, for t > 0, a; = 0, a; = L, (1.6) 

ox 

associated with equation (1.5), imply that the region is truly closed. Hastings [57], also 
studied a single species model with large diffusion rates and established a necessary and 
sufficient condition for stability. 


Nallaswamy and Shukla [98], consider a dispersal prey-predator model with functional 
response. 

P 7VL / \ I \ 

(1.7) 


dNi 


dt 


^ o,i 2N2 \ d /_ dNi\ 

iVi (a, - aniV, - ^ (A-gj-j 
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dN2 

dt 




Q2iA^i >1 d I dN2 
+ aNj dx V " dx 


( 1 . 8 ) 


It is shown by analyzing this model that the functional response h^.a destabilizing effect. 
They have shown that there is no diffusive instability, even with functional response. They 
further established that if the equilibrium state is linearly stable, a subregion of the positive 
quadrant can be found in the phase plane where it is nonlinearly stable with or without 
diffusion. Shukla et al. [119] also studied the effects of dispersal in a multi-species food 
chain model. 


Timm and Okubo [129], considered the Levin-Segel model [84], and extended it to 
plankton populations with time-varying diffusivities. The local stability of uniform equi- 
libria is examined both anal 3 d;ically and numerically. They found that diffusive instability 
is less likely to occur in systems with time-varying diffusivity than those with constant 
diffusivity. 

Recently, Pao [102, 103, 104], studied the reaction diffusion equations with nonlocal 
boundary and initial conditions and, Takeuchi and Lu [128] considered a diffusive compet- 
itive Lotka-Volterra model, and studied the permanence and global stability of the system. 


1.3 Patchy Habitat 

The real habitats are both homogeneous and heterogeneous as regards to their environ- 
mental and ecological characteristics. Heterogeneity czm arise due to topographical, geo- 
graphical and other environmental conditions. Seasonal climatic changes can also affect 
the heterogeneity of the habitat. In general in a habitat, patchiness, a type of heterogeni- 
ety, occurs due to discontinuous variations in geographical, ecological and environmental 
(including climatic changes) characteristics, Verma [131]. As regards to the forest habitats, 
patchiness can arise due to industrialization, pollution and population. 
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Thus, in reality, one may visualize a habitat to be patchy if its ecological character- 
istics are constants but different in each patch. For example, the patchy distributions of 
plankton in the sea and lakes have been well-documented [ Cassie [21], Steele [124] ]. The 
mechanism which maintains this patchiness is still, however, a subject of controversy. A 
few of the various possibilities that have been proposed are : (1) behavioral reaction to or 
association with, temperature, salinity, and nutrient distributions, (2) food-chain associa- 
tion in predator-prey relationships, (3) aggregative behavior for breeding and feeding, (4) 
deforestation of the habitat caused by industrialization, pollution and population growth 
[100, 118]. 

Levin [82, 83, 84] studied the effects of patchiness in a heterogeneous habitat by con- 
sidering the movement of species from one patch to other as proportional to the difference 
of densities of the respective species in adjoining patches. The effect of population diffusion 
on coexistence of species has also been studied in patchy habitats by using both continuous 
diffusion and discrete diffusion models [ Allen [1, 2, 3], Freedman [28, 30, 32], Takeuchi 
[126, 127], Takeuchi and Lu [128]]. In these studies, it may be noted that the diffusion is 
considered only across the patches but with in the patch no diffusion of populations are 
considered. 

Verma [131] studied the linear and non-linear stability of uniform steady states corre- 
sponding to Lotka-Volterra type competitive species and prey-predator systems in a habitat 
with diffusion in each of the n-adjoining linear patches. She has shown that the uniform 
state which is otherwise stable in a homogeneous habitat, also stable in a patchy habitat 
but with restrictive conditions. In the case of a competition model, she has found that 
if the equilibrium state is unstable without diffusion, then it can still become stable with 
diffusion but under more restrictive conditions in the case of a patchy habitat. In the case 
of prey-predator system the stable equilibrium state with diffusion remains stable in the 
case of patchy habitat also. 
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Freedman et al. [38] proposed a single species diffusion model in a two-patch environ- 
ment. It has been shown that there exists a positive, monotonic, continuous non-uniform 
steady-state solution with continuous flux, under both the reservoir and no-flux bound- 
ary conditions, that is linearly asymptotically stable. Freedman and Wu [40], extended 
the study of single species model with three-patch environment and analyze the existence, 
piece wise monotonicity and stability of the steady states with only reservoir boundary 
conditions. In this paper they use an interesting theorem for topology, namely Topological 
Translation Theorem, which reduces a “three-patch” problem to “two-patch” problem and 
“two-patch” to “one-patch”. In the same year Freedman and Kriszti [35] also studied a 
class of single species models with diffusion in a habitat with number of patches. It has been 
shown that under certain conditions, there is a unique, positive, globally asymptotically 
stable steady state. 

It may be noted here that in the above studies, the stability of non-uniform steady 
state corresponding to systems of two or more populations having general interaction, in 
a patchy habitat with diffusion has not been studied. In view of these, in this thesis we 
study the effect of patchiness on the following systems: 

• A general logistic type single species system with harvesting and with diffusion in 
each patch. 

• A general prey-predator system (here predator is partially dependent upon a prey 
population) with diffusion in each patch. 

• A general logistic type two competing species system with diffusion in each patch. 


1.4 Effect of Supplementary Resource 


In the above we have discussed the growth and existence of population with or without 
diffusion in a habitat which may be homogeneous or patchy. However, in the real world the 
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explicit dependence of population on resource or nutrient plays an important role in the 
growth and existence of these populations. Some laboratory experiment were conducted in 
[99, 59, 134] on micro-organisms using chemostat to study competition between different 
populations of microorganisms for a growth limiting nutrient or resource. The mathematical 
analysis of competing species and several non-interacting populations which depend upon 
growth limiting nutrient in a chemostat with constant input and variable washout rate 
have been studied by taking into account Michaelis-Menten kinetics, Hsu et al.[60, 65] and 
Butler et al. [10]. Some other investigations have also been conducted for the two competing 
population when the populations are wholly dependent on a self-renewable resource without 
diffusion in a homogeneous habitat [63, 45, 91]. Recently Freedman and Shukla [37] have 
studied the effect of an alternative resource for predator on a prey-predator system with 
diffusion in a homogeneous habitat. 

It may be noted here that the above investigations been conducted either without 
diffusion or with diffusion in a homogeneous habitat and the effect of patchiness of the 
habitat has not been considered. In this thesis, therefore, we investigate the non-uniform 
steady state distribution and its stability system by considering diffusion and supplementary 
resource in two-patch habitat in the following cases: 

• A general single species model with a self-renewable supplementary resource. 

• A general prey-predator type model with a self-renewable supplementary resource for 
the prey population. 

• A general two species competition type model with a common self-renewable supple- 
mentary resource for both the species. 
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1.5 Summary 

This thesis consists of seven chapters and the organization of the thesis is as follows. 
Chapter- 1 deals with general introduction and literature survey relevant to population 
dynamics with diffusion in homogeneous and patchy habitats. In chapter-2, 4 and 6, 
the effect of diffusion in both homogeneous and two patch habitats have been studied for 
single species with harvesting, prey-predator system and the two species competition model 
using generalized growth rates and interacting functions. In chapter-3, 5 and 7, the role of 
supplementar}' resource in a patchy habitat have been modelled and analyzed in the three 
respective cases of chapter-2, 4, 6. 
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Chapter 2 


A Single Species Harvesting Model 
with Diffusion in a Two-Patch 
Habitat 


2.1 Introduction 

Biotic populations are usually distributed heterogeneously in their habitat, and the dis- 
tribution is often patchy due to patchiness of the habitat which arises from a variety 
of mechanisms and processes under various conditions including deforestation in the case 
of a forest habitat. It would, thus, seem natural to study the population dynamics of 
a single species by including diffusional effects, in a patchy habitat. Many investigators 
[4, 9, 11, 12, 13, 14, 15, 17, 18, 20] have shown that in a homogeneous habitat, the diffusion, 
increases the stability in a system, but this may not always be true, if the habitat is patchy 
[7, 8, 10, 16]. 

A model of a single species population living in two patch habitats with migration 
between them across a barrier was proposed by Freedman and Waltman [11]. The model 
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was extended by [6, 8] to include the case where species (animals) leaving one habitat does 
not necessarily reach the other habitat, the existence of a positive equilibrium as a function 
of barrier strengths was examined. 

Also Freedman et al.[10] studied a single species diffusion model by assuming that 
the habitat consists of two patches and has shown that there exists a positive, monotonic, 
continuous nonuniform steady state solution that is linearly asymptotically stable under 
both reservoir (or Dirichlet) and no-fiux boundary conditions but it is not globally stable 
because of patchiness. 

Keeping the above [6, 10, 11] in view, in this chapter, we consider the dynamics of a 
diffusing single species population undergoing harvesting in two adjoining patches with a 
continuous flux matching condition at the interface. The existence of positive, monotonic, 
continuous nonuniform steady state solution with continuous flux, under both reservoir and 
no-flux boundary conditions are studied. Also the stability of both linear and non-linear 
systems are discussed. This model is proposed keeping in view the patchiness caused by 
many man made projects in the Doon Valley (India) [19]. The model is also applicable to 
the Burwash caribou herd, which lives on both sides of the Shakwak Trench in the Kluane 
mountains of the Yukon territories, Canada. The dynamics of this herd has been discussed 
in [3]. 


The organization of this chapter is as follows. In next section we describe the general 
mathematical model for two patches. In section 2.3 we analyzed our model under reservoir 
and no-flxix boundary conditions in subsection 2.3.1 and 2.3.2 respectively, and the imiform 
steady state case is discussed in subsection 2.3.3. Finally in section 2.4 we have numerically 
studied the the unsteady and steady state solution of the system with reservoir boundary 
conditions. 
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2.2 Mathematical Model 

Consider a habitat consisting of two distinct adjoining or patches labeled by the running 
variable z = 1, 2. There are no restrictions on the way the habitat is partitioned into these 
patches. They may be identical or different in size; they may represent identical regimes of 
climate and of soil or nutrient factors, or they may represent distinct such environments, 
or might simply represent a partitioning of a homogeneous environment. The population 
dynamics of a single species of logistic type is given by the form of a system of parabolic 
partial differential equations in the two patches as follows: 

= N,{s, t)g,{N({s, t)) - Kp,(Ni(s, t)) + (2.1) 

z = l,2 0 <s<L2 

Here the growth rate function gi{Ni) and the harvesting function Pi{Ni) are such that 

(H.l): giiNi),pi{Ni)eC^[0,oo), 

gi(0)>0, andViVi>0, g-(iVi) < 0 

Pi(0) = 0, andViVi>0, P-(iVi) > 0. 

When the habitat has carrying capacity Ki in the i-th patch respectively, then gi{Ki) = 0. 

Further, we assume that 
(H.2): 3 ii:; > 0 such that - PiiK^) = 0 

It follows that K* < Ki and the equality only holds when Pi{Ni) = 0. 

Since \/Ni # K*, 

gi(jVi) - pi(jVi) gi(jVi) - Si{K;) vm) - Pim __ ,,, , 

n,-k; N<-Kt PM<i<'>. 

where min{Ni,K^} ^ max{Ni,K*}. Hence, 


(A - - Pi(JV()l < 0, VNt ^ k;. 


( 2 . 2 ) 
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At the interface s = Li, the continuous flux matching and continuity conditions are 
viuitten as 


D, 


ds = - 02 -^ 57 ^ and = 


(2.3) 


The model is studied under two sets of ( namely, reservoir and nofiux ) boundary 
conditions. In the case of reservoir boundary conditions, we take 


iVi(0,t) = ir* and N 2 {L 2 ,t) = K; 


(2.4) 


In the case of no-flux boundary conditions, we take 

dN-,{0,t) _ dN2iL2,t) 

ds ds 


(2.5) 


Finally, the model is completed by assuming some positive initial distribution, that is, 


Ni{s, 0) = Xi(s) > 0, Li_i <s<Li, i = 1, 2, 


( 2 . 6 ) 


such that A^i(0, 0) = K*, 
follows: 


and A2(T2,0) = K^. An example of Xi(s) naay be consider as 


, , {K^ - Kl) 

Xi{s) = Kl -h 


The existence of non-negative, smooth solutions of the above type of system in a 
homogeneous habitat is shown in [18]. We also note that in the case of no diffusion, the 
behavior of solutions of our model is well known (see [5], Chapter 1). At each s in the z-th 
patch, Ni{s,t) — )■ AT* as t — >• oo. Clearly, ^ then a constant solution of (2.1) 

in s is impossible and the problem must be analyzed by using the techniques of two point 
boundary value problems [1, 2]. 
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2.3 Analysis of the Model 

2.3.1 Analysis under Reservoir Boundary Conditions 

In this section we show the existence of non-uniform steady state solution of our model 
(2.1), (2.3), (2.4), and (2.6), following Freedman et al. [38]. This steady state solution is 
given by 

+ Uigi(ui) - UiPiiui) = 0 i = 1, 2 (2.7) 

The reservoir boundary conditions are given by 

ui(0) = A*and U2{L2) = (2.8) 

At the interface s = Li, we have, 

dUi{^L\) dU2{Lii) fr\ f T \ n^ 

Di — — — = A — — U2[Li) (2.9) 

We use the result of Freedman and Kriszti [7] to show the existence of the positive 
steady-state solutions in the form of the following theorem. 


Theorem 2.3.1 The steady-state problem (2.7), has a unique positive solution u. I 

Since the functions, Gi{ui) = gz(ui) - Pi(itt) are such that Gi(0) > 0, G((ui) < 0, 

for Ui > 0 and Gi{K*) = 0, therefore, Gi(u) < 0, if u > max and Gi{u) > 0, 

if 0 < M < min {K*,K 2 }. Thus the positive unique steady-state solution u of (2.7) is 

respectively concave or convex on connected subsets of 
{s e ( 0 ,L 2 ) : 5 7 ^ 0, Iri,L 2 ; u{s) > max 

or {s G (0, A) : s 7 ^ 0, Li,L 2 ’, 0 < u{s) < min 

Hence the unique solution u of (2.7) satisfies 

min {Jf*, < u{s) < max {K^, 
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We now consider without loss of generality 0 < Ki < ii'|. Therefore < u < K^. 

Let pi{s, Qij), Li-i < s < Li, are the unique solution of the equation (2.7), for i=l,2 
and pi{s, ctj) are such that 

-^(0, ai) = Oi, Pi(0, q:i) = 

dp 

-^{L2,a2) = a2, p2{L2,ot2) — K^. 

Multiplying both sides of (2.7) by 2dui/ds, and integrating from 0 if i=l and from L 2 
if i=2, we get 

r T 2 

(j'lj . 0 rui(s) 

-al = -—J^^ 6(s)[gi(^t(s)) -Pi(^i(s))]d^i (2.10) 

Using (2.10), we now prove the following lemmas. 


Lemma 2.3.1 If ai> 0 , then 


dpi{s,ai) 

ds 


> ai on 0 < s < Li- 


Proof: Prom (2.10) , we get, 

' dpi{s, ai) ~ 
ds 



2 rpi{s,ai) 

jr L &[gi(fi)-Pi«i)Ki 

1J\ JKl 


(2.11) 


since 

then there exists Si > 0 such that pi{s, ai) > iiTJ on 0 < s < Si. If not, let Sq, 0 < Sq < Li, 
be the first positive value, if it exists, such that pi(so, ai) = Kf. Then by the mean value 
theorem there exists s such that 0 < s < so and dpxis, a\)/ds = 0; that is. 


9 /*Pi ( 5 , 0 : 1 ) 

ei[gi(a)-Pi(ei)Ri. 

Di Jki 


( 2 . 12 ) 
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Now since pi(s, a^) > fov 0 < s < s and from (2.2), we have, gi(fi) - pi(5) < 0. 
Hence the right hand side of (2.12) is negative, giving a contradiction. Therefore pi(s, ai) > 
K* and dpi{s,ai)/ds > 0. Now, noting that {dpi{s,ai)/ds)^ is an increasing function of 
Pi , the lemma follows. I 


Lemma 2.3.2 If Q <p 2 < and a 2 > 0, then 


dp2is,a2) 

ds 


> Ot2, 


Li ^ S 1/2- 


Proof: From (2.10),we have. 


dp2{s,a2) 


-l2 


ds 


= q| 


2 fP2(.s,a2) 


Do Jk. 


rp2(s,a2) 

/ 6 [§ 2 (^ 2 ) — P 2 ( 6)]<^^2 

•/ KX 


and since 0 < p2 < hence from (2.2), 6[g2(6)“P2(6)] > 0- Hence dp 2 {s, Oi 2 )/ds > 0 . 2 , 
Li < s < L2 for 0 < p2 < ^ 2 - ® 


Lemma 2.3.3 Define Fi{ai) by Fi{o:i) =pi{Li,ai). Then there exists dj > 0, such that 

Fi : [0,di] [Kl K;] 

F 2 : [0,6:2] -t [K* 2 , K*] 


Proof: Note that Fi(ai), if it exists, is continuous. The existence of 6:1 > 0 follows from 
the following, we note that Fi(0) = K^, Fi(a:i) > iif* if Oi > 0 (by Lemma 2.3.1), and that 
Fi(oo) = 00. Hence we choose ojj > 0 to be the least value of Oi such that jPi(a!i) = K^. 
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To show that dt^ exists, we note that ^ 2 ( 0 ) = K^. By continuity, exists for 

sufficiently small 012 , 0:2 > 0. From Lemma 2.3.2, F 2 [{K^ - Kl)/{L 2 - Li)] is either less 
than Kl or does not exist. Therefore, there exists at least an 0 : 2 , such that da = iience 
the lemma. I 


Theorem 2.3.2 There exists a continuous, monotonic solution of system (2.7) with con- 
tinuous flux at Li. 


Proof; Lemmas 2.3.1 and 2.3.2 follow that any solution we construct must be monotonic. 
By Lemma 2.3.3, for each 0 < 0:2 < 0 : 2 , we can find an such that 0 < Oi < ai for which 
Pi{Li,ai) = p 2 (Li,a 2 ). Hence ai can be solved as a function of aa, <xi = h{a 2 ), to give a 
continuous solution of (2.7) with (2.8) and (2.10). 

Let 


^ r. dpi{Li,h{a2)) ^ dp2{Lua2) 
5{a,) = D, D, g- 

Clearly Q{ot 2 ) is continuous on 0 < 0:2 < da- Then we have 


^(0) = A 


9p\ (Li, di) 

ds 


> 0 , 


and 


Hence the theorem. 


Q{&2) = -D 2 


dp2{Li,a2) 

ds 


< 0 . 


I 


Now we find the conditions for asymptotic stability of the non-uniform steady state 
in both linear and non-linear cases. 


Theorem 2.3.3 The steady-state, continuous, monotonic solution of the system (2.1), 
with continuous flux matching condition at the interface (2.3), and under reservoir boundary 
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conditions (2.4), is locally asymptotically stable if 

d 

(«2[g2(“2) - P2(«2)]) <0, K; <U2< 


(2.13) 


Proof: Linearizing (2.1) by using 

Ni{s,t) = Ui{s) +ni{s,t) i = l,2 (2-14) 

We have 

f 9 /I • 

— = ni [gt(ui) + Uig'iiui) - pi{ui) - UiP-(ui)] + (2.15) 

Using (2.14), the corresponding reservoir boundary conditions and the matching conditions 
at the interface are as follows : 


ni(0,t) 


7^2 (A, t) = 0 


Di 


ni{Li,t) 


(2.16) 


Ss ' ” ' ‘da 

Now we consider the following positive definite function {Lq = 0), 

fLi 1 


j_l JLi-i Z 


(2.17) 


Differentiating (2.17) and using (2.15) and (2.16), we get 

ri^i 


V{t) = nl[gi{ui) + uig[{ui) -pi{ui)-uip[{ui)]ds + Di ni-^ds 
?^2[g2(w2)+W2g2(«2)-p2(M2)-«2P2(w2)]ds + Ay^ 


Now, considering the following integral and using (2.16), we get 


J — Di / ni 2 d^ + A / ^ 

Jo os JLi 

= D,n,{L„ - D,n,(Li,t)?^(L,,t) 

i: &]'<->■ c {iff- 


^2 


ds 


(2.18) 
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Further since for 0 < s < Li, ui > K^, Gi(wi) = gi(ui) - Pi(ui) < 0, and Gi(ni) = 
g'liui) - pi(wi) < 0, the expression gi(ui) + uig[{ui) - Pi(ui) - uipi(ui) is negative. 

Hence V{t) is negative definite if the condition (2.13) holds true, thus proving the 
theorem. I 

Remark: For Li < s < L 2 , < K^-, G 2 (u 2 ) = g 2 (« 2 ) - P 2 (w 2 ) > 0, and G 2 (n 2 ) < 0. 

This implies that the expression in the third integral is negative for very close to K^. 
Since it is smooth, it remain so V u^-, Kl<U 2 < by (2.13). 

Finally we state and prove the nonlinear stability theorem. 

Theorem 2.3.4 The steady-state, continuous, monotonic solution of the system (2.1) with 
continuous flux matching condition at the interface (2.S), and under the reservoir boundary 
conditions (2.4), is non-linearly asymptotically stable in the subregion of K* < Ni, Ui < K 2 , 
for i-l,2, if, 


(NigiiNi) - Uigflui)) {NiPi{Ni) - Uijpi{ui)) 


Ni Ui 


Ni - Ui 


< 0 , 


(2.19) 


Proof: From (2.1), (2.7) and (2.14), we have 

^ = (A'igi(JVi) - - (NMNi) - + A-jy 

Now consider the same positive definite function as given in (2.17), 

2 rLi 


^ Jb' ^ L' 

Uifds = Y^ f njds 


Therefore, 


( 2 . 20 ) 


V{t) 


n? 

'{Ni)gi{Ni) - Uigi{ui) 

{Ni)pi{Ni) - UiPiiuiY 

'Li-i * 

Ni - Ui 

Ni-Ui 


ds 
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^ rjCt d^Th' 

+ Ea/ n,^ds 

^ JLi-r ds^ 

and the second integral is always negative by (2.18). 


( 2 . 21 ) 


Therefore < 0 if the conditions (2.19) are satisfied. 


Remark: Using the mean value theorem, 

(•^i)§i(Ai) tiigi('lti) (.^t)Pi(.^t) WtPi(ti'i) 
Ni - Ui Ni - Ui 


+ iiMi) - M^i) - ^iPi(6) (2.22) 


where is a fixed point for i = 1,2, lies between iV,- and Ui. Since < Ni,Ui < i^|, 
therefore < K^. Again from (H.l), Pi(^i) > 0, p-(^i) > 0 and g-(^i) < 0,V^i > 0. 

Also gi(^i) < 0, only when > Ki, therefore the conditions (2.19), for i = 1 is always true 
if both ui,Ni > Ki > K^- Further from (2.22) for i = 1, 2, 


0 > -'^igi(Ci) +Pi(^i) + ^iPite) > &(^i) 


is the another form of the conditions (2.19). 


2.3.2 Analysis under No-Flux Boundary Conditions 

Similarly, we analyze the steady state solution under no-flux boundary conditions of our 
model (2.1), (2.3), (2.5), and (2.6). The steady state problem takes the form, 

rPn- 

Di-^ + Uigi{ui) - UiPi{ui) = 0 i = l,2 (2.23) 

The no-fltix boundary conditions are given by 

E«i(0) = 0 = -^«2(A) (2.24) 

as as 

And at the interface s = Li, we have 

(2.25) 

as as 
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Let gi{s, Pi), Li-i < s < Li, are the unique solutions of the equation (2.23), for i=l,2 
and qi{s,pi) are such that 

dqi 


ds 

dq2 

ds 


(0,A) = 0, gi(0,A) = /5i 

P2) = P2, Q2{L2, P2) — P2 


It may be noted here that it is possible to show the existence of a positive, monotonic, 
continuous with continuous flux solution, if we can show that there exist Pi , P 2 such that 

dqi{Li,Pi) _ dq2{Li,p2) , a\_„(T ^ n 

jJi 5 - — D2 , qi{Li, Pi) — q 2 [Li, P2) 


ds 


ds 


Multiplying both sides of (2.23) by Idui/ds and integrating from 0 if i=l and from 
L 2 if i=2, we get 

= ^ ^ in [gi(?i) - Pi(6)] d^i (2.26) 

A JPi 

By using (2.26), the following lemmas are proved analogous to their counterparts in the 
previous section. 


duj 

ds 


Lemma 2.3.4 If Pi > K*, then 


dpi{s,Pi) 

ds 


> 0 , 0 < s < Li 


Proof: Prom (2.26), we get. 


dpi{s,pi) 


ds 


rpi{s,Pi) 

I 6 gi( 6 ) -Pi( 6 ) 

Js-i 


Dll Hi 


(2.27) 


Since dpi{^,pi)/ds = 0, Pi(0,,5i) = Pi there exists si > 0 such that pi{s,Pi) > Pi on 
0 < s < Si- If not, let So, Q < sq < Li, he the least positive value , if it exists, such that 
Pi ( 50 , Pi) = Pi- Then by the mean value theorem there exists s such that 0 < s < sq and 
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dpi{s, Pi)/ds = 0 ; that is, 

2 rP\{s,h) 

0 = - 75 - / 6 [gi(ei) - Pi( 6 )] di, (2.28) 

L>n J0i 

But the right hand side of (2.28) is non-zero, since Pi{s,Pi) /?i on 0 < s < s, giving 
a contradiction. Hence pi{s,Pi) > f3i, for all value 0 < s < Li. Also /3i > Ky, implies, 
p(s, /?i) > Ki, on 0 < s < Li- Hence by ( 2 . 2 ) and (2.27), the result follows. I 


Lemma 2.3.5 If 0 < P 2 {s,p 2 ) < ^2 < K 2 u)hich implies p 2 {s,P 2 ) < P 2 , for 

L\ ^ S I'2- 


Proof: Since 0 < P2{s,/32) K 2 > from (2.2), we get, 

P2 [g 2 (P 2 ) - P 2 (P 2 )] > 0 

Now, since /?2 < A’l, then P 2 g 2 (P 2 ) - 92 P 2 (P 2 ) > 0 , V /32,P2- If possible let p 2 > P 2 , then 
the right hand side of (2.26) is negative, which is impossible. Therefore P 2 {s, h) < ^ 2 - B 


Lemma 2.3.6 Define by 'Hi(A) =Pi{Li,Pi), such that 


and 


-Hi : [KiJi] -4 [kik;], 
H 2 : [k,K;] ^ [kik;]. 


Proof: Since pi{Li,Pi) is continuous function, therefore if Hi{/3i) exists, then it must 
be continuous. We choose Hi(Kl) = this is possible since pi{0,Ki) = and 


2.3 Analysis of the Model 


34 


dpi/ds{0,Kl) = 0. Now jErom Lemma 2.3.4, if /?i > Kl and since ?^i(/3i) is 

continuous, therefore 'Hi{cg) = oo. Hence we can find to be the first value of pi 

such that 'HiiPi) = K^. 

Again, choose ?^ 2 (-^ 2 ) = P 2 (-L 2 , = ATj and from Lemma 2.3.5, if ^2 < K 2 , then 

P 2 (Li, ^ 2 ) < 1^2, and hence A{i(/? 2 ) < P 2 < K^. Then there exists P 2 such that 71202) = Ki- 
Hence the result. I 


Theorem 2.3.5 There exists a continuous, monotonic solution of system (2.23) with con- 
tinuous flux at Li . 

Proof: By Lemma 2.3.6, for each P 2 ^ p 2 K^, we can choose a such that AT* < /5i < 
for which pi(Li,/3i) = p 2 {Li,P 2 )- Hence (3i can be solved as a function of P 2 in the form 
Pi = fiPz) to give a continuous solution of (2.23). We now define, 

1(A) =£>.^(ii,/{A)) - D,^(U,h). (2.29) 

Then 


102) = -D2-£{LuP2) < 0 , TiK;) = Di-^{LuP,) > 0 . 

Now by continuity, there exists /? 2 , P 2 < P 2 < such that X(/? 2 ) = 0. 

Since f(p 2 ) = Pi > Lemma 2.3.4 shows that our solution is monotonic on 
0 < s < Li- By Lemma 2.3.5, we note that on Li < s < L 2 , the solution p 2 (s,/? 2 ) < P 2 - 
11 P 2 {s,P 2 ) is not monotonic over this interval, then there exists s, Li < s < L 2 , such that 
dp 2 {s, P 2 )/ds = 0. But then, by (2.26), ^ 3 ( 5 ) > 0, ^ 2 ( 5 ) > P 2 , giving a contradiction. This 
proves the theorem. i 


Finally, the stabilitj' conditions under no-flux boundary conditions are the same as in 
the case of reservoir boundary conditions and proofs are similar. 
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Theorem 2.3.6 The steady-state, continuous, monotonic solution of the system (2.1) with 
continuous flux matching condition at the interface (2.3) and under the no- flux boundary 
conditions (2.5), is locally asymptotically stable if 

d 

-r- (W2[g2(w2) - P 2 (« 2 )]) < 0, Kl <U2< (2.30) 

au2 


Theorem 2.3.7 The steady-state, continuous, monotonic solution of the system (2.1 ) with 
continuous flux matching condition at the interface (2.3), and under the no-flux boundary 
conditions (2.5), is non-linearly asymptotically stable in the subregion of < Ni, Ui < K^, 
for i = l,2, if, 


(NigflNi) - Uigflui)) (NiPiiNi) - Uipflui)) 


Ui 


Ni-Ui 


< 0 , 


(2.31) 


Remeirk: The effect of harvesting is to decrease the equilibrium level and modify the 
conditions for stability ( see [10] ). It is to be noted further that harvesting enhances 
stability, see conditions (2.30) and (2.31). 

2.3.3 The Case of Uniform Equilibrium State 

We show that uniform steady state of the system (2.1) -4 (2.6) is globally asymptotically 
stable in the homogeneous habitat, under both set of boundary conditions. In our model, 
there is a uniform steady state N{s,t) = iir,„ 0 < s < L 2 , t > 0, where K is the common 
carrying capacity in both the patches. 


Theorem 2.3.8 Let Kl = K 2 = K. Then the unique steady state solution to (2.7) and 
either (2.8) or (2.24), N{s,t) = K is globally asymptotically stable. 
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Proof: Let V (N) = V (t) be the positive definite function about N = K, given by 




Ni-K- Khx^ 
K . 


v{N) = v{t) = y / 

1=1 

Differentiating with respect to t, and using (2.1) we get, 

{NiSiiNt) - 


ds 


Ni^Kd^Ni 




ds 


Using (2.2), the first integral on right hand side becomes negative, for all values of Ni ^ K. 
Since under both set of boundary conditions. 


Hence V{N) < 0, and V{K) = 0. Therefore V{N) is negative definite over > 0 with 
respect to N = K. Hence the theorem. I 


Remark: On comparing the results of this section and the previous two sections, it can 
be concluded that the nonuniform steady state is not globally asymptotically stable due to 
patchiness of the habitat. 

In the following section, we will show numerically the existence of the unsteady state 
solutions as well as the positive, monotonic, continuous steady state of reaction-diffusion 
system of equations with reservoir boundary condition and flux matching condition at the 
interface. 


2.4 Numerical Examples 


In this section we analyze our model (2.1), (2.3), (2.4) and (2.6) for the particular growth 
rate functions 


1 


Ki. 


and Pi(iVi) = aiNi 


giiNi) = r.- 


(2.32) 
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numerically by using finite difference method. From (2.32) we get, 


where 


Oi(Ni) = gi(Ni) - pi(N,) = r, (l 


Aj , z — 1,2. 

Ti + aiKi 


KJ' 


(2.33) 


Without loss of generality let Ki < K 2 . Our aim is to solve both unsteady state and 
steady state solutions of the system. First we study the unsteady state solutions of the 
following problem, in a two-patch habitat, 


-(s, t) = riNi{s, t) 1 1 - I + t), i = 1, 2, 


with the initial and boundary conditions: 

Ni{s, 0) = Ki + and 

^2 

Wi(0,t)=Wi, N2{L2,t) = K2. 

The continuity and flux matching conditions at the interface s = Li are given by 


(2.34) 


(2.35) 

(2.36) 


N,{L„t) = N2{L2,t), Di^iLut) = D2^{L^,t). 


(2.37) 


The corresponding steady state solution (distribution) for this is given by the following 
equation in 0 < s < L 2 , 

with the reservoir boundary conditions (2.36) and continuity and the flux matching condi- 
tions at the interfaces (2.37). 

We numerically slove the system (2.32) to (2.38) for the unsteady and steady state 
solutions, with the stabihty conditions (2.13) and (2.19), in this case are given by, 

U2 > for Ki<U2< K2- 

Zi 


(2.39) 
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and for i = 1, 2, 

Ni + Ui> Ki, for Ki < Ni, Ui < K 2 - 

for linear and non-linear stability conditions respectively. These results 
various set of parameters in Fig.2.1 - Fig.2.5. 

Fig.2.1 and Fig.2.2 show that the unsteady state solutions are monoton 
to s and t. From Fig.2.1 and 2.2 we further note that as t becomes large, uns 
tends to steady state solution for the set of parameters satisfying stability con 
and (2.40). 

It is seen from Fig. 2. 3 - 2.5 that the steady state solution (distribution) it 
tonically increasing or decreasing from the starting point Ki to the final point j 
also shows the effect of diffusion coefficients on the steady state distribution, for i 
Dj, i = 1,2 (i.e. Di 0), the steady state distribution tends to the equilibri 
carrying capacity) level in the respective patches for all s except very close to the 
terface but for large A’s its behavior is linear. 
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Figure 2.1: The unsteady state and steady state solutions, for Di = 0.7, = 0.6, Ki = 

100, K 2 = 200, ri = 0.2 and ra = 0.3 
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Figure 2.2: The unsteady state and steady state solutions, for Di = 2.0, D 2 = 2.0, Ki = 
100, K 2 = 200, Ti = 0.2 and r 2 = 0.2 
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2.5 Summary 

A dynamical model of a single species population in a two patch habitats with harvesting 
and diffusion has been studied in this chapter. For application point of view, one may think 
of two adjoining forest stands whose carrying capacities are different due to deforestation, 
soil and other environmental conditions. 

By analyzing the model, it has been shown that there exists a unique positive, mono- 
tonic, continuous steady state solution under both reservoir and no-flux boundary con- 
ditions. The stability of this non-uniform steady state is discussed for both linear and 
nonlinear cases. It is shown that the corresponding uniform steady state is globally asymp- 
totically stable. Further, the comparison between the uniform and the non-uniform steady 
state solutions indicates that the patchiness decreases the stability of the system and the 
effect of harvesting is to decrease the level of equilibrium. These results are also proved by 
numerically solving the model under both steady and unsteady state conditions. 
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Chapter 3 


A Single Species Model with 
Diffusion and Supplementary 
Resource in Homogeneous and 
Two-patch Habitats 

3.1 Introduction 

Degradation of a forest habitat caused by industrialization, pollution and increasing pop- 
ulation, is of great concern to all humankind. A typical example in this regard is the 
depletion of the forest resources in the Doon Valley located in the foot hills of Himalayas, 
Utter Pradesh, India. Here the depletion of forest has been caused mainly by limestone 
quarrying, paper and other wood based industries and associated population migration 
[12]. Recently Shukla et al. [15] have proposed a mathematical model for forest depletion 
caused by resource independent industrialization (population) by considering the spatial 
distribution of both the forest biomass as well as the density of industrialization. By study- 
ing the behavior of uniform steady state solution, they have shown that if industrialization 
increases without control, the forest biomass may not last long. 
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This study implies that when forest stand or land is used for industrialization dis- 
tributed spatially, patchiness is caused in the forest habitat leading to resource deple- 
tion. Similarly, the associated population, distributed nonuniformly and dependent upon 
resource biomass, is also found to be causing patchiness. Therefore it is important to 
study the effect of patchiness on the growth and existence of diffusing populations. It 
is noted that little effort has been made to study such problems using mathematical 
models [2, 4, 5, 15], though there exist several studies related to population diffusion 
[1, 6, 7, 8, 9, 10, 11, 13, 14, 16]. Freedman et al. [4] studied a single species diffusion 
model by assuming that the habitat consists of two adjoining patches and shown that there 
exists a positive, monotonic, continuous non uniform steady state solution that is linear 
asymptotically stable under both reservoir and no-flux boundary conditions. Freedman 
and Kriszti [2] have also consider the same problem in three patch habitat and shown that 
the positive steady state solution is pice-wise monotonic. 

Further, when resource is depleted in the forest habitat due to industrialization etc., 
the species survival is threatened. One must therefore also study the effect of explicit 
dependence of population on resource in both homogeneous and patchy habitats. In this 
regard, Freedman and Shukla [3] have studied the effect of an alternative resource of preda- 
tor population on predator-prey systems by considering diffusion in a homogeneous habitat. 
They, however neither considered the effect of patchiness in the habitat nor the behavior 
of the non-uniform steady state solution of the system. 

In this chapter, we therefore, study the effect of supplementary self-renewable resource 
on a single species population with diffusion in a two patch habitat. The aim here is to 
show qualitatively that nonuniform steady state distribution of the species population is 
positive, continuous and monotonic through out the habitat and the level of distribution 
is higher then the case without supplementary resource [4]. We also obtain the conditions 
for stability of the system in both linear and nonhnear cases. 
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3.2 Mathematical Model 

We consider a linear forest habitat 0 < s < L^, consisting of two adjoining patches Lj-i < 
s < Li, for i = 1, 2, where Lo = 0 and Li is the interface of the two patch. Let Rt{s, t) is 
the density of the nondiffusing self-renewable supplementary resource and Ni(s,t), be the 
densities of the species population bilinearly dependent on Ri{s,t), at location s and time 
t in the i-th patch, for i = 1,2. It is assumed that Ri{s,t) grows logistically in both these 
patches with the different intrinsic growth rates (rj) and the same carrying capacity (C). 
Further we assume that the growth rates gi{Ni) of Ni{s,t), i = 1, 2 is general logistic type, 
in absence of the supplementary resource and is different in each patch. Keeping in view 
of the above, the model governing the system can be written as follows: 


- § 
II 

1 

1 

(3.1) 

dN- d^N- 

^ = NigiiNi) + BiPiRiNi + 

(3.2) 

A, ^i>0, 0 <s<L 2, f = l,2. 



where Pi and 6i are respectively depletion rates of resource and the conversion rates of 
resource biomass in the corresponding patches and Di are the dijffusion coefficient of Ni{s, t) 
in the i-th patch respectively, for i = 1, 2. 

The growth rate gi(iVi) functions are such that 
(H.l): gi{Ni) E C^[0, oo), gi(0) > 0, and m > 0, g'i{Ni) < 0. 


We note that when the habitat has carrying capacity Ki in the i-th patch respectively, 
then giiKi) = 0. 


Further we assume that 
(H.2): 3 > 0 such that 

gi{K*) + 9iPiR* = 0 and nil - - piK* = 0. 
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The existence of (R*, K*) are shown by graphically in Fig.3.1 an it is noted from this figure 
that the equilibrium level of population increases as the level of supplementary resource 
density i.e. its carrying capacity increases. Further it is also clear both analjdically and 
graphically, that < C and > Ki. 

The model (3.1), (3.2) is studied under both reservoir and no-flux boundary conditions 
separately. In reservoir boundary conditions, we assume, 

Ni{0,t) = K* and N 2 {L 2 ,t) == K; (3.3) 

and in the case of no-fiux boundary conditions, we have 

dN^{0,t) ^ dN2{L2,t) 

ds ds 

To study the behavior of non-uniform solutions, we also assume the continuity and 
flux matching conditions for Ni{s,t) at the interface s = Ti, V t > 0, as 

Ni{L^,t) = N2{L,,t\ = D2^(Li,t) (3.5) 

Since the resource Ri{s,t) is non-diffusing, these types of boundary and matching 
conditions are not required as there is no diffusion term in (3.1). It is further noted that at 
s = Li, Ri{s, t) will be continuous or discontinuous according as different levels of depletion 
rates in the two patches. 

Finally the model is completed by assuming some positive initial distribution for both 
forest resource biomass density and population density, that is. 


Riis, 0) = Xi(s) > 0, j 

Li-i <s<Li 

(3.6) 

Ni{s,0) = Si{s)>0, 1 

-'i— 1 <i S <. Li 

(3.7) 


(3.4) 
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Our main aim is to study the behavior of the non-uniform steady state solutions of 
the above model, in the patch 0 < s < L 2 , as this will through light on the effect of 
supplementary resource on the species survival in the habitat. 

Before we study our original model, in the following we state some important results 
for the corresponding model in a single homogeneous habitat 0 < s < L 2 without and with 
diffusion. 


3.3 Analysis of the Model in a Homogeneous Habitat 

3.3.1 Model without Diffusion 

and species population are 
Our model reduces to the 

(3.8) 

(3.9) 

From equation (3.8) and (3.9) it follows that there are four equilibria, namely Bq = 
[0,0], Be = [O', 0], Bk = [0,K] and B* = [i2*,iir*], where R*, K*, from (3.8) and (3.9) are 
given by, 

r(l - ^) - /3K* = 0, and g{K*) + epR* = 0. 

o 

For r > I3K, the existence of B* is shown in Fig.3.1. It is also noted that, R* < C and 
K* > K. The nature of R* and K* with respect to various parameters is shown in the 
table 3.1, which can be checked analytically. 


In this case both the supplementary forest resource biomass 
uniformly distributed for all 0 < s < 1^2 without diffusion, 
following form: 

dN 

V = Ng{N) + e^RN 
dt 
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Table 3.3: Effects of various parameters on R* and K*. 


(i) 

r t 

R*^ 

A"* t 

(ii) 

at 

R* t 

An 

(iii) 

ct 

R* t 

An 

(iv) 



An 

(v) 

9t 


An 

(Vi) 

Pt 

R* > C/2 

R* i 

An 


The notation f and I means increase and decrease respectively. 

It is also observed that the most interesting interior (non-zero) equilibrium E*, if it 
exists (i.e. r > PK), then by using the following positive definite functions, 

ViR,N) = d{R-R*f + {N-K*f, (3.10) 

V{N, R)=(^R-R* - +0(n-K*- (3.11) 

where d and d are positive constants, and by using Lyaponov direct method, we conclude 
that the equilibria E* is locally as well as globally asymptotically stable for suitable choice 
of d and 9. 

3.3.2 Model with Diffusion 

Now we study the behavior of the uniform steady state solution of the model (3.1) - (3.7) 
with diffusion in a homogeneous habitat. Here Ri{s,t), Ni{s,t), a,-, Ki, Di and pi become 
R{s,t), a, K D and P respectively, Vs € Then the model in this case can 

be written as, 

^ = tR{1 - - I3RN (3.12) 

^{s,t)=Ng(N)+epRN+D^{s,t) 


(3.13) 
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0 < ^ < 1, yS > 0, 0 < s < La- 

The boundary and initial conditions are respectively, 


N{0,t) = K* = N{L 2 ,t) (reservoir boundary conditions) (3-14) 


or 


and 


dN{L2,t) 

ds ds 


(no-flux boundary conditions) 


(3.15) 


R{s, 0) = x(s) >0, 0 < s < L 2 (3.16) 

N(s, 0) = Sis) >0, 0 <s<L2 (3.17) 

It can be easily verified that by using the same positive definite functions as given in (3.10) 
and (3.11), the following theorem is true. 


Theorem 3.3.1 The uniform steady state of the system, i.e. (R*, K*) is locally as well 
as globally asymptotically stable if r > ^K. 


Therefore if an equilibrium is stable in nondiflrusing case it is always stable with 
diffusion. Hence there is no diffusion instability of the system. 

Now we go for our original model, and study the behavior of the non-uniform steady 
state solutions of the system. 


3.4 Analysis of the Model with Diffusion in a Two- 
Patch Habitat 

3.4.1 The Non-Uniform Steady State under Reservoir Boundary 
Conditions 

We consider the steady state problem of the system (3.1) - (3.7). Our aim is to show that 
there exists a non-uniform positive, continuous and monotonic steady state distribution 
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of the species population (wt(s)) under reservoir boundary conditions and with continuous 
flux matching at the interface s = Li. We will note later that the steady state distribution 
of the supplementary forest resource biomass {wi{s), i = 1, 2) in each patch of the habitat is 
positive, continuous and monotonic but there may exist a continuity or jump discontinuity 
at the interface s = Li according as the relation between growth rates and interaction 
coefficients. The steady-state system takes the form, 


nwi ^1 - ^ j - piWiUi = 0 

(3.18) 

dP’Ui 

Di ^^2 + + diPiWiUi = 0 

(3.19) 

The reservoir boundary conditions are, 


= U2{L2,t)=Kl 

(3.20) 


The continuous and the matching conditions at the interface are, 

ux(Ti) = U2{Li), (3.21) 


From (3.18), we can solve lOj as a function of Ui, as follows 


Wi = C 




i = 1,2. 


(3.22) 


From (3.22) it is clear that if Ui are continuously monotonically increasing (decreas- 
ing) then Wi are continuously monotonically decreasing (increasing) with continuity or 
discontinuity at the interface s = Li, according as 


Pi P 2 Pi , P2 

— = or 

ri r2 ri T 2 


(3.23) 


Now by substituting the value of Wi from (3.22) in (3.19), we get 
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where Gi(ui) = gi(ui)-f and Gj(0) > 0, G-(ui) < 0, \/ui ^ 0. Therefore, 

for i = 1, 2, Gi{ui) are logistic functions. 

As in chapter 2, we can show that, there exists an unique steady state solution u of 
(3.24) satisfying, 

min {KIK^} < u{s) < max {K^K^} 

We now consider without loss of generality 0 < < u,- < z = 1,2. 

Let pi{s,ai), Li-i < s < Li, are the unique solution of the equation (3.24), for 
z = 1, 2, such that 

^(0,q;i) = tti, pi{0,ai)=Kl 
-^( 1 / 2 ) 0 ^ 2 ) = Oi2, P2{L2, OL 2 ) = K 2 

To show the existence of the monotonic solutions it is sufficient to show that there 
exists cti and Q !2 such that 

Pi(Li,ai) =p2{Li,a2), = D2^{L^,a2) 

Multiplying both side of (3.24) by 2dui/ds and integrating from 0 if z = 1 and from 
1/2 if z = 2, we get 

= “■ A Jk- (5.25) 

Hence by using equation (3.25) and proceeding in a similar manner as in the chapter 
2, the following Lemma and Theorem can be proved. 


dui 

ds 


Lemma 3.4.1 If ai > 0 , then 

dpi{s, ai) 


> ai on 0 < 


s < Li- 


ds 
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Lemma 3.4.2 IfO<p 2 <K 2 and 02 > 0, then 


dp2is,a2) 

ds 


> 0 : 2 , 


1j\ ^ ^ 2 - 


Lemma 3.4.3 Define Fi{ai) by Fi{ai) = pi{L\,OLi). Then there exists dj > 0 such that 

: [0, di] [K*, K*2] 

F 2 : [0,d2] [K^, Kl] 


Theorem 3.4.1 There exists a continuous, monotonic solution Ui of specie 
with the no-flux boundary conditions (3.20) and continuous flux at Li. 

Proof; From Lemmas 3.4.1 and 3.4.2, it follows that any solution we construct must be 
monotonic. By Lemma 3.4.3, for each 0 < 0:2 < <3:2, 'we can find an Oi such that 0 < Oi < di 
for which px{Li,a\) = ^ 2 (^ 1 , 02 )- Hence a\ can be solved as a function of a. 2 , Oi = h{o’ 2 ), 
to give a continuous solution of (3.24) with (3.20) and (3.21). 

Let 

\ dpi{Li,h{a2)) r, 5P2(Li,Q:2) 

G(a2) = D, gj D, 

Clearly 0 ( 0 : 2 ) is continuous on 0 < 0:2 < d 2 . Then we have 

G(0) > 0 , 

as 

and 


Hence the theorem. 
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Remsirk: Prom theorem 3.4.1 it is noted that if Kl > K 2 (or < i^l) then the steady 
state solution Ui is monotonically increasing (or decreasing), for 0 < s < L 2 and accordingly 
from (3.22) Wi is also continuously decreases (or increases) in both the patches, Q < s < Li 
and Li < s < L 2 . Further the steady state distribution Wi of the renewable supplementary 
resource is continuous and monotonic at each point of both the patches but continuous or 
discontinuous at the interface of the two adjoining patches according as (3.23). 

Now we find the conditions for asymptotic stability in both linear and non-linear case. 


Theorem 3.4.2 The steady-state continuous, monotonic solutions of system (3.1) and 
(3.2) with continuous solutions and flux at the interface s = Li is locally asymptotically 
stable if 

Bi = gi{ui) -f Uig[{ui) -f- Oi^iWi < 0 (3.26) 

and 

ai^[wi — OiUif < AAiBi (3.27) 

where 

A = ri (^1 - - PiUi = < 0, for i = l,2. 


Proof : Let the steady-state solution of system (3.18) be 

w{s) = 

And let the steady-state solution of system (3.19) be 


zi;i(s), 0 < s < Iri 
W 2 {s), Li < s <L 2 


u{s) = 


_ j ui{s), 0< s < Li 


U2{s), Li < S <L2 


Lineariazing (3.1) and (3.2) by using, 


Ri{s,t) = Wi{s) + mi{s,t) 


(3.28) 
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Ni{s, t) = -P ni(s, t) 


(3.29) 


We have, 


dvrij 

dt 


= nrii 


i (l 


TliPiWi 


dlXi T%i 

-^ = + ^iSiiUi) + OiPiWi] -f TTlieiPiUi -f Di-^ 


(3.30) 

(3.31) 


Using (3.28), the corresponding boundary and matching conditions can be obtained as 
follows 

ni(0,t) = K^, n 2 iL 2 ,t) = Kl 
ni(Li,t) =n 2 (Li,t), Di^(Lx,t) = D 2 ^(I^i, i) 

Now we consider the following positive definite function, 

A 1 


X •'Li-i ^ 


from which we get. 


T>/4^ V- I" 

By using (3.30), (3.31), (3.32) and (3.33), and integration by parts, we get. 


ds 


(3.32) 

(3.33) 

(3.34) 

(3.35) 


= IZ / mjAids-'^ f UimilPiWi - 9il3iUi]ds + "^2 [ nfBids 

I Li-i j J Li-i ]L 


2 ri, 


2 rLi 


where 


Ai = ri(l- ^ j - PiUi and Bi = gi{ui) + Uig-(ui) + OiPiWi 
Again, by using (3.22), for i = 1,2, we get 

jO \ 

TiWi 


^ 2riWi , Pi ^ 

A = -— + r, (!--«, 


< 0 


and for i = 1, 


Bi = gi(«i) + uig'x(Mi) + diPiC ^1 — 


(3.36) 


(3.37) 


(3.38) 




3.4 Analysis of the Model with DiflFusion in a Two-Patch Habitat 


60 


since ui > > Ki, therefore gi(tti) < 0 and gi(ui) < 0, V ui. Hence the condition (3.38) 

is a very weak condition, later we will show in an example that for simple logistic function 
= tii(l — uifKi), the condition is automatically satisfied. 

From (3.36) it is noted that V{t) is negative definite, if the conditions (3.26) and 
(3.27) axe satisfied, for i = 1,2. 

In a similar manner, the following nonlinear stability result can be proved by using 
the same positive definite function as in theorem 3.4.2. 

We consider the case, when Kl < then R\>B^. 


Theorem 3.4.3 The steady-state continuous, monotonic solutions of the nonlinear system 
(3.1) and (3.2) with continuous flux at the interface s = L\ is asymptotically stable in the 
subregion o/R; < Ri,Wi< < Ni,Ui < provided the following conditions 

are satisfied: 


(i) 


Ni-Ui 


(w) 


•^2i = n 



Ri + Wi\ 

C ) 


- 0iNi < 0 


(Hi) 


01 [wi - 6iUif < ATiiT^2i 


3.4.2 The Non-Uniform Steady State under No-Flux Boundary 
Conditions 

In this section, consider the steady state problem of the system (3.18) and (3.19) with 
no-fiux boundary condition. 


(3.39) 
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and continuity and flux matching conditions at the interface (3.21). 


In this case also, there exists a non-uniform positive, continuous and monotonic steady 
state distribution of the species population density (ui(s)) and the steady state distribu- 
tion of the supplementary resource biomass {wi{s), i = 1,2) in each patch of the habitat 
is positive, continuous and monotonic but it may be continuous or discontinuous at the 
interface s — Li, because of the following relation. 


Wi = C 



i = l,2. 


(3.40) 


and the condition (3.23). 


In this case also, it is clear from the above that if Uj are continuously monotonically 
increasing (decreasing) then Wi are continuously monotonically decreasing (increasing) with 
continuity or discontinuity at the interface s = Li, according as (3.23). 

After substituting the value of Wi from (3.40) in (3.19), we get 

* 

Di-^ + UiGi{ui) = 0, 2 = 1, 2 (3.41) 

where Gt(ui) = gj(ui) -1- 9iPiC{l - PiUi/ri). 


Let Pi(s, Pi), Li-i < s < Li, are the unique solutions of (3.41), for i = 1,2, such that 

=0 . Pi(0,A)=A (3.42) 

^^^^ = 0 , P2(i2,A) = A (3.43) 

The existence of the monotonic solutions follows, if we can show that there exist Pi and 
P 2 , such that 

dpi{Li,Pi) _ ^ dp2{Li,P2) 


PiiLi,Pi) =P2(Li,p2) , Di 


(3.44) 
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From (3.41), multiplying both side by 2dui/ds and then integrating with respect to s from 
0 if i = 1 and from L 2 if i = 2 , we get 

‘'Pi 

Then by exactly in similar manner as in previous section, we can state the following 
theorem. 

Theorem 3.4.4 There exists a continuous, monotonic solution Ui of species population 
with the no-flux boundary conditions (S.S9) and continuous flux at Li. 

We now state the linear and non-linear asymptotic stability theorems of the non- 
uniform steady state solutions of the system (3.18), (3.19) with (3.39) and (3.21) by using 
the Lyaponov direct method. 


duj 

ds 


Theorem 3.4.5 The steady-state continuous, monotonic solutions of system (8.1) and 
(3.2) with continuous solutions and flux at the interface s = Li is locally asymptotically 


stable if 

Bi = gi(ui) -f Uig((ui) -1- 9ipiWi < 0 

(3.46) 

and 

ai^[wi - OiUi]"^ < 4AiBi 

(3.47) 

where 




Ai = ri(l- - PiUi = < 0, for i = l,2. 


Proof: Same as Theorem 3.4.2 and only difference is in the condition (3.32), which becomes 


dni 

ds 


(0,t) = 0 = 


dn2 

ds 




now. 


(3.48) 
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In a similar manner, the following nonlinear stability result can be proved by using 
positive definite function as in theorem 3.4.3, for Kl < K^. 


Theorem 3.4.6 The steady state continuous, monotonic solutions of nonlinear system 
with continuous flux at the interface s = Li is asymptotically stable in the subregion of R; 
{R^ < Ri, Wi < R^] Ni < Ni, Ui < provided the following conditions are satisfied: 


(i) 


Tu = — < 0 

- Ui 


(ii) 


^2i = fi 



Ri Wi\ 

C ) 


-PiNi<0 


(ui) 


/?? [u)i - OiUif < 4J’ii J‘2i 


The above theorems imply that the system will settle down to the steady state dis- 
tribution in the two patches under certain conditions, the magnitude of the steady state 
distribution of resources biomass density being lower than it’s original value but the density 
distribution of population being correspondingly higher in each patch. 

3.4.3 The Uniform Steady State under Both Sets of Boundary 
Conditions 

Now we consider the case when the species population is uniformly distributed between 
[0,L2]) be- when Ni{s,t) = N 2 {s,t) = I*, s E [0,1^] and V t > 0. Since the depletion 
of supplementary resource {Ri{s,t)) due to population is different in different patch, i.e. 
/?i 7 ^ ^ 2 - Let Ri{s,t) = RJ, 0 < s < Li, V f > 0 and R 2 {s,t) = i?;, Li < s < L 2 , V t > 0. 
Without loss of generality, we choose ^ 2 > Pi, this implies from (3.22) that RI > R^. Now 
we go for the globally stability of the system. 
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Theorem 3.4.7 If ui = U 2 = I* , Wi = and W 2 = R^i Then the system is globally 
asymptotically stable. 


Proof: Let V{x,y) be the positive definite function about Hi = R*, Ni = I*, given by 

V{x, y) = E (Ri-H*- ds + jz {^i -I* - ds (3.49) 

Difi'erentiating (3.49) with respect to t, and using (3.1) and (3.2) we get, 

„ r rdi 1 r^i 


ns,t) = E r (Ri-Rt? -y ds + Y: f ' {Ni-rf 

jZD- r 


'ka 


ds 


+ 


Ni ds'^ 

Since under both set of boundary conditions, 


ds 


[x^{o,t) - = ML2,t) - Rn^iR2,t) = 0 


Then by using (3.51), we get 

2 Ni - r d^Ni 


(3.50) 


(3.51) 


z'" J'* (my, 


Hence from (3.50) and (3.51), F < 0 VRi ^ R^, \fNi ^ /*, and V{Rl, I*) = 0. 
Therefore V{Ri, N) is negative definite, proving the theorem. I 

It may be remarked from the above nonlinear steady state analysis for non-uniform 
and uniform cases that patchiness destablizing the system, we can see by comparing the 
statement of the theorem 3.4.7 with theorem 3.4.2, and theorem 3.4.3 or theorem 3.4.5 and 
theorem 3.4.6. 
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3.5 A Particular Case 


For clarity and simplification, we consider an example by takeing gj(iV’i) = Oj (1 — Ni/Ki). 
Then the corresponding system can be written as, 


BNj 

dt 


d'^N- 

aiNi (1 - Ni/Ki) + OipiRiNi + Di-^ 


the existence of non-negative equilibrium {R*,K*) is given by. 


R* = 


aiC{ri - l3iKi) ^ riKijoj -f pjOjC) 

rai + l5'/eiCKi ’ ' “ nai -h ^feiCKi ’ 


provided Vi > PiKi. 


(3.53) 

(3.54) 


From above it is clear that R* < C and K* > Ki. This shows that the effect of 
supplementary resource is to increase the equilibrium level of the species population. Also 
we note that in absence of interaction (i.e. Pi = 0)^ R\ = R^ = C , the biomass is uniformly 
distributed in the entire habitat. 


The model is completed by taking same boundary, matching and initial conditions as 
in the general case (3.3) — )■ (3.7). 

Exactly in a similar as in the general case we study the above system (3.53) and (3.54) 
with (3.3) — >• (3.7) in both homogeneous habitat (0 < s < L 2 ) without and with diffusion 
and in two patch habitat. 


3.5.1 Analysis of the Model in a Homogeneous Habitat 

Case 1: Model without diffusion: 

In this case both the supplementary forest resource biomass and species population 
are uniformly distributed for all 0 < s < L 2 . Our model reduces to the following form: 

r tR(1 -^}-/3RN 


(3.55) 
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dN N 


(3.56) 


From equation (3.55) and (3.56) it follows that there are four equilibria, namely 
Eo = [0, Qi],Ec = [0, C\,Ek = [K, 0] and E* = [R*,K% where 


R* 


aC{T - PK) ,,, rK(a + /39C) 

ra + /3^eCK ’ ro + I5^9CK ’ ’’ > 


(3.57) 


It is also noted from (3.57) that, R* < C and K* > K. 


The variational matrix of the system is given by 
r (1 - 2R/C) - I3N 


[M] = 


9I3N 


-/3R 

a (1 - 2N/K) - epR 


(3.58) 


From (3.58) and the standard stability theory, we note the following obvious observa- 
tions. The equilibria Eq is always unstable. The equilibrium points Ec and Ek are stable 
or saddle point according as o -f- OpC and r — /3K are negative or positive respectively. It 
is also observed that the most interesting interior (non-zero) equilibria E* is locally stable, 
if it exists (i.e. r > ^K). 


Alternatively if r > PK, by using the same positive definite functions as in (3.10) and 
(3.11), and by using Lyaponov direct method, we conclude that the equilibria E* is locally 
as well as globally asymptotically stable. 


Case 2: Model with dififusion: 

Now we study the behavior of the uniform steady state solution of the model (3.53) 
- (3.54) with diflrusion in a single homogeneous habitat. In this case we consider that both 
the populations are spatially distributed. Here Ri{s,t), Ni{s,t), o^, Ki, Di and /?,■ becomes 
il(s, t), A’(s,t), a, K D and P respectively, Vs € [0,1/2]. Then the model in this case can 
be written as. 


f-rR(l-§)-^RN 


(3.59) 
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ON N d'^N 

-^{s.t) = oJV(I - -) + ej3RN + D—{3, t) 

o<«<i. /3>o, a<s<U- 

The boundary and initial conditions are respectively, 


N{Q,t) = K* = N{L 2 ,t), or 


aiV(0,t) 

ds 


0 = 


dN{L2,t) 

ds 


(3.60) 


(3.61) 


and 


R{s, 0) = x(s) >0, 0 < s < L 2 , and iV(s, 0) = d(s) > 0, 0 < s < L 2 (3.62) 

It can be easily verified that under same positive definite functions given in (3.10) and 
(3.11), the following theorem is true. 


Theorem 3.5.1 The uniform steady state of the system, i.e. (R*, K* ) is locally as well 
as globally asymptotically stable if r > PK. 


Therefore if an equilibrium is stable in nondifPusing case it is always stable with 
difiiision. 

Now we go for our original model (3.53), (3.54), and study the behavior of the non- 
uniform steady state solutions of the system. 


3.5.2 Analysis of the Model with Diffusion in a Two-Patch Habi- 
tat 


Case 1: Analysis imder Reservoir Boundary Conditions 


Consider the following steady state system. 



— PiWiUi = 0 


(3.63) 
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(fui , Ui 




(3.64) 


with (3.20) and (3.21). 


Again from (3.63), we can solve Wi as a. function of Ui, as follows 

A 


= (7 1 1 -Uij , i = 1, 2. 


(3.65) 


From (3.65) it is clear that if Ui are continuously monotonically increasing (decreas- 
ing) then Wi are continuously monotonically decreasing (increasing) with continuity or 
discontinuity at the interface s = Li, according as 


Pi P2 Pi , P2 

— = — , or — ^ . 

n r2 n r2 


(3.66) 


Now by substituting the value of Wi from (3.65) in (3.64), we get 

P^Ui _ Ui^ 




Kt 

I > 


0, i = l,2 


(3.67) 


where Ci = ai + OiPiC. 


Exactly as same in general case (chapter 2), there exists an unique steady state 
solution u of (3.67), which satisfies 

min {KI,K;} < u{s) < max {KIK;} 

We now consider without loss of generality 0 < < Ui < K^- 

Let pi{s, ap, Li-i < s < Li, are the unique solution of the equation (3.67), for i=l,2. 
And pi{s, ai) are such that 

^(0, ai) = tti, Pi(0,q;i) = AT* 

^{L2,a2) = a2, p2{L2,a2)=K; 

as 
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Multiplying both side of (3.67) by 2dui/ds and integrating from 0 if i=l and from L2 
if i=2, we get 



Hence by using (3.68), the following theorem holds. 


Theorem 3.5.2 There exists a continuous, monotonic solution u, of species population 
with the no-flux boundary conditions (3.20) and continuous flux at Li. 


Now we find the most interesting results, the conditions for asymptotic stability in 
both linear and non-linear cases. 


Theorem 3.5.3 The steady-state continuous, monotonic solutions 
(8.64) with continuous solutions and flux at the interface s = Li is 
stable if 

K* < min. 


K: 


2 i 


Fi -\- (Ff - AEiGiyF 

2Ei 


of system (3.63) and 
locally asymptotically 

(3.69) 


and 


max. 



r2K2(a2 + 02020) F 2 - (Fj - 4E2G2)^^^ 
2o;2r2 + 0202^^2 2 E 2 


(3.70) 


Proof : Let the steady state solution of system (3.63) be 

, , f u>i(s), 0 ^ s ^ Li 
= I L,<s<L, 

Also let the steady-state solution of system (3.64) be 

0 < s < Li 

~ \ U2(S), Li<S< L 2 
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Lineariazing (3.1) and (3.2) by using, 

Ri{s, t) = Wi{s) + Tni{s, t) 

A/" t (s, t) — Uj(s) “b Tij (a, i) 
dirii f / 2w., 

K' 

2ui 


We have, 


at 


TJT-i r (^1 - - PiUi 


- ni/3iWi 


drii 

dt 


= rii 


aiil- ~]+ 9i/3iWi 

Ki 


O^Tl' 

+ rmeijSiUi + A-^ 


(3.71) 

(3.72) 

(3.73) 

(3.74) 


Using (3.71) and (3.72), the corresponding boundary and matching conditions can be ob- 
tained as followsO 

= n2{L2,t) = K*^ 


Now we consider the following positive definite function. 


V(t)=t£_^llml + ,^]ds. 


from which we get, 


[Li 

drrii 

dni 

JLi-i 

rat 

•"'“arj 


ds 


(3.75) 

(3.76) 

(3.77) 

(3.78) 


By using (3.73), (3.74), (3.75) and (3.76), and integration by parts, we get. 


^i't) — i rnlAids — ^ - Oi^iU^ds + ^ / nfBids 

^ J Li 1 Li 1 J Li — 1 


2 rU 


2 rl'i 


2 .Li 

JLi^l 


drij 

ds 


ds 


(3.79) 


where 


Ai = ri(l- - Pm and Bi = aj ^1 - + OipiWi 

Agjiin, by using (3.65), for i = 1,2, we get 
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B, = (Qi + 0iAC')- 


Q^i + OiPiC ^ Oii 

Wr 


Ui 


(3.81) 


Let 


A/f‘ ry * 

Mi = ai + 9ipiC, A/i = -jM -|- -~ 


Now from (3.81), we note that for Bi < 0, we have 


Mi riKi{ai + OifSiC) . , ^ 

Ui> — = r -, i = 1, 2 


(3.82) 


Mi 2airi + di^fCKi 
Since < Ui < K^, it is clear from (3.82) that ui > Mi/Mi is automatically satisfied. 
Therefore this condition is required only for U 2 , for Bs < 0. 

From (3.79) it is noted that 'F(t) is negative definite, if (3.82) is satisfied along with 
the following condition, for i = 1, 2 


0^[wi - diU^^ < 4ABi 

Now from (3.83), by using (3.65) and solving for Uj, i = 1,2, we get 


Kl < 


Fi - {Ff - AEiGif!^ 


2Ei 


<Ui< 


Fi + (if _ AEiGif!^ 

2Ei 


<k; 


where 


(3.83) 


(3.84) 


Ei = /3f +4ftVi. Fi = 2ft C+4lriMi+/3iKi), Gt = fffCP+iuKi. 

Since < Ui < K^, and by using (3.82), (3.84), we get the conditions (3.69) and (3.70). 
Hence the theorem. i 


The following numerical example shows the feasibility of the conditions (3.69) and 
(3.70). 

Example: Let 


r = 0.2, ai = 0.2, a2 = 0.25, 
(7= 1000, Ki=40, K2 = 70, 


/?i = 0.002, /32 = 0.0025 
$1 = 0.2, 02 = 0.5 
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we get, 

= 66.6667, = 78.1395 

= 77295.45, 
= 0.1146. 

Hence in this case, (3.69) and (3.70) becomes 


- (Ff _ 4E^Giy/y2Ei = 0.1098, + {F^ - 4EiGiy/y2Ei 

F2 + (i^l - ^E^Gyyi'^ I 2 E 2 = 50648.8, F 2 - (F| - AE2G2yf^l2E2 


Kl = 66.6667 <Ui< K; = 78.1395, i = 1, 2. 

In a similar manner, the following nonlinear stability result can be proved by using 
positive definite function as in theorem 2.3, for Ki < K^. 


Theorem 3.5.4 The steady state continuous, monotonic solutions of nonlinear system 
with continuous flux at the interface s = Li is asymptotically stable in the subregion o/R; 
{R^ < Ri,Wi < Ry, Ni < Ni,Ui < N 2 }, provided the following conditions are satisfied: 

(i) Fii = az (1 - + ^iPiRi < 0 

(u) F2i = ri (1 - m < 0 

{in) Pf [wt - diUif < 4TxiT2i 


Proof: By using (3.72) and (3.71) the nonlinear model (3.1) and (3.2), we get 


dui 

~di 


= Hi 


■t\i 


I riiiOiP^Ui ”1" 


d^rii 


(3.85) 


dmi 


= m,- 


r,- (1 - - m] - riiPm 


dt 


(3.86) 
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We consider the following positive definite function, 

A 1 




from which we get, 

2 rU 


Ni + Ui 


v(t) = E r "? h (i - I + 


Ki 


ds + 


^ rLi 

YL I Pi (^i ~ 

x=l 


+ 


tr 


mj 


Ti 1 


C 


‘ n.. 

i=l 


ds^ ’ 


Hence V{t) is negative definite if the conditions (i), (ii) and (iii) holds. 
Note: The conditions (i) and (ii) are always satisfied if, 


ai < min 


d" a a -D 

oii — QiPiRi 


Ti < min 


Ti- 


Ki 
Ri + Wi 


1< 


Ki ai '• 


C 


^PiNi 


C Ti 


Now using (3.89) and (3.90), we get the condition (iii) is always satisfied if 


(3.87) 


(3.88) 


(3.89) 

(3.90) 


+ Bi^Kf - 2diR^Ki\ < ^QT^iQT^i 

[l 

Ki ai 

1 

1 

(M 

1 

1 — { 

1 



(3 


,91) 


Hence the system is nonlinearly asymptotically stable if (3.89) — t (3.91) are satisfied. 


Case 2: Analysis under No-Flux Boundary Conditions 

The same analysis and results as in reservoir case (i.e. Case 1), are true for no-flux 
boundary conditions also. Therefore the nonuniform steady state is positive, monotonic 
and the system is linearly as well as nonlinearly stable under same set of condition as in 
theorem 3.5.2, 3.5.3 and 3.5.4. 
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3.5.3 The Uniform Steady State under Both Sets of Boundary 
Conditions 

Now we consider the case when the species population is uniformly distributed between 
[0,1/2], i.e. when Ni{s,t) = N 2 {s,t) = 7*, V s € [0,1-2] and V t > 0. Since the depletion 
of supplementary resource {Ri{s,t)) due to population is different in different patch, i.e. 
Pi 7 ^ P 2 , let Ri{s,t) = i?*, 0 < s < Li.'i t > 0 and R 2 {s,t) = R^, Li < s < L 2 , V t > 0. 
Without loss of generality, we choose Po > Pi, this implies from (3.22) that Rl > R^. Now 
we go for the globally stability of the sj'stem. 


Theorem 3.5.5 If ui — U 2 = I*, Wi = R^ and W 2 = R 2 , Then the system is globally 
asymptotically stable. 


Proof: Let V{x,y) be the positive definite function about Ri = R*, Ni = I*, given by 
V{x, y) = j^ei (Ri - R* - j ds + j^ ^ [Ni -I* - ds (3.92) 

Differentiating (3.92) with respect to t. and using (3.1) and (3.2) we get. 


^ rLi ^ rLi 

V{s,t) = + (AT.--/*) 


'Ki 


ds 




Nj - r d^Nj 

ds^ 


ds 


Since under both sets of boundary conditions, 


[xi{0,t) - = [x 2 {L 2 ,t) - K^]-^{L 2 ,t) = 0 


Then by using (3.94), we get 

2 rU Ni - r ePNi 


(3.93) 


(3.94) 


Ni-r^N. 2 Li j* [dNiV, 

Hence from (3.93) and (3.94), T> < 0 V7?i Rl ^Ni / 7*, and V{Rl I*) = 0. 
Therefore V{Ri,Ni) is negative definite, proving the theorem. ■ 
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3.6 Numerical Examples 


In this section we calculate and study, the numerical solution of the steady state solutions 
the above particular case (3.63) and (3.64) with reservoir boundary conditions, and compare 
the results with the case discussed in chapter 1, for the following three set of values of the 
parameters, 


Parameters 

Figure 3.2 

Figure 3.3 

Figure 3.4 

Dx 

0.6 

0.6 

0.6 

D2 

0.7 

0.7 

0.7 

ai 

0.08 

0.08 

0.08 

(l2 

0.08 

0.08 

0.08 

Ki 

100 

50 

50 

K2 

50 

100 

100 

n 

0.4 

0.4 

0.4 

r2 

0.8 

0.8 

0.6 

C 

40 

40 

40 

A 

0.0003 

0.0003 

0.0003 

/32 

0.0006 

0.0006 

0.0006 

e 

0.5 

0.5 

0.5 


Using these values of parameters, we get 


Parameters 

Figure 3.2 

Figure 3.3 

Figure 3.4 

Kl 

106.9 

53.6 

53.6 


57.18 

113.72 

113.3 

R\ 

36.79 

38.39 

38.39 


38.28 

36.58 

35.46 


The results are shown in Fig.3.2, 3.3, 3.4. It can be noted that, in presence of the 
supplementary resource the level of the steady state distributions are higher at each point 
of the habitat. Again, in Fig.3.4, we consider the case when A/^'i 7^ it is shown that 

the steady state distribution of the renewable supplementary resource has a discontinuity 
at the interface of the patches. But in the first two figures, since = yd2/7’2, hence, the 
resource is continuous at the interface. 
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3.7 Summary 

In this chapter, a mathematical model is proposed to study the role of supplementary 
renewable resource on a single species population model in a two patch habitat In the model 
it is assumed that the density of resource biomass is governed by the logistic equation with 
the different intrinsic growth rate but the same carrying capacity in the entire habitat. It 
is further assumed that the densities of species population is also governed by the logistic 
equations in both the patches but with different growth rates and carrying capacities. The 
rate of depletion of resource hiomass density due to population is also considered to be 
different in the two patches. 

It is shown that there exists a positive, monotonically increasing, continuous steady- 
state solution with continuous flux, in the case of both reservoir and no-flux boundary 
conditions, for both forest resource biomass and species population, that is asymptotically 
stable in the both linear and non linear cases under some conditions. It has been shown that 
forest resources remains at steady-state but at a lower level then it’s carrying capacity the 
magnitude of which depends upon the level of population, but the marnitude of the steady 
state distribution of the population is always higher in presence of resource in comparison 
to the case of without resource. It is also noted that the patchiness in the habitat has a 
destablizing effect on the system. 
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4.2 Mathematical Model 

We consider a dynamic model of two logistically growing species with prey-predator type 
interaction and diffusion in a two-patch forest habitat by assuming that the second species 
(i.e. predator population) uses the first species (prey population) as an alternative resource. 
In such a case the rate of change of density of the first species decreases due increase in 
the density of the second species, but the density of the second species increases due to the 
increase in the density of the first species in both the patches. Let Xi{s,t) and yi{s,t) be 
the densities of first and second species in the i-th patch respectively. Then the model can 
be written as the following system of autonomous partial differential equations: 

= xds, t)gi t)) - a(,, t)p( t)) + Du (4.1) 

t)f, (y,(s,i)) + ()Pi t)) + (4.2) 

0 < s < L 2 , z = 1 , 2 , 7 t > 0 

where the i-th patch is assumed to lie along the spatial length Li_i < s < Li (Lq = 0 ), 
ii(xi) and fi(yi) are the respective specific growth rates, Pt( 2 :i) is the interaction rates, and 
Dii and I? 2 i are the diffusion coefficients of Xi and yi in the two patches, 7 ^ > 0 be the 
jrowth rate coefficient of yi due to Xi in the i-th patch respectively. 

In writing the system (4.1) and (4.2) it is assumed that the growth rates and carrying 
apacities of two species are smaller or larger in the first patch than their corresponding 
alues in the second patch depending upon whether the first patch is more or less depleted, 
^’hus we assume gi{xi) and fi(yi) are logistic type functions and satisfying Hi: 

Ii: gi{xi),fi{yi)eC^[Q,<^) 

gi(0) > 0, for Xi > 0, E,'iixi) < 0, 
fi( 0 ) > 0 , for yi > 0 , f/(yi) < 0 
gi(JTi) = 0 , fi(Mi)= 0 , z = l, 2 . 
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3.7 Summary 

In this chapter, a mathematical model is proposed to study the role of supplementary 
renewable resource on a single species population model in a two patch habitat In the model 
it is assumed that the density of resource biomass is governed by the logistic equation with 
the different intrinsic growth rate but the same carr 3 dng capacity in the entire habitat. It 
is further assumed that the densities of species population is also governed by the logistic 
equations in both the patches but with different growth rates and carrying capacities. The 
rate of depletion of resource biomass density due to population is also considered to be 
different in the two patches. 

It is shown that there exists a positive, monotonically increasing, continuous steady- 
state solution with continuous flux, in the case of both reservoir and no-flux boundary 
conditions, for both forest resource biomass and species population, that is asymptotically 
stable in the both linear and non linear cases under some conditions. It has been shown that 
forest resources remains at steady-state but at a lower level then it’s carrying capacity the 
magnitude of which depends upon the level of population, but the marnitude of the steady 
state distribution of the population is always higher in presence of resource in comparison 
to the case of without resource. It is also noted that the patchiness in the habitat has a 
destablizing effect on the system. 
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Chapter 4 


A Prey-Predator Type Model with 
Diffusion in Homogeneous and 
Two-Patch Habitats 


4.1 Introduction 

An interesting problem in mathematical ecology is to study the growth and co-existence of 
species with diffusion in both homogeneous and heterogeneous habitats. A diffusion process 
in the habitat tends to give rise to a uniform density of population in space. As a conse- 
quence, it is expected that diffusion, when it occurs, plays the role of increasing stability 
in a system of interacting populations [1,4 - 4 - 10, 13 -4 26,28,29 and their cross references]. 
Levin [14, 15] has given elaborate survey of models with diffusion in both homogeneous 
and heterogeneous environment. McMurtrie [16], has also surveyed the literature related 
to populations model with diffusion and reported the effects of dispersal and spatial het- 
erogeneity on stability of both single species and for Predator-Prey system. Hastings [9], 
has studied the global stability in n-species Lotka-Volterra systems with diffusion. He also 
established a necessary and sufficient conditions for stability of a single species model in the 
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case of large diffusion rates [9]. Nallaswamy and Shukla [18], considered a prey-predator 
model with functional response and diffusion and have shown, that if the equilibrium state 
is linearly stable, a subregion of the positive quadrant can be found in the phase plane 
where it is non-linearly stable with or without diffusion. Freedman and Shukla [4], have 
studied the effect of a predator resource on a diffusive Predator-Prey system, showing the 
stabilizing role of diffusion. There is an important exception, known as ’’diffusive insta- 
bility” [19, 20, 24, 31, 32], however which might not be a rare event especially in aquatic 
systems. 

In most of the above mentioned studies the existence of the uniform steady state 
has been assumed and the corresponding model is analyzed in some particular cases. In 
a few cases the stability of non-uniform steady state has also been studied in a homoge- 
neous environment [33]. It may, however, be noted that in real ecosystems the habitats 
are often patchy due to resource depletion or environmental changes and populations are 
usually distributed heterogeneously in their habitats. It may be pointed out that here that 
little effort has been made to study such systems with diffusion in a patchy habitat using 
mathematical models. Freedman et al. [5] and Freedman and Wu [6], have studied a single 
species diffusion model by assuming that the habitat consists of respectively two and three 
adjoining patches. In [5] Freedman et al. have shown that there exists a positive, mono- 
tonic, continuous non-uniform steady state solution that is locally asymptotically stable 
under both reservoir and no-flux boundary conditions. They, however, did not study the 
stability of non-linear system, nor the case of two species interaction in a patchy habitat. 

In this chapter, we, therefore, study a logistically growing two species prey-predator 
type model with diffusion in a homogeneous (non-patchy) and two-patch habitats and 
discuss the stability of both the linear and non-linear systems, under both the reservoir 
and no-flux boundary conditions. The model is proposed by keeping in view the depletion of 
forest resources by industrialization and population causing patchiness in the Doon Valley 
situated at the foot hills of Himalayas in India, [25]. 
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4.2 Mathematical Model 

We consider a dynamic model of two logistically growing species with prey-predator type 
interaction and diffusion in a two-patch forest habitat by assuming that the second species 
(i.e. predator population) uses the first species (prey population) as an alternative resource. 
In such a case the rate of change of density of the first species decreases due increase in 
the density of the second species, but the density of the second species increases due to the 
increase in the density of the first species in both the patches. Let Xi{s,t) and yi{s,t) be 
the densities of first and second species in the i-th patch respectively. Then the model can 
be written as the following system of autonomous partial differential equations: 

= Xi{s,t)gi ixi{s,t)) - yi{s,t)j)i {xi{s,t)) d- (4.1) 

{yi{s,t)) + 7iyi(s,i)Pi (a:i(s,t)) + (4.2) 

0 < s < A, i = 1, 2, 7i > 0 

where the i-th patch is assumed to lie along the spatial length Li_i < s < Li (Lq = 0), 
gj(a:j) and fi{yi) are the respective specific growth rates, Pi{xi) is the interaction rates, and 
Du and D2i are the diffusion coefficients of Xi and yi in the two patches, 7, > 0 be the 
growth rate coefficient of yi due to Xi in the i-th patch respectively. 

In writing the system (4.1) and (4.2) it is assumed that the growth rates and carrying 
capacities of two species are smaller or larger in the first patch than their corresponding 
values in the second patch depending upon whether the first patch is more or less depleted. 
Thus we assume gi{xi) and fi{yi) are logistic type functions and satisfying Hi: 

Hi : gi(2;i), fi(yi) € 0 ^( 0 , 00) 

, gi(0) > 0, for Xi > 0, g[ixi) < 0, 

fi(0) > 0, for yi > 0, f/(yi) < 0 
gi(A) = 0, fi(Mi)=0, z = l,2. 
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Where Ki, Mi are the carrying capacity of the first and second species respectively 
in the i-th patch. 

In (4.1) the functional response (predator response function) Pi(a:i) satisfies the fol- 
lowing conditions, 

H2 : Pi(2;i) € C 2 [ 0 , 00) 

Pi(0) = 0, for Xi > 0, p((2:i) >0, i = 1, 2. 

Several types of gi(n;i), f,(yi) and Pi(a;i) have been catalogued in [10]. For examples, 
gi(2:i) = Tifl - {xi/Ki)], gi(xi) = ri{Ki - Xi)/{Ki -j- CiX,-) or gi(xi) = ri[l - {xi/Ki)% 
1 > c > 0. Similar forms are also valid for ft(yi). The predator response function, Pi(xi) 

is assumed to be bounded above and examples for which are, Pi(xi) = aiXi/{ai +Xj), 

Pi(xi) = ctjXj or Pi(xi) = Q!i(l - e~“*'^‘). 

We further assume that, 

H3: 3 X,* > 0, y* > 0 such that 

gt (^i ) - Vi Pi ) = 0 (4.3) 

y* [ft iy* ) + 7iPi (a;,* )] = o (4.4) 

From (4.3) and H2, we have gt(2;*) > 0. From Hi, we note that giiKi) = 0 and 
g'(xj) <0, V Xj > 0. This implies that gi(x*) > 0 = gi{Ki). Hence x^ < Ki. Also from 
(4.4), (Hi) and (Hz), we get y* > Mi. 

The continuous matching conditions at the interface s = Li for the system (4.1) and 
(4.2), are assumed to be 

xi{Li,t) = X2(Ti,t), and = y 2 {Li,t). (4.5) 

The continuous flux matching conditions at the interface s = Li are written as. 


^ dxiiLi,t) dx2{Li,t) 

dl = ds 


(4.6) 
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dyi{Li,t) dy 2 {Li,t) 

ds as 

(4.7) 

The model is studied under two set of boundary conditions i.e. 
In the case of reservoir boundary conditions, we take 

reservoir and no-flux. 

Xi{0,t) = Xl , X2{L2,t) = xl 

(4.8) 

yi(0,t) = 2/i , y2{L2,t) = yl 

(4.9) 

In the case of no-flux boundary conditions , we have 

dxi{0,t) dx 2 {L 2 ,t) 

ds ~ ds 

9yi{0,t) dy 2 {L 2 ,t) 

ds ~ ds 

(4.10) 

(4.11) 

Finally the model is completed by assuming some positive initial distribution, that 

is, for i = 1, 2, 

a:i(s, 0) = Xi{s) > 0, Ti-i <s<Li (4.12) 

^ 1 ( 5 , 0) — ^i(^) ^ 0, Li^l 5 <C Li 

(4.13) 

We first study the existence and stability behavior of the system (4.1) and (4.2) in a 
homogeneous habitat without patchiness, the effect of patchiness will be investigated later. 


4.3 Analysis of the Model in a Homogeneous Habitat 

4.3.1 Model without Diffusion 

In this case Xi = x, yi = y, gi{xi) = g{x), fi{yi) = f(2/), Pi(a;i) = p(x), i = l,2. Thus the 
sj^stem (2.1) and (2.2) reduces to the following form: 

^ = xg{x) - yp{x) 


(4.14) 
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= y[f (y) + rp{^)] (4.15) 

This is a Prey-Predator type model, where the density y of predator species, partially 
depend upon the density x of prey species. In absence of prey, predator grows logistically, 
i.e. predator has an alternative resource. The functions g(a;), f{y) and p(2:) satisfy the 
same type of assumption as mentioned in Hi and H2. And the assumption H3 can be 
written as, 3 or* > 0 and y* > 0 such that. 


a;*g(n:*) - y*p{x*) = 0 

(4.16) 

y*[f(y*) + 7P(3;*)] = 0 

(4.17) 


From equations (4.16) and (4.17) it follows that there are four equilibria, namely (i) 
Eo = [0,0], (ii) Ek = [1^,0], (iii) Em = [0,M] and (iv) E* = [a:*,y*]. The existence of 
E*, the interior equilibrium can be noted from the Fig.4.2. See also [3]. The nature of x* 
and y* with respect to various parameters is shown in the table 4.1, which can be checked 
analytically. 


Table 4.3: Effects of various parameters on x* and y*. 


(i) 

g(0)t 

X* t 


(ii) 

f(0)t 

a:* t 

y* 4- 

(iii) 

Kt 

X* t 

y* t 

(iv) 

Mt 

X* 4- 

y* t 

(v) 

7 t 

X* 4 

y* t 


In particular it is noted from the table that if the carrying capacity of predator species 
(i.e. M) increases then x* decreases and y* increases, (see also Fig.4.2). 

The variational matrix in the general case is given by 
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[M] = 


g{x) + xg'{x) ~yp'{x) 

■yyp'ix) 


-p(x) 

f(y) + yf'(y) + 7P(2:) 


(4.18) 


Now calculating the variational matrix from (4.18) for each equilibria and noting the 
hypotheses Hj — >■ H 3 and the standard stability theory of ordinary differential equations, 
we note the following obvious remarks: 


The equilibrium point Eq is always unstable. The equilibrium point Ek is a stable or 
saddle point according as f(0) + 7 p(iir) is negative or positive. Similarly the Em is a stable 
or saddle point according as g(0) — Afp'(O) is negative or positive. In general, there is no 
obvious remark to be made about the stability of the most interesting non-zero equilibria 
E*, if it exists. 


Therefore, our aim of this section is to obtain local stability and instability as well as 
global stability conditions for E*. 

We now state main results of this section in Theorem 4.3.1, 4.3.2 and 4.3.3 and Lemma 

4.3.1. 


Theorem 4 . 3.1 The equilibrium E* is locally asymptotically stable, if H* < 0, where 


H* 


x*g'(x*) + g(a:*) - y*p'{x*) 


(4.19) 


Proof: Consider a positive definite function V(x,y) with respect to {x*,y*) in the positive 
quadrant, 

V{x,y) = ^ 


7 


.P(^*)> 


{X - xr + piy - y*Y 
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We now compute the derivative of V vdth respect to t along solutions of system (4.14) and 
(4.15), then expand V about {x*,y*), we get, 




(4.20) 


7{x - + (y - y-fe(y-) 


p(x*) 




+ H.O.T. 


= Jix - + {y- y*)^f'iy*) + H.O.T.* 

p[x*} 


(4.21) 


where H* = x*g'(x*) + g(a:*) - y*p'(x*). 


Therefore V(x,y) is negative definite in the neighborhood of E* provided H* < 0. 
Hence the theorem. I 


The next result shows that E* will be unstable if H* is sufficiently positive, in fact 
we can find an explicit lower bound for H* by using Chetayev instability theorem. 


Theorem 4.3.2 If 


then E* is unstable. 

Proof: We now use the instability theorem of Chetayev to prove the theorem. Let 

U {x, y) = \ [(3^ - - Cl (y - y*f] (4.23) 

where Ci > 0. It is noted that U{x,y) > 0 in the interior of a cone C, with principal axis 
parallel to the x-axis and vertex at {x*,y*). As before we now compute U{x,y) along the 
' solution of the system in the neighbourhood of E* to get 


7 p(a:*)p'(a;*) 

H* > J, \ - > 0 
-f'Cr) 


(4.22) 
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The final result of this section, Theorem 4.3.3 gives us the criteria for E* to be glob- 
ally asymptotically stable. First we prove a lemma which establishes a region of attraction 
of the system. 


Lemma 4.3.1 

rant is 


The region of attraction for all solutions initiating in the non-negative quad- 
A = {(x, y):0<x<K,0<y< 


Proof: Prom system (4.14) and (4.15), and hypotheses (Hi)-(H 3 ), we get that x(t) < 
x(t)g(x(t)) < g(0)x(t)[l — x(t)/IC], then lim supt-^oo^(t) < K, for a;(0) < K. Similarly we 
can show that lim supf_^^x{t) K, for a:(0) > K. We now consider y{t) = y{t)[f{y{t)) -\- 
7 p(r(t))]. Then, y(t) < y(t)[f(0)(l - y{t) /M) -hPoo]- Hence 

lim supt_,^y{t) < M[f(0) -f- 7 PDo]/f( 0 ) 

We note that solutions initiating on dA H int 'RA enter into int A. Hence A is an 
attracting set. ■ 


Theorem 4.3.3 Let the condition 

(x — x*)[xg{x) — y*p(a;)] <0, Vi 7 ^ i* 
hold. Then E* is unique and is globally asymptotically stable in A. 


Proof: Let us consider the positive definite function 

r p(0 - P(^*) 


y(s,t) =7 [ 

Jx* 


p(0 


d4+ [y-y* 


y*ln— 

y% 


(4.27) 


(4.28) 
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then V{x,y), along the solutions of the system, is given by 


_ ^p{x) - ^{x*) dx ^ y - y* dy 
V^s,t) - 7 ^ + — ^ 

^ - yp(2^)] + {y- y*)[f(2/) + 7P(a;)] 

= 7— - p|^y ^ [xg{x) - y*p(x)] - - y*)p(2:)] 

+ (y - y*)[f(y) + tpC^:*)] + 7(y - y*)(p(a^) - P(a;*)) 


^P^)_-P(:r-) _ 

x-x* M p(a;) 


. y-y 


Now since 


pW , -p(^ , *) > 0, Vx,^x-andMll£(!Q<o. Vj, / y- 

X — X* y ~ y* 


Therefore, V (x, y) is negative definite if 


{x — x*) — y* <0, Vxy^a;*. 

. PW J 


(4.29) 


The condition (4.27), implies that as long as the portion of the prey curve y = 
xg(x)/p(a;), 0 < :r < X* lies above the line y = y* and that the portion of the prey curve 
y = xg{x) /pix), X* < X < K lies below y = y*, E* is unique [see Fig.4.2]. 
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Now in the following example we show that the existence of E* and the condition 
(4.19) and (4.27) for the local and the global stability are feasible. 

Example: Take g(x) = r{l-x/K), f(y) = s(l -x/M), and p(x) = ax, where 0 < a < 1. 

Then the the equilibrium point E*{x*,y*) is given by 

* sK (r — aM) , rM (s + a'yK) . , ■ 

^ ~ ^ ’ y — ; — provided r > aM. 

rs + a^^KM rs + a^'yKM 

It is also clear from the expression of x* and y* that the results shown in the previous 
table is valid, where g(0) = r, f(0) = s and rest are same. Further if a; increases then x* 
decreases and y* increases, provided x* > Kjl. 

Now we note that the conditions (4.19) and (4.27) are automatically satisfied, as 
E* = -Tx*jK < 0 and 


(x - x*)[xg(x) - y*p(x)] == 


T 


x(x — X) < 0, Vx ^ X*. 


Hence E* is locally as well as globally asymptotically stable. 


4.3.2 Model with Diffusion 

Now, we wish to consider the model (4.14) and (4.15) with diffusion and analyze the uniform 
equilibrium E* under both reservoir and no-flux boundary conditions. The mathematical 
model in this case can be written as, 

^ = xg(x) - yp(x) + (4.30) 

^ = y[f(y)+7P(:r)]-Hi^20 (4-31) 

•' 0 <s<L 

where A > 0 are diffusion coefficient and 7 > 0 is the conversion of biomass constant. The 
reservoir and no-flux boundary conditions are respectively. 


x(0,t) = X* = x{L,t) and y(0,t) =y* = y{L,t), 


(4.32) 
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^ Q ^ t) dy{0, t) ^ Q ^ dy{L, t) 
ds ds ds ds 


(4.33) 


Now if we take the following positive definite function 

W{x,y)= [ V{x,y)ds, 

Jo 

around the equilibrium point then using similar arguments as in previous section 1. , 
we get the following theorems for local and global stability of the system. 


Theorem 4.3,4 The equilibrium E* is locally asymptotically stable, if H* < 0, where 


H* = x*g\x*) + g{x*) - y V(a:*)- 


I 


It is pointed out that even if iJ* > 0 this uniform equilibrium may become stable 
with diffusion under a condition as given in the following theorem. 

Theorem 4.3.5 Let H* > 0. Then the equilibrium E* is locally asymptotically stable, if 
H* < 


Proof of Theorem 4.3.4 and Theorem 4.3.5 : Consider a positive definite function 
W(x,y) with respect to {x*,y*) in the positive quadrant. 


W{x,y)= [ 
Jo 


i 1 
2 


7 




{x-xr + ^iy-y*? 


ds 


(4.34) 


We now compute the derivative of W with respect to t along solutions of system (4.14) and 
(4.15), then expand V about (x*,y*), we get, 


Wix,y) = [ 

Jo 


'Y{x — X*) 


p'(a:*)^ dx {y-y*)dy 


TTT + 


p(:z ;*) ) dt y* dt 

rL 


ds 


= f C{x- x*fH*ds + f {y- y*fr{y*)ds 
Jo Jo 

+ / CDi{x-x*)^ds+ D,-^^ds + H.O.T. 

Jo os^ Jo y- y* os^ 


(4.35) 
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where C = 7 ^ and H* = x*^{x*) + g{x*) - y*p'ix*). Again 


By using Poincare’s Inequality, we get 
-CDi 


t 2 




ds (Using Integration by Parts) 

(4.36) 


1" 

d{x — X*)' 

lo 

ds 


ds>—CDi^f (x — x*)'^ds 
Jo 


(4.37) 


Therefore W{x,y) is negative definite in the neighbourhood of E* provided H* — 
< 0. Hence the theorem 4.3.5. I 

Theorem 4.3.6 Let the condition 

{x - a:*)[a;g(x) - y*p(x)] <0, Vx x* 
holds. Then E* is unique and is globally asymptotically stable in A. 

Proof: Here we consider the following positive definite function 

W{x,y)=[ V{x,y)ds, (4.38) 

Jo 

around the solutions space of the above system, where V{x,y) is same as in section 1.4.1. 
Differentiating (4.38) with respect to t, we get 


W{s,t) = 7 [ 

Jo 


i p(x) - p(x*) 5x , f^y-y*dy 


p(x) dt 


ds + 


i: 


y dt 


ds 


/•i p(x) - p(x*) _.Jxg(x) 

= Vo X-X- ■ * + ^ 


f{y)-{{y*) 


y-y* 


ds 


+ 


Jo p(x) OS Jo y os 


p(x) os~ Jo y 

How integrating by parts and using (4.32) or (4.33), we have 


(4.39) 


h 


= 7D1 f 
Jo 


^ P (^) ~ P(^*) 

p(x) ds^ 


di 


= -7A / 
Jo 


^ p(x*) . ( dx\ 

2? iP ( Fi ) ^ 

p2(x) \ds J 
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and 


Hence the result. 




X 


^y -y* d^y 
y ds‘ 


ds = 


-D2 [ 
Jo 


i as , 


ds < 0. 


Thus the uniform equilibrium E* of the system with diffusion is locally as well as 
globally asymptotically stable under the same set of conditions, as in the previous case. 
Further by Theorem 4.3.5, an equilibrium which is unstable without diffusion can become 
stable with diffusion. 


4.4 Analysis of the Model with Diffusion in a Two- 
Patch Habitat 

In the first subsection we study the model (4.1)-(4.13), in the case of uniform steady state. 
Then in next subsections, we will study the non-uniform state state case. 

4.4.1 The Uniform Steady State under Both Sets of Boundary 
Conditions 

The main purpose of this section to show that uniform steady-state in the two patch habitat 
is globally asymptotically stable. In this case, it is clear that, under both sets of boundary 
conditions, there is a uniform steady-state, = K*, 0 < s < L 21 1 > 0, yi{s,t) = M*, 

0 < s < L 2 , t >0, (i=l,2) where K* and M* are the common uniform equilibrium of first 
and second species respectively in the two patches. Here = K* and 2/1 = 2/2 ~ 

By using similar arguments and proofs as in section 1.4, the following theorems for stability 
can be proved for the system (4.1) -4 (4.7) with (4.8) and (4.9) or (4.10) and (4.11). 
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Theorem 4.4,1 The equilibrium (K*,M*) is locally asymptotically stable, if 

for i = 1,2, (4.40) 

and the following conditions are satisfied, for i — 1,2 

YriM*v'i{K*) - PiiK*)f < (4.41) 

where 

H* = + gi[K*) - M*p\{K*) (4.42) 


Proof: Lineariazing the system (4.1)-(4.2), by asing 


we get, 


Xi{s,t) = K* + ni{s,t) 
yi{s,t) = M* + mi(s,t) 


71 ’ 

= ni[gi{K*) + K*g!,{K*) - M%{K*)] - m.Pi(A'*) -h Dir ‘ 






dt 

We considering the following positive definite function, 

1 ^ r^' 


ds^ 


V(i) = ^Er [ixi-K*f + iyi-Mr] 

Z ^ J Li^i 


where d,, i = 1, 2 are positive constants. 
Differentiating (4.45) and using (4.43)- (4.44), we get 

^ r^i „ ~ /■t'i 


y = 


E f ' nfi?;ds-f-E f ‘ m^M^fl{M*)ds 

^ J Li^i I J Li^\ 

+ hr riimi[-pi{IC) + 'riM*p',iK*)]ds 

, d^Uij , ^ 


(4.43) 

(4.44) 

(4.45) 


(4.46) 
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Again using integration by parts, for both 

4^ f^‘ ^ 


ds 


ds^ 

t r 


2 rLi 


ds^ 


type of boundary conditions, we get. 



Therefore, from (4.46), V is negative definite, if the conditions (4.40) and (4.41) are true, 
for 2 = 1,2. I 


Moreover, 


Theorem 4.4.2 Let H* > 0. Then the equilibrium (K*,M*) is locally asymptotically 
stable, if the conditions ( 4 - 41 ) along with 

fori = l,2 

[Li — Li-iy 

hold. 


Proof: By using Poincare’s Inequality, we get 


f^i 

Du / 

JLi-i 


( dui ' 
ds 


ds < Du 


TT 


{Li - L 


— r 

i_l)2 


n,-ds 


Therefore from (4.46), we get, 
i, rLi 


^ rX/i ^ rLi ( Qrrr * \ ^ 

+ 12 nimi[-pi(ir*) -h7iM*p-(iir*)]ds-^r»2i^ V~^/ 


V < 


Hence the theorem. 


I 
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We note here that the stability of uniform equilibrium depends on two respective 
conditions on and H^- Further, the behavior of HI, i = 1,2, with respect to K* is 
similar to H* as shown in Fig.4.1. Hence it may be concluded here the patchiness has a 
destabilising effect if the change in any parameter of the system caused by patchiness due 
to industrialization, decreases K*. 

We now state the global stability of the uniform steady state. 

Theorem 4.4.3 The uniform steady-state (K* , M*) is globally asymptotically stable if 

{xi - K*)[xigi{xi) - M*pi(xi)] < 0, Vxi ^ K*. (4.47) 


Proof: Let V(x,y) be the positive definite function about x = K*, y = M*, given by 

vi.,y) = /; (4-48) 


Differentiating (4.48) with respect to t, and using (4.1) and (4.2) we get, 

dx, 
dt 




^ fi-i ( y i-M* \ 


r - K') 


a;igi(xi) 


L Pi(^i) 


Vi 


M* 


dt 


ds 


tl-Du jf 


+ 


+ 


ds 


yi - M* 

Li Pi{Xi)--pi{K*)d^X^^^^ 2 


Pi(a:i) ds 

Now consider the integral, 

2 rLi pi{xi) - Pi{K*) cf^Xj 


? a? ^ ^ ^ 


h 




Pi (a:,) 


ds^ 


ds 


= E^iAi 

1 


Pi(a:i) - Pi{K*) dxi 


Pt(2;i) 


ds 


• -f 

JLi-i 


Li Pi(ii:*)p;(a;i) f dxf 


Pii^i) 


ds 


ds 
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Y pf{Xi) \ds J 

4- n^^fr -f) L ~ P i(-^*) p2(a;i(I'i,i)) — P2(-K’*) 

= -El, A. f' 

Y JLi-i Pti^i) \ds J 


because under both set of boundary conditions, 

dx-i 
ds 


[x,(o,t) - ir*l^(o,t) = = 0 


(4.50) 


(4.51) 


Similarly, 




Hence V{x,y) < 0, and V{K*,M*) = 0. Therefore V{x,y) is negative definite over 
x> 0,y> 0 with respect to a:* = K*, y* — M*, proving the theorem. I 


This result generalizes the result of [9, 11] for a patchy habitat. 


4.4.2 The Non-Uniform Equilibrium State 


Now we consider the model (4.1)-(4.13) and assume the existence of the unsteady state 
solution of the system [see [22](Theorem 1, page 111 and Theorem 3, page 123)] and also 
[33]. Our aim here is to show that there exists a positive, monotonic, continuous steady- 
state solution for the each species separately, with continuous flux at the interface under 
both reservoir and no-flux boundary conditions and derive the conditions for asymptotic 
stability of non-uniform equilibrium state. 

There are fom possible cases: 

(1) xl > xl and yl>yl 

(2) X 2 > xl and y 2 < Vi 

(3) rc 2 < xl and yl < yl 

(4) xl < rcl and yl > yl 
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Without loss of generality, (3) can be reduced to (1) and (4) can be reduced to (2). 
Therefore we study the cases (1) and (2) separately under two set of boundary conditions. 

We assume the following conditions on Xi and yj for i = 1, 2. 

- 2 :*)[ 2 ;igi(a;i) - yiPi(a:i)] < 0, Vrcj x* and < Vi < (4.53) 

iVi - yi)[yi^i{yi) +7t1/iPi(^i)] < o, Vyi ^ y* and min.{xl, x^} < x,- < max.{x\,X 2 } (4.54) 
The analysis in the remaining part of this chapter is valid under these two conditions. 

4.4.3 The Model under Reservoir Boundary Conditions: When 

x *2 > x\ and > yl- 

We first consider the steady-state problem and show that there exists a non-uniform positive 
mono tonic solutions Ui{s), u,(s) under reservoir boundary conditions and with continuous 
flux at the interface. 

The steady-state system takes the form 


Du ^ -l-Uigi(ui) UiPi(wi) = 0 

(4.55) 

D2i ^ + VifiiVi) -b 'nViPiiUi) = 0 

(4.56) 

The reservoir boundary conditions are 


Ui{0)=x{,U2{L2)-X2 

(4.57) 

vi{0) = yl,v2{L2) = y2 

(4.58) 

The continuous solutions and flux matching condition at the interface 


Ui{Li) = U2{Li);vi(Li) = V2{Li) 

^ dui{Li) „ du 2 {Li) 
ds ds 

(4.59) 

(4.60) 
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(4.61) 


n dvi{Li) dv2{Li) 

^21 -j- — — JJ-n — j 

as as 

Let pi(s, Q!i) and gi(s, /?i),0 < s < Li be the unique solutions of (4.55) and (4.56) respec- 
tively, for i=l, such that 

5pi(0,a:i) 


99i(0,)di) 


ds 


= /3i, 9i(0,/?i) = yi 


(4.62) 


(4.63) 


Let P 2 (SjQ! 2 ) and q 2 {s, P 2 ),Li < s < L 2 he the unique solutions of (4.55) and (4.56) 
respectively, for i=2, such that 

dp2iL2,a2) 


ds 

dq2{L2,^2) 

ds 


= 02, P2{L2,a2) = X*2 


= 1^2, 92 (^ 2 , /?2) = 2/2 


(4.64) 

(4.65) 


We will have shown the existence of the monotonic solutions if we can show that there 
exists Oii,a 2 ,Pi and P 2 such that 


Pi(-Li,Q:i) =P2(-Li,q:2) , qi{Li,/ 3 i) = 92(^1,,^) 

n dpi{Li,ai) _ ^ dp2{Li,a2) ^ dqi{Li,^i) dq2{Li,02) 

;;; = ^12 , 4^21 ;;; = ^ 12 - 


ds 


ds 


ds 


ds 


(4.66) 

(4.67) 


From (4.55) and (4.56), multiplying both side by 2dui/ds and 2dvi/ds respectively 
and integrating with respect to s from 0 if i=l and from L 2 if i=2 , we get. 


2 

a? - — / [nigi(ui) - ViPi{ui)]dui{s) 
JJli Jx\ 

■ ^ 2 /■"<(«) 

z=0f-— [Uifi(Ui) -H 7i^^iPi(44i)]«i^'i(s) 
442* 4yr 


(4.68) 

(4.69) 


Keeping in view of the conditions (4.53) and (4.54), we note here that pi{s,ai) and 
9i(s, A) also satisfy the same conditions. 


In order to construct our required results we need some preliminary lemmas. 
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Lemma 4.4.1 Ifa^ Pi > 0, then 
dpi{s, ai) 


ds 


> q;i, and > A, on 0 < s < Li- 

cfs 


Proof: Prom (4.68) , we get, 


api(s, O i) 

ds 


l2 


2 rPi{s,ai) 




11 


rpi{s,ai) 

/ bigibi) ~ 9iPi(Pi)]^^Pi (4.70) 

Jxl 


Now using (4.53), it follows from (4.70), 


dpi(s, ai) 


ds 


> a\. 


Since 


Hence 


dpi 


Pi(0, Oi) = Xj, -^(0, Oi) = Oi and Pi(s,q:i) > xL for 0 < s < Li. 

OS 


dpi{s,ai) 

ds 


> Oi, VO < s < Li. 


Similarly from (4.69) , we get. 


dqi{s,Pi) 


ds 


and by using (4.53), we get 


^ 0 

= Pi - / [gifi(9i) + 7i9iPi(Pi)]rf9i (4.71) 

J-J 21 Jyl 


dqi{s,Pi) 

ds 


> pi, VO < 5 ^ Li- 


Hence the lemma. 


Lemma 4.4.2 If OC 2 , P 2 > ^) ^ < P 2 < ^2 0 < 92 < ^ 2 ? then 


3P2{s.o:2) ^ i,<s<L,. 


ds 


ds 
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Proof: Prom (4.68), we have, 


dp2{s, OL 2 ) 


t 2 


ds 


(Xn 


rP2 (5,0:2) 

/ \P2S2(j>2) - q2P2{P2)]dp2 (4.72) 


2 fP2{s,a2) 

Di2 Jx 


Hence by using (4.53) we have, dp 2 {s-,Ci 2 )/ds > Oi 2 ^ Li < s < L 2 for 0 < ^>2 < ^ 2 - Similarly 
from (4.69) and (4.54), implies dq 2 {s, h)/ds > P 2 , Li < s < L 2 for 0 < q 2 < y^- ■ 


Now we define four continuous maps by the following lemmas, proofs are analogous to 
lemma 2.3.3 of chapter 2. 

Lemma 4.4.3 Define Fii(ai) by Pii(a;j) = pfiLi^ai). Then there exists di > 0 such that 

Pii : [0,0:1] — f [xijX^j 

Fi2 : [0,0:2] 


Lemma 4.4.4 Define F 2 i{fii) by P 2 i(A) = qi{Li,Pi)- Then there exists fii > 0 such that 

F21 ■ [0, a ] — >■ [2/1,2/21 
F22 ■ [0, ^ [2/21^1] 


In the following theorems we have shown the existence of non-uniform steady state 
solutions. 

Theorem 4.4.4 There exists a continuous, monotonic solution u, of system (4-55) with 
continuous flux at L\ (i.e. under (4-d6) and (4-d7) ). 
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Proof: Prom Lemmas 4.4.1 and 4.4.2, it follows that any solution we construct must be 
monotonically increasing. Again by Lemma 4.4.3, for each 0 < Q 2 ^ d; 2 , we can find an 
Oil such that 0 < oii < di for which = P 2 {I'i,ci 2 )- Hence Oi can be solved as a 

function of 0 : 2 , ai = h{a 2 ), to give a continuous solution of (4.55) with (4.57), (4.59) and 
(4.60). 

Let 

= (4.73) 

Clearly 0 ( 02 ) is continuous on 0 < 0:2 < 0 : 2 . Then we have 

G(0) = Ai— ^ >0, 

and 

G(aa) = < 0. 

Therefore 3 0 : 2 , 0 < 0:2 < d; 2 , such that 0 ( 0 : 2 ) = 0. Hence the theorem. I 

Exactly similar result is also holds true for the second species. 

Theorem 4.4.5 There exists a continuous, monotonic solution Vi of system (4-56) with 
continuous flux at Li (i.e. under (4--66) and (4-67) ). I 


We now study the linear and non-linear asymptotic stability of the steady state (4.55) 
to (4.61) with respect to the system (4.1) — >• (4.13). 

Theorem 4.4.6 The steady-state, continuous, monotonic solutions of the linearised sys- 
tem ( 4 . 55 ) and (4-56) with continuous flux at the interface s = L\ is asymptotically stable 
provided the following conditions are satisfied: 

(z) ^li(tti, Vi) = Uig'iiUi) -b gi(Ui) - ViPi(Ui) < 0 , 

(ii) Q2i{Vii Vi) = Vifl(Vi) -b ft(z^t) + TiPiC^i) — 
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(«*) [Pz(^ii) - < 4guQ2i 

for x*^<Ui< xl, yl<Vi< yl and i = l, 2 . 


Proof: Let the steady-state solution of system (4.55) be 

^ / ^i(s)> 0 < s < Li 
^ ^ ( U 2 (s), Li < s < L 2 

And let the steady-state solution of system (4.56) be 

_ / ^ 1 ( 5 ), 0 < s < Li- 
^ ^ 1 ^^ 2 ( 5 ), Li<S<L2 


Lineariazing (4.1) and (4.2) by using , 


Xi(s,t) = Ui(s) + ni(s,t) 

(4.74) 

yi(s,t) =Vi{s) + mi{s,t) 

(4.76) 


We have, 


dUlihA = Ui {gi{ui) + Uigi{ui)) - niVip'i{ui) - miPi{ui) -f 


dt 


ds^ 


= rni (fi(ui) + Viflivi)) + maiPiiui) + niUi7iP-(ui) + 


(4.76) 

(4.77) 


Using (4.74) and (4.75) the corresponding initial, boundary and matching conditions can 
be obtained as follows 


ni{s, 0) = Xi(5) - Ui{s), mi{s, 0) = 5i{s) - Vi{s) 


ni{ 0 ,t) = 0 = n2iL2,t) , Tni(0,t) = 0 = m2{L2,t) 
= n2{Li,t) , Tni{Li,t) = m2{Li,t) 
Dn^{Li,i) = D,,^{Lut) 


ds 
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Now we consider the following positive definite function, 

From which we get, 

By using (4.76) and (4.77) , we get. 


ds 


(4.78) 


j [gt(wi) + Mig((ui)]ds-^ / n\vip[{ui)d. 

1 I JLi^i 


2 rLi 


^ fU 2 i. 

^ niTTiiPi (Ui)ds + Y^ m2 [fi (Ui) + Vifl (u,)] ds 

/ ’7r,miUip((uf)ds-l- ^7i f m?p,(uj)ds 

1 I •'Li-1 

n ^ T-^ 7^* 


+ 


5mi 


(4.79) 


Therefore V'(t) is negative definite, if conditions (i), (ii) and (iii) are satisfied. Hence the 
theorem is proved. | 


Theorem 4.4,7 The steady-state, continuous, monotonic solutions of non-linear system 
(4-55) to (4-61) with continuous flux at the interface s = Li is asymptotically stable in 
the subregion o/R: {rcj < Xi, Ui < xl, yl < yi, Vi < y|, for i = 1, 2} provided the following 
conditions are satisfied: 


(1) ~ - t^i(p»(a^t) - Pt(Mi)) ^ ^ 


Xi — Ui 




Vi - Vi 

UiPiiXi) Pi(^i) “ Pi(“i) ' 


^^2 


Xi 


Vi 


Xi — Ui 


UiVi 
XiPi 


( 3 ) 
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Proof: From (4.1), (4.2), (4.74) and (4.75), we have 

-^z=(ui + ni)gi{ui 4- Tii) - Uigi{ui) - {vi ■+- mi)pi(ui + rij) + ViPi(ui) + (4-80) 

dlTl' 3"^7Tt' 

= (Ui -f mi)ii{vi -f rrii) - Vifi{vi) +7i(t)i + Tni)pi{ui + Ui) - 7 iniPi(ui) -f (4-81) 

We now consider the same positive definite function (4.78) as in Theorem 4.4.6. Differen- 
tiating (4.78) with respect to t and after using equations (4.80) and (4.81), we get 

— - r {uj 4- ni)gi{ui + Tij) - Uigiiuj)' 


m = E r — 

^ JLi^i [_ 


rii 


ds 


- ?/:. 


n: 


Vi 


'Pi{ui + ni) -piiui)' 


rii 


ds 


mf 


(t>i 4- rni)fi{vi 4- mj) - Vifijvi) 


^ ^ I'Ll ^ ^ rLi 3^Tli 

- E / nimiPi(ui 4- ni)ds 4- E-^ii J 

- r 

j hi^i 

+ Ylli ! 

J Li^i 


rrii 

^Pijui + nj) -piiui)' 


4- C£7iPi(ui 4- Ui) 


ds 


4- 


d'^rrii 


E -^21 / m 


ds 


Now, 


ds^ 

Pi{ui -4- rii) - Pi(^i) 


> 0, for all values of rii. 

Ui 

Also by using boundary, matching conditions and dui/ ds > 0, we have 


(4.82) 


(4.83) 


(4.84) 


Similarly 


Therefore V{t) < 0 if the conditions (1),(2),(3) holds true. Hence the theorem is proved. 


I 
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It is noted that if we linearize the conditions of theorem 4.4.7, we get the conditions 
of theorem 4.4.6. 


Remark: When Du —¥ 0, for i — 1,2. In such a case, from (4.55), we have, u,- = 
Wtgt(^it)/Pi(wi). Then the stability condition (i) of theorem 4.4.6 modified as follows, 

Si{u<) = uig-fe) + 

We note that for fixed s, the behavior of the function Qi{ui) with respect to u,- is same as 
that of H*{x*) relative to x*. Hence Qi{ui) also becomes less negative if Ui decreases [see 
Fig.4.1]. In particular if Ki decreases then x* decreases, therefore for a fixed s, ui decrease 
and Qi{ui) becomes less negative. Thus, the effects of those parameters such as Ki, Mi, 
which decrease Ui, are destabilizing. This implies that if carrying capacity of the prey re- 
source in a particular patch decreases due to depletion of that patch then it is destabilizing. 


4.4.4 The Model under Reservoir Boundciry Conditions: When 
> xl and yf > 

It can be checked that for the monotonicity of steady state solution is similar as 

the above case and for yl > y 2 , the monotonicity is just the reverse of the above case, [see 
also Fig.4.5]. 

In this case also the linear and non-linear asymptotic stability conditions for the 
steady state system under reservoir boundary conditions are similar as in Theorem 4.4.6 
and Theorem 4.4.7. 

The above theorems imply that the system will settle down to the steady state dis- 
tribution in the two patches under certain conditions, the magnitude of the steady state 
distribution of the prey density being lower than it’s original value but the density distribu- 
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tion of the predator species being correspondingly higher. The above analysis also suggest 
that patchiness causes destabilization of the system. 

4.4.5 The Model under No-flux Boundciry Conditions: When 

> xl and > yl 


The steady state system becomes, 


r. (Puiis) 


+ Uigi{ui) - ■yiPt(ui) = 0 


D 


CpVi(s) 
ds^ 

The no-flux boundary conditions are, 


+ Vifi(vi) + 'riViPi{ui) = 0 


dui{0) du2{L2) 

ds ds 

dvijO) _ Q _ dv 2 {L 2 ) 
ds ds 

The continuous and the matching conditions at the interface are. 


Ul{Ll) — U2 {Li), Vi{Li) — V2 {Li) 


Diiu[{Li) = Di2U2{Li), D2iv[{Li) — D22V2{Li)- 


(4.86) 

(4.87) 

(4.88) 

(4.89) 

(4.90) 

(4.91) 


Let pi{s,ai) and q\{s,0i), 0 < s < Xi are the unique solutions of (4.86) and (4.87) 
respectively, for i = 1, such that 

, Pi(0,ai)=ai (4.92) 

as 

= o ■ ft(0,A)=A (4.93) 

Let p 2 {s, 02 ) and q2is,0), Li < s < L 2 axe the unique solutions of (4.86) and (4.87) 
respectively, for i = 2, such that 

= 0 , p2(L2,a2) = a2 (4.94) 
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’ 9'(i2-A)=A («5) 

We will have shown the existence of the monotonic solutions if we can show that there 
exists tti, a 2 , Pi and P 2 such that 


P\{Li,oli) — P2{Li,(X2) , (ii{Li,Pi) — q2{Li, P2) 


dpi{Li,ai) dp2iLi,a2) 

— Fs — = — Fs — ’ 


^ dqi{Li,Pi) dq2{Li,P2) 

Dn g; = -Du— gj 


(4.96) 

(4.97) 


From (4.86), multiplying both side by 2dui/ds and integrating w.r.to s from 0 if i=l and 
from L 2 if i=2 , we get 


duj 

ds 


2 



njPi(ui)]d;ui(s). 


(4.98) 


Similarly from (4.87), we have 

= “TP / [^;ifi(ui) -|-7iUiPi(ui)]diii(s). (4.99) 

lJ2i JPi 

Here also we need similar type of lemmas as in the case of reservoir boundary conditions. 


dvj 

ds 


Lemma 4.4.5 If ai > x^, Pi > Pi, pi > and qi > yl, then 




ds 


ds 


Proof: From (4.98), we get, 


dpi{s, ai) 


ds 


^ 2 rp^ 

Dll •'ai 


7>i(«,ai) 


\pigi(pi) ~ qiPi{pi)]dpi 


(4.100) 


Since dpi{0,ai)/ds = 0, Pi(0,ai) = ai there exists si > 0 such that Pi(s,Q;i) > ai on 
0 < s < Si. If not, let so, 0 < sq < Li, be the first positive value , if it exists, such that 
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Pi(so, Q:i) — Oi\. Then by the mean value theorem there exists s such that 0 < s < so 
dp\{s, ai)fds = 0; that is, 

0 = / [pigiCpi) - giPi(pi)] dpi (4.101) 

Ull Joti 

But the right-hand side of (4.101) is nonzero, since pi > cti on 0 < s < s. Also cti > a:*, 
therefore, pi > xl and since qi > y*. This imphes that, Pigiipi) — giPi(pi) < 0, giving a 
contradiction. Hence pi(s,Q;i) > cti, for all value 0 < s < Li. This implies, p{s, ai) > 
on 0 < s < Ti. Hence by (4.47) and (4.100), the first part of the lemma. Similarly, other 
part of the lemma follows. I 

Lemma 4.4.6 If 0 < ^ 2 ( 5 , eta) < X 2 and 0 < ^ 2 ( 3 , /If) < then aa < implies 
P 2 {s, af) < OKa and ^ implies g 2 {s, ^f) < for Lx<s< I-a- 

Proof: Since 0 < p 2 {s,<X 2 ) < xf then from (4.53), we get, P 2 g 2 (P 2 ) — y 2 P 2 (P 2 ) > 0. Now 
Since 0 < 52 ( 5 , A 2 ) < and 02 < xf then P 2 g 2 (P 2 ) — g 2 P 2 (P 2 ) > 0, for all 02 , P 2 - Hence, 
from (4.98), P 2 < oia- Again by using (4.99) the follows other part of the lemma. I 

As in the case of reservoir boundary conditions, 3 four continuous mapping as follows: 
Lemma 4.4.7 Define Guioi) by Gii{ai) =pi{Li,ai). Then 3 a,- such that 

Gii : [a;i,Q:i] -4 [x\.,x!^, 

Gi 2 : [0:2, 0:2] 

I 


Lemma 4.4.8 Define G 2 i{Pi) by (? 2 t(A) — gi(Li, ^i). Then 3 0i such that 

Gai : [2/1 j a] [yiiV^]’ 

G22 ■■ [h,y^ [yl.yl]- 
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Proof: Same as lemma 2.3.6, (chapter 2). I 

The existence theorems are given as follows: 

Theorem 4.4.8 (i) There exists a continuous, monotonic solution of system 
continuous flux at Li. 

(ii) There exists a continuous, monotonic solution of system 
continuous flux at Li. 

Proof: Analogous to theorem 1.5.8 and theorem 1.5.9. I 

Now we state the stability conditions of the system with no-flux boundary condition. 
Proves are similar to the proofs of theorem 1.5.10 and theorem 1.5.11 . These conditions 
are the same as in the case of reservoir boundary conditions. 

Theorem 4.4.9 The steady-state, continuous, monotonic solutions of system (4-86) and 
(4-87) with continuous flux at the interface s = Li is asymptotically stable provided the 
following conditions are satisfied: 

(i) Vi) = Uigi{ui) -f gi(ui) - ‘Uip'(ui) < 0 , 

(w) g 2 i {ui, Vi) = vfli (vi) -f fj (ui) -f 7 , Pi (ui) < 0, and 

(Hi) IPi(ui) - 7iVtP'i(^)f < ^Q\iQ2i 

for xl<Ui< X 2 , vl < Vi < yl and i = l,2. 

Theorem 4.4.10 The steady-state, continuous, monotonic solutions of non-linear system 
(4-86) to (4-91) with continuous flux at the interface s = Li is asymptotically stable in the 
subregion o/R: 


(4-86) with 
(4-87) with 


{xl < Xi, Ui < X 2 , yl < yi, Vi < y^, for i = 1, 2} 
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provided the following conditions are satisfied: 

(1) ~ ~ Vi{Pi{Xi) — Pi(ni)) 

Xi — Ui 

(2) J^2i = + 7iPi(a:i)] - + 7iPi(xi)] 

Vi -Vi 

(3) [ifiPiM _ - Pi(^») 

Ui Y Xi Iti 


In this also the similar results for the steady state distributions are valid as in the 
case of reservoir boundary conditions. 

4.4.6 Both the Species have Uniform Equilibrium State in the 
Second Patch 

In this case, the steady-state solutions of both the species are variable in the first patch 
and constant in second patch, i.e. U 2 = X 2 and V 2 = 1 / 2 , Li < s < L 2 , t > 0. As shown 
in the general case we note that here also the steady state solution is positive, continuous 
and monotonic in the first patch [ see Fig.4.8 ]. For stability analysis we use the positive 
definite function, 

y(t) = M [^\xi - uifds + [ \x 2 - X 2 fds+ [ \yi- vifds + f (y 2 - y^fds 
2 I Jo Jo JLi 

and analyze in the similar manner as in Theorem 4.4.6, and can prove the following theorem. 


l2 


^ . V-iVi ^ ^ 


Theorem 4.4.11 Let U 2 = x^ and V 2 = yl- Then the steady-state, continuous, mono- 
tonic solutions of non-linear system (2.1) to (2.4) and with either (2.5) or (2.6) and with 
continuous flux at the interface s = Li is asymptotically stable in the subregion of R: 


{xl < xi, ui < X 2 , y{ < yi, vi < ya} 
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provided the following conditions are satisfied: 

(*-^) •wi) = ■uigi(ui) + gi(ui) - 'yip'i(ui) < 0 , 

(^•^) ^2 = ^*2S2i^l) + g2{xl) - y;p'2i^*2) < 0, 

(^0 = ■WifiK) +fi(vi) + 7iPi(iti) < 0, 

(***-a) [Pi(ui) - 7i^iPi(«i)]^ < ^Qi^i 

[P2(^2) - 72y|P2(a:;)]^ < 


It is noted here that these condition are less in number then the conditions of the 
theorem 4.4.6 and more than theorem 4.4.3. Thus the distributions in this case are more 
stable than the general case but less stable then the uniform distribution in the two patches. 

4.5 Numerical Examples 


In the following the steady state solutions of couple-reaction difiusion system in a two-patch 
habitat with reservoir boundary condition and fl^lx matching condition at the interface are 
obtained numerically to show there existence. For this we consider the following particular 
form of functions: 

gi{ui) = ru{l-Ui/Ki), fi{vi)=T2iil -Vi/Mi), Pi(ui) = ; , i = l,2. 

Clt I OiUi 

Then the steady-state system (4.55) and (4.56), i=l,2, can be modified as follows : 

Du^ = - ^/Ki) - (4.102) 

ds^ Oi + biUi 
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= r2iVi{l - Vi/Mi) + jiVi- 


The reservoir boundary conditions are: 


Uj I 


(4.103) 


ui(0)— aji, U2{L2)=x*2, ui(0) = y*, V2{L2)=yl (4.104) 


The continuous flux matching conditions at the interface are 

dui{Li) du2{Li) ^ dvi{Li) ^ dv2iLi) 

= U2(Li), vi(Li) = n2(Li) 


(4.105) 

(4.106) 


We solve system of equations (4.102) to (4.106) numerically by using finite difference 
method in various cases using different set of values for parameters satisfing the stabil- 
ity criteria. The corresponding graphs are shown in Figures 4.3, 4.4, 4.5, 4.6 and 4.7. 


Prom the graphs it can be noted that in all cases the steady state solution is monotonic. 
In particular from Fig.4.3-4.6, we observe that for all possible combination of x*’s and 
y*'s, the steady state solutions are always monotonic. In Fig.4.7, the effect of diffusion 
coefficients on the non-uniform steady state distributions is shown. It is noted from this 
figure that as respective diffusion coefficients in the two patches increase the steady state 
solutions tends towards linear distributions. In Fig.4.8, it is shown that if the stecidy state 
solutions of both the species have uniform equilibrium state in the second patch, then the 
steady state distribution is monotonic in the first patch for both the species. 
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Figure 4.7: Effects of diffusion coefficients on the stead 
v'alues of D^-’s. 

3ofid curves for Dij's = 0.04, dotted curves for Aj’s = C 
= 2.00, for bi = l.OOOE - 6, 7i = 1.0, z = 1, 2 and 
ri = 0.01, ii:i = 50, ai = 2.000£^ - 05, ai = 0.5 

r2 = 0.02, ii:2 = 100.5, aa = 2.000E' - 05, 02 = 0.5 

51 = 0.02, Ml = 100, ai = 2.000E’ - 05, oi = 0.5 

52 = 0.03, M 2 = 150, a2 = 2.000E - 05, 02 = 0.5 
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Interface ■ 



Igure 4.8: The Steady State distributions when the both the species have uniform equi 
ibrium state {K 2 = 2^2 = 120, M 2 = yl = HO) in the second patch. 
ri = 0.174, ii:i = 58, ai = 2.0000^^ - 04, oi = 0.2 
Si = 0.316, i\\ = 79, as = 2.0000E’ - 05, oi = 0.2 
1: = 07 Tt L-i 1 1 0.6, hi = l.OOOO^; - 05, 7i = 1.0, Vi, j = 1, 2. 
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4.6 Summary 


In this chapter, a dynamical model for two logistically growing species with prey-predator 
tJTP® interaction with diffusion in a two-patch environment is studied by assuming the rate 
of change of density of the first species decreases as the second species increases and the 
growth rate density of the later increases with the increase of density of the former. This 
model is proposed by keeping in view the situation in Doon Valley (Uttar Pradesh, India) 
where the depletion of forest biomass has increased patchiness in this valley by grc wth of 
biomass dependent industrialization (population) [see Shukla et al. [ 25 ], Munn et al. [ 17 ] 
for detail information]. 

To understand the role of diffusion and patchiness the proposed model has been 
analyzed for both homogeneous and two patch habitats. In absence of patchiness in the 
aabitat an equilibrium which is stable locally as well as globally without diffusion under 
certain conditions remains stable under the same set of conditions with diffusion. However if 
m equilibrium is unstable it may become stable with diffusion in the homogeneous habitat 
inder a certain condition. In the case with diffusion in two patch habitat the uniform 
iteady state remains stable under two similar conditions as in the case of homogeneous 
labitat with diffusion. 

In the case of non-uniform steady state, it is shown that there exists a positive, mono- 
onic, continuous steady state solution with continuous flux at the interface of two adjoining 
latches, under both the reservoir and no-flux boundary conditions, for each species, that is 
symptotically stable in the linear and non-linear cases under some conditions. Also when 
he steady state solutions of both the species are non-uniform in the first patch and uniform 
1 the second patch, then the stability of the system increases (decreases) as compared to 
le case of non-uniform (uniform) steady state in both the patches respectively. Thus it is 
Dncluded that patchiness, destabilizes the system. 
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Chapter 5 


A Prey-Predator Type Model with 
Diffusion and a Supplementary 
Resource for the Prey in a 
Two-Patch Habitat 

5.1 Introduction 

As mentioned in chapter 4, an interesting problem in mathematical ecology is to study 
the growth and co-existence of species with diffusion in both homogeneous and patchy 
habitats. As noted before that diffusion, when it occxns, plays the role of increasiug sta- 
bility in a system of interacting populations [5, 7, 19, 20, 21, 23]. Levin [12, 13] has given 
elaborate survey of models with diffusion in both homogeneous and heterogeneous envi- 
ronment. McMurtrie [14], has also surveyed the literature related to populations model 
with diffusion and reported the effects of dispersal and spatial heterogeneity on stability of 
both single species and for Predator-Prey system. Nallaswamy and Shukla [17], considered 
a prey-predator model with functional response and diffusion and have shown, that if the 
equilibrium state is linearly stable, a subregion of the positive quadrant can be found in 
the phase plane where it is non-linearly stable with or without diffusion. 



5.1 Introduction 


130 


It may be noted here that in the above study the role of alternative or supplemen- 
tary resource on equilibrium levels of populations as well as on their stability has not been 
discussed, although the study of resource-based interacting population biology is an in- 
teresting area of research in population d3Tiamics. Some experimental investigations on 
micro-organisms using the chemostat [8, 18] have been conducted and perhaps the best 
laboratory idealization of nature for population studies has been described in [24]. Sev- 
eral mathematical models of such systems, involving competition and other types of non- 
interacting populations which depend upon growth limiting nutrient in a chemostat with 
constant input and variable washout rates have been studied [1, 10, 11]. Also some other 
mathematical investigations related to two competing populations which are wholly de- 
pendent on a self-renewable resource in a habitat without diffusion have been presented 
[6, 9, 15]. But very little attention has been given in the resource-based prey-predator sys- 
tem with diffusion, Freedman and Shukla [2], have studied the effect of a predator resource 
on a diffusive Predator-Prey system, showing the stabilizing role of diffusion. 

In this chapter, we, therefore, study a logisticaUy growing two species prey-predator 
type model with a self-renewable supplementary resource for prey population and diffusion 
in a two-patch habitats and discuss the stability of both the linear and nonlinear systems, 
under both the reservoir and no-flux boundary conditions. The model is proposed by 
keeping in view the depletion of forest resources due to increased forest resource dependent 
or independent industrialization and population causing patchiness in the Boon Valley, 
situated at the foot hills of Himalayas in India and affecting its biodiversity, Shukla et al. 
[ 22 ]. 

This chapter is organized as follows, flrst we write the prey-predator model with a 
self- renewable supplementary resource for the prey in a two-patch habitat. Next section, 
we study our main model in a two-patch habitat for both non-umform and uniform steady 
state cases under both reservoir and no-flux boundary conditions. 
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5.2 Mathematical Model 

We consider a dynamic model of two logistically growing animal (such as deer and wolf) 
species with prey-predator type interaction and diffusion in a two-patch forest habitat by 
assuming that the predator species uses the prey species as an alternative resource. In 
such a case the rate of change of density of the prey species decreases due increase in the 
density of the predator species, but its density increases due to the increase in the density 
of the prey species in both the patches. We also consider that there is a supplementary 
self-renewable resource for the prey species such that its density increases as the density 
of the resource biomass increases but the density of the resource biomass correspondingly 
decreases. Let Ri{s^ t), t) and yi(s, t), i = 1, 2 be the densities of resource biomass, prey 

and predator species in the i-th patch respectively. Assuming that the resource biomass is 
non-diffusing, the model governing the system can' be written as follows: 

^ = aiRi{l - - ctiRiXi (5.1) 

dx' 

yiPi(^i) d” dOiiRiXi H- (5.2) 

^ = yMVi) + 7iyiPi(a^i) + -^2.^ (5.3) 

0 < s < L 2 , and i = 1, 2. 

where the i-th patch is assumed to lie along the spatial length Li-i < s < Li (Lq = 0), 
Ci (i — 1,2) is the carrying capacity of the supplementary resource in the i-th patch and 
0 is the conversion rate of biomass density by the prey population assumed as constant. 
The functions gi(xi) and fi(yi) are the respective specific growth rates of prey and predator 
populations, Pi(xi) is the interaction rate (predator response functions) and Du, D^i are 
the diffusion coefficient of a:,- and yi in the i-th patch respectively. The constants oci, i = l,2 
are positive interaction rate coefficients of the prey species wdth the supplementary resource 
and 7 i, i = 1, 2 are conversion rate coefficients. 
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We assume the following assumption for gi{xi), fAyA, and 
AHi: gi(xi), fi(yi), pi(xi) GC2[0 ,oo) 

gi(0) > 0, fi(0) > 0, pi(0) = 0 

For Xf > 0 , g'i{xi) < 0, p;(xi) > 0 

For yi > 0 , f/(yi) < 0 

When the environment has a carrying capacity Ki and Mi respectively for prey and 
predator populations in the i-th patch, then gi{Ki) = 0, fi(Mi) = 0, for i = 1, 2. Further 
we assume that, 

AH 2 - 3 R*,x**, y** > 0, such that 
= Ci[ai - aixl*]/ai 

xrgiixr) - yrM^r) + o^iR^xr = o 

fi(yr) + 7zPi(a;r) = 0 

The model is studied under two set of boundary conditions i.e. reservoir and no-flux. 
In the case of reservoir boundary conditions, we take 


xi(o,i)=xr, 

X2{L2i t) — 3^2 

(5.4) 

yi(0,t) = yr. 

y2{L2,t) = yT 

(5.5) 

and in the case of no-flux boundary conditions, we consider, 


axi(0,t) _ Q 
ds 

_dX2{L2,t) 

ds 

(5.6) 

5yi(0,t) _ 
ds 

_ dy2{L2,t) 
ds 

(6.7) 


We also assume the continuity and flux matching conditions at the interface s = Li. 
The continuity conditions at the interface s = Li are, 

Xi(I/i,t) = X 2 (I'i,t), yi{Rut) = y2{Li,t) and Ri{Li,t) = R2{Li,t) (5.8) 


I 
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and the continuous flux matching conditions at the interface s = Li 

for Xi{s,t) and yi(s,t) 

are, 


P 3xl{Ly,t) aX2(Lut) 

” ds as 

(5.9) 

n ^ViiLut) dy 2 {Li,t) 

“ ds as 

(5.10) 

Finally the model is completed by assuming some positive initial distribution of each 

species, for z = 1, 2, that is. 


0) = Xi{s) > 0, Li_i <s<Li 

(5.11) 

yi{s, 0) = (ii(s) > 0, Li-i <$<Li 

(5.12) 

Ri{s, 0) = Roi{s) > 0, L,_i <s <Li. 

(5.13) 


5.3 Analysis of the Model in Two-Patch Habitat 


Our aim is to analyze the long time behavior of the system in both uniform and nonuniform 
cases in the following sections. 

In next two subsections we will study the model (5.1)-(5.13), in the case of nonuniform 
steady state and the uniform state state respectively. 

5.3.1 The Non-uniform Steady State: Under Both Sets of Bound- 
ciry Conditions 


Let Ui, Vi and Wi are the steady state solutions of the prey populations (x,), predator pop- 
ulations {ui) and the supplementary resource {Ri). Then the steady state system becomes: 
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Du 


(fui 

+ Uigi(«t) - Vij>i{ui) + 9aiWiUi = 0 


Ai 


dPvi 

ds^ 


+ + 7it;iPi(tii) = 0 


(5.15) 

(5.16) 


Now substituting the value of Wi from (5.14) into (5.15) and (5.16), we get, 


CpUi 

i ^^2 + ni^i(ni) - ViPiiui) = 0 

(5.17) 

d^Vi 

(5.18) 


where, 

C- 

Qi{ui) = gi(zii) + 9ai—{ai - cmui), i = 1, 2 

ai 


Since 5i(0) = gi(0) + Cidon > 0 and ^((uj) = gf(tii) - < 0, for i = 1 , 2 . Hence 

the behavior the steady state system (5.17) and (5.18) with same set of boundary and 
matching conditions are identical with the steady state solutions of the previous chapter, 
i.e. the case when there is no supplementary resource for the prey populations, since the 
assumption H 3 of the chapter 4, now modified as, 

3 x**,y** > 0 , such that, 

^rsi(xT)-vrPiK)=o, ( 5 . 19 ) 

fi(!(r)+7iPi(^r)=(l- (5.20) 

Also conditions (4.53) and (4.54) in this case modified as foiiows, 


{Xi - a:**)[a:i^i(a:i) - yiPi(xi)] < 0, \/xi ^ x,** and min{yl*,y;^} <yi< max{yr,yl*} 

( 5 . 21 ) 

ivi - yi*)[yMyi)+'YiyiPii^i)] < ^ yV and min{xl*,xl*} <Xi< max{x*j*,x;*} 

(5.22) 

The analysis in the remaining part of this chapter we assuming these two conditions. 
Remark: Since we are only interested on the positive steady state of the system. Therefore, 
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from (5.14), Ui < ai/oa and hence > gi[ui), V Ui. Now, from (5.19) and (5.20), 
get, 


we 


.. a;re?.-(rcr) 

Vi = —v iJv ■ ■ > vl anda:r>a:t, 


Pt(2:.-*) 

where x\ and y* satisfy the following equations of chapter-4 [see (4.3) and (4.4)], 


- yrPi(x*) = 0, (5.23) 

^iiVi) + 7iPi(a;i) = 0. (5.24) 

Hence in presence of a supplementary resource for the prey population, the level of steady 
state solutions of both the species are higher at each location in the habitat. It is shown 
numerically in Fig.5.1 [see section 5.4]. 


Now exactly in similar manner as the statement and proof of Theorem 4.4.4, 4.4.5 for 
reservoir boundary conditions we can state the following theorem without proof. 


Theorem 5.3.1 (i) There exists a positive, continuous, monotonic solution of system 

(5.17) with continuous flux at Li. 

(a) There exists a positive, continuous, monotonic solution of system 

(5.18) with continuous flux at Lx. 

Now we consider the stability analysis of the system (5.1)-(5.3), (5.8)-(5.13) with 
reservoir boundary conditions (5.4) and (5.5). First we state the local stability of the 
system by the following theorem. 

Theorem 5.3.2 The steady-state, continuous, monotonic solutions of the linearized sys- 
tem (5.1)-(5.S) with reservoir boundary conditions and continuous flux at the interface 
s = Lx is asymptotically stable provided the following conditions are satisfied: 


(z) A'i < 0, yi <0, Zi<o 


(5.25) 
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where, 


(u) Uf < 4Xiyi 

(Hi) XiyiZi < y.y^f + 


(5.26) 

(5.27) 


.s; - &(«i) + j,,. ^ ^ 

_ f-i 2w,-\ I 

i ai (1 j - aiUi, Ih = -[7it;.-p;.(ui) - p,(u,)], W, = 

for Xi < Ui < , Vi < Vi < y^* , where x^* and y’^* are given by ( 5 . 19 ) and ( 5 . 20 ), 

for i = 1, 2. 


Proof: Linearizing (5.1), (5.2) and (5.3) by using, 

Ri{s,t) = Wi{s)-i-ri{s,t) 
Xi[s, t) = Ui{s) + n,(s, t) 
yi{s, t) = Vi(s) + mi(s, t) 

we have. 


dri 

It 


= Tf 


■ . (. 2Wi\ 

0-i ^ j CXiUi 


- niOiWi 


dui 


- ^i[gi(^ii) + Uig'i{ui) - ■Uip.(ui) + OaiWi] - miPi{ui) + riOociUi + Du 
= mi[fi(ui) + Uif/(t;f) + 7iPi(tti)] + THJiVip'iiui) + 

os^ 


d^Ui 

1?' 


dt 


(5.28) 

(5.29) 

(5.30) 

(5.31) 

(5.32) 

(5.33) 


Using (5.28), (5.29) and (5.30) the corresponding boundary and matching conditions 
can be obtained as follows 


ni( 0 ,t) = 0 = n2{L2,t), mi{ 0 ,t) = 0 = 7712(^2, i) 


ni (Li ,t)=n 2 {Li ,t), mi {Li , t) = m2 (Li, t) 
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Now we consider the following positive definite function, 

DiflFerentiating (5.34)w.r.t. t, we get, 

VTi) = V f^' L 5r,-l 

By using (5.31), (5.32) and (5.33), we get, 

^ ? X._: ”*■ - ^»PiK) + e^iWilds 

^ + Vifl(vi) + 'YiPi{ui)]ds 


(6.34) 


+ E 


nimi[~pt(ui)+iiVip[{Ui)]ds + '^ f ' nir,a,leu, - wMs 

' 2 


T-v /* * (9 77-2 T~> d^TTl' 

? “ /i,_. "■ s 


Therefore, 


2 rii 

X/ ^ + -Zjr? + TUiTiimi + 2>ViTin,](is 


(5.35) 


where the functions Af, 3^^, Z/,-, and Wj are as follows, 

= S,i{ui) + «igK«t) - ■«,?■(«») + dotiWi, yi = fi(t;i) + t;if/(u,) + jiPi{ui), 

^ \ 1 Of- 

= ai\l- —j - OiUi, Hi = -[jiVip'iiui) - Pi(«i)], Wi = -jldui - Wi] 

Hence V is negative definite if the conditions (5.25), (5.26) and (5.27) stated in theorem 
are satisfied, for both i = 1,2. I 


In the following we state the nonlinear 


stability conditions for the system. 


Theorem 5.3.3 The steaiy-state, contmuous, monotonic oolnUons of nonlinear system 

(5.1}-(5.S), (5.S)-(5.1S) with resen,o.r boundary conditions (S.fHS.S) is asymptotically 
stable m the subregion ofR: {.;• < ,, ^ ^ 

the following conditions are satisfied: 


where 


(0 Ki <0, fifyi< 0 , < 0 

(n) J^i < 

KififyiKi < NyiMl, + 




A/;* 


Ki 

Ki 

•^wi 


Xigijxi) Uigi{ui) _ Pi(a;,-) - Pj(u^) 

Xi-Ui Xi-Ui 

yi^ijyi) (Vi ) 


yi-Vi 


+ 7iPi(«i) 


/- Ri + Wi\ 

Pi(^t) ~ Pt (Ui) 


im- 


1 
2 

^[Oui - iZij. 


Xi — Ui 


Pi(Wi) 


(5.36) 

(5.37) 

(5.38) 


(5.39) 


Proof: By using (5.28), (5.29) and (5.30), we get from (5.1), (5.2) and (5.3), 


0>i ^1 ^ J OjUiJ UiOCiRi 


drii 


= ni 


Xigiixi) - Uigi{ui) .. Pi(a:i) - Pi(ui) . ^ „■ r s. . 


(5.40) 

d^m 


drrii 

dt 


= rrii 


vMyi) - Viii{vi) 


Vi -Vi 


+ 7iPi(^i) 


+ ni 


ityi 


Pi{xi)-Pi{ni) 


(5.41) 

ff^nxi 

+ P2i-^ (5.42) 
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Here also we consider the same positive definite function as in the case of linear stability. 

nt) = tl^_^lln‘i + m1+rf]ds 

From which we get, 




By using (5.40), (5.41) and (5.42), we get, 

_ Xj Uj 

yMvi) - Vift(vi) 


ds 


m = v, /■ 


Xi — Ui 


ds 


m? 


I ''■^1— JL L 

, 2 ft Ri Wi^ 

^ A ,_, V ^ 

^ r^i r D.fx,') — D.-l'u.-'l 

1 Ai-i 


ds 


ds 


+ 


Pi{Xi) - Pi{Ui) ^ , 

d^rii 


Tiim, 


Jl ft.. «J 

^ d^Th' ^ rL* f^TD’ 


^ f-Li 

ds + y TiTiiailOui — i?i]d!s 

2 JLi^i 


Therefore, 


y(t) 


2 rU 

y I [f^xin-j + + J^zirj + 2Muinim,i -I- 2Mv,irini]ds 


where the functions Myu J^zU and U^i are given by (5.39). Hence V is negative 
definite if the conditions (5.36)-(5.38) holds, for both i = 1,2. 


It is noted that if we linearize the conditions of theorem 5.3.3, we get the conditions 


of theorem 5.3.2. 
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the 


Same results as diseased above for reservoir bounds^ conditions are also true for 
system (5.1)-(5.3). (5.8M5.13) with no-flmc boundary conditions (5.6) and (6.7). 


5.3.2 


The Uniform Equilibrium State: Under Both 
ary Conditions 


Sots of Bound- 


The main purpose of this section is to find the conditions for local and global stability 
of the uniform equilibrium state, i.e. xAs,t) = K'. ^ M' and fli(s.f) = C\ for 

0 < s < I 7 , t > 0, of the system, under both sets of boundary conditions. 

Theorem 5.3.4 The equilibrium {C-,K\M') is locally aeymptoticallp stable, if the fol- 
lowing conditions are satisfied, 


H* + doiiC* <0, i = 1, 2, 


and 


- Pi(ir*)]2 < 4[ir; + eaiC*][M*ti{M*)], i = i, 2 


where H* is given by 


H* = Zi{K*) + 


(5.44) 


(5.45) 

(5.46) 


Proof: Lineariazing the system (5.1)-(5.3), by using 


Ri{s,t) 

= C* + ri{s,t) 

(5.47) 

Xj (s, t) 

= K* + ni{s,t) 

(5.48) 

yi{s,t) 

= M* + mi{s, t) 

(5.49) 

II 

r aiC*^ 

(5.50) 


we get, 
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dui 

dt 


os^ 


os^ 

Considering the following positive definite function, 

“ 5 ? H, ~ + (yi~M'f + di(lu- c-)=] 

where dj, i = 1, 2 axe positive constants. 

Differentiating (5.53) and using (5.50)-(5.52), we get 

^ rU 


(5.52) 


(5.53) 


diaiC* 


Ci 


^ + ^o;iC'*]ds + ^^ nn^M*i'i{M*)ds + Y^ r. 

^ ^ rLi 2 

+ zl J + 'riM*p[{K*)]ds + Y^ I ' niri[ai{9K* - diC*}]ds 

^ f d'^rii ^ f^i d^rrii t 

Again, using integration by parts, for both type of boundary conditions, we get 


ds 


(5.54) 




ds. 


Choosing di, for t = 1, 2, such that, the coefficients of Uiri become zero, i.e. di = d 2 = 
9K* jC*. Therefore, from (5.54), V is negative definite, if the conditions (5.44) and (5.45) 
axe satisfied. I 

Moreover, 


Theorem 5.3.5 Let Rl-VQoLiC* > 0. Then the equilibrium {C*,K*,M*) is locally asymp- 
totically stable, if the conditions (5.45) along with 


HI + eoLiC* < Du 


TT^ 


{Li-Li.iY’ 


fori = 1,2, hold. 
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Habitat 


Proof: By using Poincare’s Inequality, we get 

n , _ _ 




Therefore from (5.54), we get, 


TT rLi 


2 rL, 

E / n]x[H* -f- OaiC* - Du ],. , A , 

Y A- ,2r 2 ^ - 

1 A._, *[ Ci + 

P‘L (!?)■'• 


Hence the theorem. 


We now state the following theorem for 


the global stability of the uniform steady state. 


Theorem 5.3.6 The uniform steady-state (C*,K%M*) is globally 


A(x,) = - M^Mx,) 

Xi - K* 


asymptotically stable if 


+ 6aiC*<0, ^Xi^K*. 




Xi - K* 


diM ^ ff.(M*)-f;(y,) 


Vi - M* 


Vi - M* 


Xi [ yi - M* 


2(a _ D \2l , 2L 




Ci Xi 


(5.55) 


(5.56) 


. (5.57) 


Proof: Let us consider the following positive definite function. 


Vi -M*- M*ln^ ds 
M*J 


(5.58) 
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Differentiating (5.58) with respect to t, and using (5.1)-(5.3) we get, 

4C(V)’-S- 


I JLx^i Xi 


+ E /. ivi - y*f 


1 


Xi-K* 
Vi - M* 


+ eaiC* 

2 -n/ 


ds 




7i 




Xi - K* 


ds 


+ E ^ “ ^*){yi - M*)ai \e - i?,-] ds 

1 JLi^i 


+ tOu2‘ «- 

I Xi OS j JLi^i yi 


ds^ 


ds 


Now consider the integrals, and using both set of boundary and flux matching condi- 
tions at the interface, 


,dxi 


dxi 


[xi(0,t) - = [X 2 (l 2 ,i) - = 0 


ds 


and 


we get, 


and 


[5,(0. () - M‘)^( 0 .t) = [n(h,t) - M‘]^{L2,t) = 0 


1 ' 


I'i Xi - K* d^Xj 
Xi ds^ 

yi-M*d\ 


ds 




i.- ir* 

ds , 


V 4-1 yi ds^ 1 4-1 y! 


1‘i M* (dyi 


ds 


ds 


(5.59) 


(5.60) 


(5.61) 


(5.62) 


\ds_ 

Now if the conditions (5.55) — > (5.57) are satisfied, then V{x, y) < 0, and ^((7*, K*, M*) — 
0. Therefore V{x, y) is negative definite over il > 0, x > 0, ^ > 0 with respect to Rl = C*, 
x\ = K*^y^ — M*, proving the theorem. ® 
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On comparing the analysis of this chapter with chapter-4 we may conclude that 

the role of supplementary resource is to increase the level of non-uniform Lady st tL 
each point of the habitat. ^ 


5.4 Numerical Example 


Now, we discuss a numerical example to study the behavoir of the steady state solutions 
of the above system. The results of this section 


predator system as discussed in chapter-4. For this 


are compared with the case of a prey- 


form of functions: 


we consider the following particular 


n (^1 _ 5 . (^1 - , P^(u.-) = e.Ui, i = l ,2 

also for simplicity take, Ci = C 2 = C. Then the steady state system (5.17) and (5.18), 


becomes 


n ^ 

+ Ui 

ds^ 


r. n _ ^ I r N 

1 ~ ) H (a,- — a, -It,) 


Ki 


Oi 


CjUfUj — 0 


d^Vi 


with reservoir boundary conditions 

^1 (0) = xj*, M 2 (i' 2 ) = X 2 *, and 
^i(0) = jir,^2(T2)=j/2**, 

and the continuity and flux matching conditions at the interface s = Li, 


dui 


dU2 


dvi 


dvo 


Dn^iU) = Dn'^iU), i^ 2 i^(T,) = D 22 ^(X 0 , and 


ds^~^'^ ds^~^' -^^ds 

y-i{Li) = U2 {Li), Vi{Li) = U 2 (Zri). 


(5.63) 

(5.64) 

(5.65) 

(5.66) 

(5.67) 

(5.68) 



Steady Slate Solutions 
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The equations (5.63)-(5.68) are numerically solved, for the Mowing set of parameters, 
T, = 10, Li = 20, Du = 0.8, = 0.9, = 0.8, = 0.9, r, = 0.03, rs = 0.025, 

ai = 1.0, 02 = 1.0, Si = 0.03, S 2 = 0.01, ^ = b.OE - 05, ej = 5.0£; - 05, = 100, 

K 2 = 125, Ml = 75, M 2 = 50, oi = 5.0£’ - 05, q ;2 = 2.0JE; - 05, C = 60, 71 = 0.4, 
72 = 1.0 and 9 = 0.7, and the solutions are shown in Fig. 5 . 1 . By using the above values 
of the parameters, we get, xf = 87.2 > xj = 80.77, x^* = 118.58 > x^ = 109.09, 
y” = 118.60 > yi = 115.385, and = 168.58 > = 159.09. It is noted from the 

figure that, in presence of a supplementary resource for the prey, the level of steady state 
distributions of both the prey and predator species are higher at each points of the habitat. 


5.5 Summary 

In chapter-4, we have discussed the existence and stability of the non-uniform steady state 
corresponding to a prey-predator system, where the predator species partially depends upon 
the prey species, in a two patch habitat with diffusion. In this chapter the same ecological 
problem is modelled and analyzed when there is a nondiffusing renewable resource for the 
prey population in both the patches. It has been assumed that both the prey and the 
predator population follow a general type of logistic growth in each of the patches as in the 
previous chapter. Also the self-renewable supplementary resource is logistically growing. 

It is shown, in a similar manner as in chapter-4, that there exists a positive, monotonic, 
continuous steady state solution with continuous matching at the interface for both the 
species. We have also obtained the criteria for asymptotic stability both linear and nonlinear 
cases under both set of boundary conditions. It has been further shown that in presence of 
supplementary resource the level of steady state distribution of prey population increases 
and hence the level of steady state distribution of predator species also increases. It is 
shown that the effect of patchiness is destabilizing even in presence of a self-renewable 
supplementary resource for the prey species. 
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Chapter 6 


A Two Competing Species Model 
with Diffusion in Homogeneous and 
Two-Patch Habitats 


6.1 Introduction 


The mathematical modelling of competing populations in a habitat has become quite rigor- 
ous after the work of Lotka and Vblterra. It is only during the last few decades that the envi- 
ronmental and ecological ('Ih'cts have been considered on the evolutionary processes of inter- 
acting speci(!K. rii(*se effects have been taken into account by identifying species mi gration 
with dispersive and convective processes which might have continuously varying properties 
due to environmental as well as ecological gradients in the habitat, and may depend upon 
densities as well as space co-ordinates of the species, [1, 2, 5, 9, 12, 13, 14, 18, 19, 20, 21]. 
In particular Kapur [13] has shown that, for a general class of competition models, diffu- 
sive instability cannot arise. It is known that the opportunities for migration and habitat 
diversification provided by the spatial changes in environment make coexistence of species 
possible which could not otherwise survive. It has been shown that an unstable equilibrium 
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may become stable with dispersion under certain conditions in a homogeneous habitat [9]. 

It may, however, be noted that the habitats in general, are heterogeneous due to vari- 
ety of factors such as geographical, climatic, ecological etc.. If the habitat is heterogeneous 
and favorable, different combinations of species are likely to be favored in some regions and 
maintained elsewhere principally by dispersal from more favored regions, and this process 
enhances the overall species richness. But opposite would be the case if the habitat be- 
comes unfavorable to the species due to climatic changes, resource depletion such in a forest 
habitat which might have been degraded due to industrialization leading to patchiness. 

It may be noted here that one of the methods for study the stability of interacting 
species model with space dependent properties is to divide the habitat in a number of 
patches with interaction and diffusion for each species is different but constant in each 
patch [3, 4, 8, 19]. The effect of patchiness on the stability of interacting population can 
also studied by using the same method [19]. Freedman et al., [4] studied a single species 
diffusion model by assuming that the habitat consists of two patches and shown that the 
steady-state solution is monotonically increasing and linear asymptotically stable under 
both reservoir and no-flux boundary conditions. Later they [3, 8] discussed the single 
species diffusion model in multi-patch environment and also studied the global stability for 
the system. In the above studies the effect of patchiness on two interacting populations 
has not been investigated. 

Keeping the above in view, in this chapter, we study the coexistence of two com- 
peting species with logistic type growths and diffusion in both homogeneous and the two- 
patch habitats in the similar manner as in case of prey-predator system already studied in 
chapter-4. This chapter is organized as follows, in section 6.2, the complete model of two 
competing species in a two-patch habitat is written. In section 6.3, we study the model 
in a homogeneous habitat without and with diSusion and compare the stability behavior 
of both the system. In section 6.4, we study our main model ( described in section 6.2 ) 
in a two-patch habitat for both uniform and non-uniform steady state cases under both 
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reservoir and no-flux boundary conditions. In the last section s ■ 

so.ut.ons 01 steady state solutions with „ bolClirr ““ “ 


6.2 Mathematical Model 


we consider a genera, nrode. for two conrpeting and disusing species in a two-patchy 

a . a . et x. (s. f) and s,.(s. t) be the densities of two competing species in the i-th patch, 

or t , . The system ,s then governed by the foliowing autonomous partial differential 
equations: 


( 6 . 1 ) 


Oxi(^s t'j 

= a;i(s, t)gi {xi{s, t)) - yi(s, t)p,. t)) + 

ds^ 

dyi{s,t) 2 

= Viis, t)fi {yi{s, t)) - Xi{s, t)qi {yi(s, £)) + (6.2) 


0 < s < L2 i = l,2 


where the i-th patch is assumed to lie along the spatial length Li_i <s< Li (Lq = 0), 

gi(Xi) and are the respective specific growth rates, p,(xi), afo,) are interaction rates 

of Xi(s,() and J(,(s,t), and Du and Bs. are the diffusion coefficient of x. and ft in the i-th 
patch respectively. 


We assume the following assumption for g,(xi), pi(xi) and qi(y,-) 


■^ 1 - St(^t)) ^'(^1)1 Pi{^i)i G C^[ 0 , 00) 

gi(0) > 0, fi(0) > 0, Pi(0) = 0, qi(0) = 0 
For Xj > 0 , g'i{xi) < 0, Pi(xi) > 0 
For Pi>0, f/(yi) < 0, q((j/i) > 0 
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When the first and second 
the i-th patch, then 


species have the carrying capacities Ki, Mi respectively in 


H2-- gi{Ki) = 0, fi{Mi) = 0, i = 1, 2. 

We also assume that, 


H 3 : 3 X* > Q, yl > 0 such that 


- vM^i) = 0 (6.3) 

y* ft {y* ) - Xi q,-(?/;) = 0 . ( 6 . 4 ) 


There are only two possibility of existence of positive equilibria (a:t, y^) (see also 
Fig. 6.1 and Fig.6.2 in next section), either 


(i) Ki > 


m 

qi(o) 


, Mi> 


g»-(o) 

p;-(o) 


or (ii) Ki < 


m 


, Mi< 


gf(0) 

pm 


Since from (6.3), we have, x*igi{x^) = 2/*Pi(a:J'), this implies that gi(x|) > 0 = gi(Ki) 
and again since from assumption Hi, gi(x,) < 0, for all > 0. Hence xl < Ki. Similarly 
from (6.4) and Hi, we get j/,* < Mi. 


We also assume the continuity and flux matching conditions at the interface s = Li. 
The continuity conditions at the interface s = la for this system are, 


Xi{Li,t) = X 2 {Li,t) and yi{Li,t) = y 2 {Li,t) (6.5) 

and the continuous flux matching conditions at the interface s = Li for Xi(s,t) and yi{s,t) 
are written as. 


dxiiLi,t) dx2{Li,t) 

^ = as 


( 6 . 6 ) 

(6.7) 
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The model is studied under two set of boundary conditions i.e. reservoir and no-flux. 
In the case of reservoir boundary conditions, we take 

X: (0, t)=x\ , X 2 {L 2 , t) = xl, and i/i(0, t) = i/* , i/ 2 (T 2 , t) = (6.8) 

In the case of no-flux boundary conditions , we have 

axi(o,t) _ dm(o,t) ate(i,,t) 

as ds = ' as 

Finally the model is completed by assuming some positive initial distribution, that 
is, 

2;i(5, 0) = Xt(5) > 0, Xi_i <s<Li, (6.10) 

yi{s, 0) = (5i(s) > 0, Li-i <s <Li. (6.11) 

We first study the existence and stability behavior of the system (6.1) and (6.2) in a 
homogeneous habitat without patchiness, the effect of patchiness will be investigated later. 


6.3 Analysis of the Model in a Homogeneous Habitat 

6.3.1 Model without Diffusion 

In this case Xi = x^ yi = y, gi{xi) — g{^), fi(yi) = f(y)> Pt(^i) = P(2^)) *=1)2. Thus the 
system (6.1) and (6.2) reduces to the following form: 

^ = a:g(a;) - yv{x) (6-12) 

at 

^ = yf(y) - ^^(y) 

at 


(6.13) 
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The functions g(c:), f (,) and p(a;) satisfy the same fype of assumption as mentioned 
m Hi and Ha- And the assumption Hs can be written as, 3 i’ > 0 and j/* > 0 such that, 

2;*g(2;*)-J^*p(a;*) = 0 (6.14) 

y*W-x*q{y*) = 0 ( 6 . 15 ) 

From equations (6.14) and (6.15) it follows that there are four equilibria, namely (i) 
Eo = [0,0], (ii) Ek = [K,0l (iii) Em = [0,M] and (iv) E* = The existence of E*, 

the interior equilibrium can be ensure by (see Fig.6.1 and Fig.6.2), either 


Case (i) K > 




Note: Now we note from the two isoclines given by (6.12) and (6.13) i.e., 


y = ^andx=^, 
Pix) q{y) 


that the gradients at {x*, y*) are respectively, given as follows 

^-^andf^) =5W, 

Now from Fig. 6.1, we have 


qfe-) ^ . sM ^ 


FT,* 


p(x*)’ 


-Ff| " p(a;»)' 


Multiplying both side by — Fr 2 P( 2 :*)) we get p(a:*)q(y*) > provived iTj < 0 and 

H 2 < 0. But from Fig.6.2, by similar manner, we get 


q(r) 

m 


< 


HI 


p{x*) 


=> p(a:*)q(y*) < 


Therefore the conditions p(x*)q(j/*) > corresponds to case (i) [Fig.6.1] above and 

other condition corresponds to case (ii) [Fig.6.2]. i 


Now, the variational matrix, in the general case, is given by 
g(x) + xg'(x) - yp'{x) -p{^) 

-q(y) f(y) + 2/f'(j/)-2rq'(2/) 


[M] 


(6.16) 





Figure 6.1: Existence of E*, i.e. {x*,y*) [for case (i)] 


Model in a Homogen..nc 



Figure 6.2: Existence of E*, i.e. {x*, y*) [for case (ii)] 


6.3 Analysi s of the Model in a Homn^eneom 


Now calculating the variational matrix from 16.161 for i -.w . 
hypotheses Hi -y Hj and the standard stabUitv th f T 

we note the following obvious remarks: ly o “r mary differential equations, 

The equilibrium point E, is always unstable. The equilibrium point is a stable or 

saddle point according as f (0) - f.q-(o) , negative or positive. Similarly the^^ ^ 

or saddle point according as gfO') - TWnVni io .t- ^ 

^ ^ ^ ^ ^ general, there is no 

“O®* mteresting non-tero equilibria 

-D 5 u It exists. 


Therefore, our aim of this section is to obtain local stabihty 
stability conditions for E*. 


and as well as global 


We now state main results of this section in the forms of Theorems 6.3.1, and 6.3.2 
and Lemma 6.3.1. 


Theorem 6.3.1 The equilibrium E’ is locally asymptotically stable, if the following con- 
ditions are satisfied, 


and 


H*<o, h;<o, 


(6.17) 


where 


P(^*W) < HiH* 


(6.18) 


Hi = + g(x*) - 


o^*g{x*) 

p{x*) 


P'i^l 


and HI = y*f'{y*) + f(y*) - 


y*i{y*) , 

^{y*) ^ 


iy*) 


■ 


Proof: By considering the following positive definite function 

V (x, y) = \[{^- + (2/ - y*?] 

the theorem can be proved.[ see proof of theorem 4.3.1]. 
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As (iisrusM'd \>v{nr(‘ jt rnay noted 

i.e. coriTspoiKis tu Honce E* ic ('i)’ 


is locally stable 

The final result uf this section, Theorem 6.3.2 ei, 


case (ii). 


ally aaymp<o(i,-ally Mai,),.. First we prove a lero !"! “ “““ 

prove a lemraa whrch establishes a region of attraction 


of the system. 


Lemma 6.3.1 Thr irqton of attmrfinr, n r . 

mracUon for all solut,ors mittafej i„ tho positive guadrant 


-'t {(l,j/) : 0 < I <ii-,o <S(< 


M} 


Proof: Itoiii ((i.l.’l .imi {f..I.3) and hypotheses {Hf] to (Hf), we get, 
m< x(0g(a:(t)) < g(o)a;(i)[i _ 

m < v{t)f{y(t)) < mvm-v{t)iM\, 

then, as t -4 cx), 

for :r(()) < Ah Urn sup x{t) < K, and for x(0) > K, ' lim sup x(t) ^ K, 

for ?/(()) < M, Itm sup y{t) < M, and for y(0) > M, lim sup y{t) -> M. 

Therefore, we m)t(« that .solutions initiating on 5Af\int %+ enter into int A. Hence A is 
the region of attraction of the given system. I 

Theorem 6.3.2 Let the conditions 

{x - x*)[xg(x) - ?/*p(x)] <0, Vx^x*, y ^ y* (6.19) 

{y - (j/) - ^*q(j/)] < O, V ar ^ x*, ^ y* (6.20) 
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'p{^) q(2/)‘ 

^ 4 

< — 

’2;g(x) - y*p(a:)1 

'yf(y)-a:*q(y)l 

_ X y _ 

xy 

X — X* 

y-y* 


p(a;)q(y) < 

’^g(2:)-y*p(x)l 

yf(y)-a;*q(y)' 


X —X* 

y-y* 


hold. Then E* (m F%g.6.2) which is locally stable is also globally asymptotically stable in 


Proof: By considering the following positive definite function 


y)= (x-x*- 3:*ln^^ + (y -y* -.y*ln^^ 
the theorem can be proved.[ see proof of theorem 4.3.3]. 


Remark: If we linearize (6.19) and (6.20), around the equilibrium point E*, we get < 0 
and i? 2 * < 0. 

In the following example we show that the conditions for the local and the global 
stability are feasible. 

Example: We take, 

g(x) = r (^1 - , fiy) = ^ ’ P(^) = = “2?/ 

Then the the equilibrium point E*{x*,y*) is given by 

- 5-^(^ ~ 0^1-^) * _ ~ 5^2-^) (6.22) 

rs — olxOi-xKM' ts — 01XO12KM 

From (6.22) it is clear that positive equilibria exists under two set of conditions, namely 
Case 1: r < aiM, s < a 2 K and rs < aia 2 KM 
Case 2: r > aiM, s > 0 : 2 ^ ’'■s > oiia 2 KM 

The nature of x* and y* with respect to various parameters can be calculated from (6.22). 
It is shown in the following table. 
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Table 6.3: Effects of varinnc 


x(i) 

rt 


OLIJ 

icieis on 

Case 1 
Case 2 

(uj 

st 


2/.- 4- 

Case 1 





Case 2 

(111) 

Kt 

xn 

Pit 

Case 1 




pU 

Case 2 

(iv) 

M t 

x^f 

1 2/U 

rCase 1 




prt 

Case 2 

(v) 

Ql t 

x*t 

Pii 

Case 1 




p;t 

Case 2 

(vi) 

^2 t 


p,*t 

Case 1 



xff 

p;t 

Case 2 


and y* 


< 0 Vx 2 :*, y*. 


Now we note that , = -rx'/K < 0, HI = ~sy-/M < 0 

teW - rp(.t)l ^ 

iC y-y* 

Therefore the condition (6.19) and (6.20) are automatically satisfied. Further both the 
local and global stability condition (6.18) and (6.21) becomes, 

TS 


<^ 10:2 < 


KM' 


which is automatically satisfied in Ca^e 2. Hence E* is locally as well as globaUy asymp- 
totically stable in Case 2. 


6.3.2 Model with Diffusion 

Now, we wish to consider the model (6.12) and (6.13) with diffusion and analyze the uniform 
equilibrium B* under both reservoir and no-fiux boundary conditions. The model (6.12) 
and (6.13) in this case can be written as. 


dx / \ n 

— =xg(x)-yp(x) + Ei— 

(6.23) 

^ = yf(y)-xq(y)] + E2^ 

(6.24) 
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0<s<L 

where A > 0 are diffusion coefficients. The reservoir and no-flux boundary conditions are 

respectively, 

a:(0, t) = x* = x{L, t) and y{0, t)=y* = y(L, t), (6.25) 

^ = 0 = ^and?fi) = o = ^. (e.2a) 

Now if we take the following positive definite function 

= / V{x,y)ds, 

Jo 

around the equilibrium point E*, then using similar arguments as in section 3.1, we get the 
following theorems for local and global stability of the system. 

Theorem 6.3.3 The equilibrium E* is locally asymptotically stable, if the conditions (6.17) 
and (6.18) are satisfied. 


Moreover by using Poincare’s inequality, it is pointed out that even if (6.18) does not hold 
true, this uniform equilibrium may become stable with diffusion under a condition as given 
in the following theorem. 

Theorem 6.3.4 The equilibrium E* is locally asymptotically stable, if the condition (6.17), 
along with 

p(x-)qte-) < (a- - ( h ; - ) (6.27) 

are satisfied. ® 


Theorem 6.3.5 If the conditions (6.19)-(6.21) are satisfied. Then E* is globally asymp- 
totically stable in A. 
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Thus the uniform equilibrium ET (in Fig.6.2) of the system is locaUy ns well as globally 
asymptotically stable under the same set of conditions, as in the previous case. Farther by 

Theorem 6.3.4, an equilibrium which is unstable without diffusion can become stable with 
diffusion. 


6.4 Analysis of the Model with Diffusion in a Two- 
Patch Habitat 

In the first subsection we study the model (6.1)- (6.11), in the case of uniform equilibrium 
state. In the next section we study the non-uniform case. 


6.4.1 The Uniform Equilibrium State under Both Sets of Bound- 
ary Conditions 

The aim of this section is to show that uniform steady-state is globally asymptotically 
stable. In this case, it is clear that, under both sets of boundary conditions, there is a 
uniform steady-state, x{s,t) = K*, 0 < s < 1 / 2 ) t > 0, y{s,t) = M*, 0 < s < L 2 , t > 0, 
(i=l,2) where K* and M* are the common uniform equilibrium point of first species and 
second species respectively in the two patches. Let x\= X 2 ~ K* and yl = y 2 = M*. By 
using similar arguments as in section 3.1, the following theorems for stability can be proved 
for the system (6.1) — )• (6.7) with (6.8) or (6.9). 

Theorem 6.4.1 The unique uniform steady state solutions (K*, M*) is locally asymptot- 
ically stable, if Hf < 0, fori,j = 1,2, and the following conditions are satisfied, 

PiiK*)qi{M*) < Hi*HT i=l,2 (6.28) 
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where 






Similarly as in Theorem 6.3.4. the above nniform steady state is locally asymptoti- 
cally stable under less stringent conditions as the foUowing theorem, 


Theorem 6.4.2 The ™i,ue nniform steady state solutions (K>, M’) is locally asymptot- 
.cally stable, .f H- - < 0, for i,j = 1 , 2, and (he follonhng conditions 


Pi(A'*)qi(Ar) < ( H\* - Du— — I lw*--n \ . 

I L,.3 ) , > = 1.2 (6.29) 

are satisfied. 


We now state the global stability of the uniform steady state, 


Theorem 6.4.3 The uniform steady-state solutions (Kf M') is globally asymptotically 
Stable, if the following conditions are satisfied, 


{Xi - ir*)[a;igi(a;i) - M*pi(xi)] <0, Vxif x*, 
iVi - M*)[yifi{yi) - K*qi{yi)] <0, V yi ^ M*, 
Pii^ihiiyi) < 



kf.(yi)-ii:*qi(y.-)l 

Xi-K* 

Vi -M* 


(6.30) 

(6.31) 

(6.32) 


Proof: Let V(x,y) be the positive definite function about x = K*, y = M*, given by 
V{t)=V(x,y) = jzffi-IC-irbai^f^ds 
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Differentiating with respect to t, and using (6.1) and (6.2) we get, 



Now by assumption (6.33), first integral of the right hand side is negative, V except 

* 1,2. Similarly by (6.34) the third integral of right hand side is negative, V yi, 

except yi = M*, i = 1,2. 

Now consider the integral. 


Under both set of boundary conditions, 


ds^ 


ds 


[x.(0, t) - K‘] 1) - K-] = 0 


ds 


Similarly, 


I2 




I J Li^l 


D2i\ 


ds < 0 


Vi j 5s2 

Hence V{x,y) < 0, and V(K*,M*) = 0. Therefore V{x,y) is negative definite over a: > 0, 
y > 0 with respect to x = K*,y = M*, proving the theorem. I 


6.4.2 The Non-Uniform Equilibrium State 

We assume that the unsteady state solutions of the model (6.1) to (6.11) exists, [see [18] 
(Theorem 1, page 111 and Theorem 3, page 123)]. 
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Our aim here, therefore, is to show that there exists a positive, monotonic, continuous 
steady state solution for the each of the species separately, under continuous flux matching 
at the interface of these patches, in cases of both reservoir and no-flux boundary conditions 
and derive the conditions for asymptotic stability for both linear and non-linear models. 


There are four possible cases: 

(1) x *2 > xl and > yl 

(2) x *2 > xl and y^ < y* 

(3) X 2 < xl and y^ < y* 

(4) xl < xl and y| > y* 

Without loss of generality, (3) can be reduced to (1) and (4) can be reduced to (2). 
Therefore we consider the case (1) and (2) separately under two set of boundary conditions. 

We assume the following conditions on Xi and yj for i = 1, 2. 

{xi - a;*)(r;,gi(a:i) - yiPi(xi)] < 0, Vr, # xl and Tnin.{yl,yl} < Vi < max.{yl, y^} (6.33) 

ivi - yi)[yMyi) - 2:zqi(yi)] < o, Vy^ ^ yl and min.{xl,xl} <Xi< max.{xl,xl} (6.34) 


The alalysis in the remaining part of this chapter is valid under these two conditions, 
(6.33) and (6.34). 

6.4.3 The Model under Reservoir Boundary Conditions: When 

X2 > and yl > yl- 

We first consider the steady-state problem and show that there exists a non-uniform positive 
monotonic solutions Ui{s), Vi{s) under reservoir boundary conditions and with continuous 
flux at the interface. 

The steady-state system takes the form 

+ Uigiiui) - ViPiiui) = 0 


(6.35) 
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^ dPvi{s) 

■^2» ^^2 + ^ift(Ut) - Uiqi{Vi) = 0 

(6.36) 

The reservoir boundary conditions are 


ui{Q)=^xlu2iL2) = x; 

(6.37) 

^l(0)=yi,U2(l2)=P2 

(6.38) 

The continuous solutions and flux matching conditions at the interface are 


^i(Ti) = ^ 2 (^ 1 ); ui(Li) = V2{Li) 

(6.39) 

dui{Li) du2{Li) 

-^11 t = JLJ 12 , 

as as 

(6.40) 

n dvi{Li) dv2(Li) 

^21 y = 1^22 y 

as as 

(6.41) 


Let pi(.s, ai) and f/i(s,/?i),0 < 5 < Li be the unique solutions of (6.35) and (6.36) respec- 
tively, for i=l , such that 


(6.42) 

- = A. ?i(0,A) = ».' (6.43) 

Let P 2 {s,a 2 ) and q 2 {s,P 2 ),Li < s < L 2 he the unique solutions of (6.35) and (6.36) 
respectively, for i=2, such that 

- = (X2, P2{L2, 0 ^ 2 ) = (6.44) 

os 

S^h lM. = fe, 5,(I,,A) = y| (6.45) 

OS 

The existence of the monotonic solutions follows, if we can show that there exist cxi, 012 , Pi 
and p2 such that 

Pi(Ti,Q!i) =P2(Ti,Q;2) , gpLi, Pi) = <l2{Li, P2) 

^ dpiiLuOi) ^ dp2{Li,a2) = 

=M2 


(6.46) 

(6.47) 
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From (6.35), moltiplymg both side by 2iu,lds and integrating w.r.to s from 0 if i=l 
and from L 2 if i=2 , we get, 


dui 

ds 


2 ^ 


(s) 


(6.48) 


Similarly from (6.36), we have. 


dvi _ 2 2 ni{s) 

_ds\ Nfi(vi)-Wiqi(ni)]dt;i(s) 


(6.49) 


In order to construct our required solution we need some preliminary lemmas. We 
prove these lemmas are similar as in chapter-4. 


Lemma 6.4.1 If ai > 0 , then 


dpi{s,Oii) 

ds 


> ai on 0 < s < Li. 


Proof; From (6.48), we get. 


api(s,o;i) 


l2 


ds 


= ar 


_2_ 

Ai 


rpi{s,ai) 

/ [Pigi(Pi) -9iPi(Pi)]dpi 


(6.50) 


since 

= Oil > 0 , Pi(0,ai) = , 

Then there exists Sj > 0 such that pi{s, oi) > on 0 < s < Si. If not, let sq, 0 < sq < Li, 
be the first positive value, if it exists, such that pi(so,a:i) = xj. Then by the mean value 
theorem there exists s such that 0 < s < So and dpi{s, ai)/ ds = 0; that is. 



2 r‘^"’“'^[p^gj(p^)_gipi(pi)]dpi. 

Dll Jxi 


(6.51) 


Nov since pi(s,ai) > xj for 0 < s < s, and from (6.33), we have, pigi(pi) - giPi(pi) < 
ce the right hand side of (6.51) is negative , giving a contradiction. Therefore 
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function of pi, the lemma follows. 


notingthat(dpi(s,ai)M)^i 


IS an increasing 

■ 


Lemma 6.4.2 If 0 < p^ < and Q2 > 0, then 

dP2{s,OC2) 

>0^2, Li<S<L2. 


Proof: From (6.48), we have, 


9p2{s,a2) 

ds 


_ 2 2 rPiisAi) 


From ( 6 . 33 ). we have, P2S2M - g2P2{P2) > 0 . Hence 
0 < P2 < 3^2- 


92P2(p2)]dp2 

's > 0(2, Li < s < L 2 for 


Similarly from (6.49) and (6.34) we have. 


Lemma 6.4.3 // A > 0 , then 

dqiis,/3i) 

ds 


> 01 on 0 < s < Li. 


Lemma 6.4.4 7/0 < 92 < 2/2 and 02 > 0, then 

592(5, 02) 


ds 


> 02 i Li ^ s < L 2 . 
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There exist four continuous functions F- fi i-i 9'\ t 

for proof see same lemma in chapter 2. 


Lemma 6.4.5 Define F,,{a,) by Fu{c,) = p,(L„a,). Then there e»t> 4 > 0 such that 

Fn ■ [0,Q:i] [a:*,x;] 

■^12 : [ 0 ,^ 2 ] -> [xj.Xi] 


Lemma 6.4.6 Define F2iil3i) by F 2 i{fii) = qfihJi). Then there exists h > 0 such that 

Fii : [Oji] -4- 

Fi 2 '■ [0,^2] [2^2, 


Theorem 6.4.4 There exists a positive, continuous, monotonic solution of the steady state 
system (6.35) with continuous flux at Li- 


Proof: From lemmas 6.4.1 and 6.4.2, it follows that any solution, we construct, must be 
monotonic. By Lemma 6.4.5, for each 0 < 02 < 0:2, we can find an ai such that 0 < cti < di 
for which Pi(Li, ai) = P 2 {F\,cx 2 ). Hence ai can be solved as a function of a 2 , oi = h{oi 2 ), 
to give a continuous solution of (6.35) with (6.37), (6.39) and (6.40). 

Let 


^ api(Li,%2)) p dp2{Li,a2) 

0(02) D,2 . 


Clearly 0(03) is continuous on 0 < 02 < 02. Then we have 


0(0) = Dn 


dpi[Li, ai) 


> 0 , 


ds 
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and 


G(d2) = 


ds 


< 0 . 


Therefore 3 ^ 2 , 0 < < da, such that G(d 2 ) = 0. Hence 


the theorem. 


Similar results are also true for the second species. 


Theorem 6.4.5 There exists a continuous, monotonic solution of the system (6.36) with 
continuous flux at Li . 


We now study the asymptotic stability of the steady state system (6.35) and (6 36 to 
(6.41). 


Theorem 6.4.6 The steady state, continuous, monotonic solutions of system (6.35) and 
(6.36) with continuous flux at the interface s = L, is asymptotically stable provided the 
following conditions are satisfied: 


W ^te(wi)]<o, 

d 

^ Kf2('yi)] < 0 > 

\lfor xl<Ui< x^. Ml <V 2 < M 2 and i = 1, 2. 


(lit) [Pt(^i) "ir <5i(Ui)] 4 


A. 

dui 


[Wigi(^i.)] 


■roof: Let the steady-state solution of system (6.35) be 


u(s) = 


wi(s), 0 < s < Li 

ti2(s), Li ^ S 1/2 
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And let the steady-state solution of system (6.36) be 

v(s) = l 0<5<Li 

I Ms), Li<s<L2 

Lineariazing (6.1) and (6.2) by using , 

Xi(s, t) = 'Ui(s) -f ni(s, t) 

yi{s, t) = Vi{s) -f m,(s,t) 

We have, 


171 


(6.52) 

(6.53) 


dni{s, t) 

Wt 


[gi M) + «igi(«i)] - niUfp;.(uf) - mipj(ui) + (6 

ds^ ^ 


.54) 


5m,(s, t) ^ 

~ mi [fi(vi) + Vif'(vi)] - miUiq^(vi) - niqi(vi) + (6.55) 

Using (6.52) and (6.53) the corresponding initial , boundary and matching conditions 
be obtained as follows 


can 


Ms, 0) = Xi(s) - Ui{s), mi{s, 0) = 5i{s) - t;i(s) 
ni(0,t) = 0 = n 2 (L 2 ,^) , mi(0,t) = 0 = m 2 (i: 2 ,t) 

MMi) = n2{Li,t) , mi{Li,t) = m2(Ai,t) 




ds 




Now we consider the following positive definite function, 


n<) = EX,)_5N+">f]* 


(6.56) 


Prom which we get, 


£, fLi 


drii drrii 



6.4 Analysis of the Model with Diffusion in a TwrvPat.r}, Habitat 


172 


By using (6.54) and (6.55) , we get, 


V{t) = 


^ [St(wi) + Wjg^(uj)] ds [ '>^iViPi{ui)ds 

- E / ’ riimipi{ui)ds - Y, D 2 i [^] ds 
1 T [ ds 

^ fl'i 2 

+ E / ds-Y r rn]ui(^{vi)ds 

1 J I JLi-i 


2 rLi 

Y I nimiqi{vi)ds - Y Ai 

1 1 JLi_i 


[ dnj 

ds 


t 2 


ds 


Since Pj(ui) > 0, 'dxi ^ 0. Therefore second integral of the right hand side of above equa- 
tion is negative. Similarly sixth integral is also negative. Hence V{t) is negative definite if 
conditions (i), (ii) and (iii) are satisfies, and the theorem is proved. I 


Theorem 6.4.7 The steady-state , continuous , monotonic solutions of non-linear system 
(6.35) to (6.4-1) with continuous flux at the interface s = Li is asymptotically stable in the 
subregion of R: {a:* < Xi, Ui < Xj, yl < yu Vi < y^, for i = 1, 2}, provided the following 
conditions are satisfied: 


(1) ~ ^ Q 


Xi — Ui 


(2) < Q ^ 


Vi - Vi 


(3) XiPi 


■uipi(xi) ^ niqi(yi) 

2 

'xigi(Xi) - UigiiUi)' 

’yMyi) - VifiiviY 

a;i Vi . 

< AuiVi 

Xi—Ui 

yi- Vi 


Proof: Prom (6.1), (6.2), (6.52) and (6.53), we have 
^ = {ui -f ni)gi{ui + ni) - Uigiiui) - {vi + Tn,)pi(tti + n,-) -t- ViPi{ui) 4- Du 


d'^rii 

ds^ 


(6.57) 
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dt 




Now consider the positive definite function, 


ds^ 


From which we get, 

After using equations {6.57) and (6.58), we get 

?£/■ 


yn) 


(uj 4 - ni)gi{ui + m) — UiP,i(u{) 


'>i™tP,(»t + t,.)* + EA, F‘ rj^ds 

1 1 -/ii-i os^ 


+ 


X r^i 


(vj + + ruj) - Vifi(vi) 


ds 


Zii 


ds 


Now, 


U. 

tlx 

y j nimiqi{vi + mi)ds + yD2i rrn^—ds 

Pi(Ui + Tti] - Pi(7i,-) 


' qi(u.- + mf) -qi(^j) ' 

rrii 

‘ • > 
2 


rii 


TOi) - qi{vi) 


mi 


>0, Vm. 


> 0, V m,-. 


Also by using (6.56), we get 

2 

I 

1 


2 ,1 


Similarly 


(as. 


■A 7-. 


ds <0 
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Therefore l>(f) < 0 if the conditions (1),(2),(3) holds true. Hence the theorem. | 

Example: Let us consider an example by choosing the following function in our model 

(6.1) and (6.2), 


g^{x^) ru (l , f,(j/.) _ ^.2. (l - , pi(a:i) = cuXu qi(y,) 

0 < Cii, C2i <1, i = l,2 


C2iyi- 


Then 


k; = and jlf- - 

rnr,, - cue, KM ruTu-CuccK,M, 

From the above it is clear that K- < Ki and M* < Mi. 

Here the conditions (1) and (2) of non-linear stability becomes, 

Now the function Gi{xi) = Xigi{xi) is decreasing from the point x^i = Ki/2. Similarly 
the function Fi{yi) = VifiiVi) is starts decreasing from t/oi = Mi/2. 

Now if we choose ru, r 2 i, cu, C 2 i are such that > Ki/2 and yl > Mi/2. Then (1) 
and (2) are satisfied in the following region: 


< Xi, Ui < xl, yl < yu Vi < yl, for i = 1, 2} 


It is pointed out here that the condition (3) of non linear stability can also satisfied 
in above region for a suitable choice of gi, f,- and Pi satisfying Hi ^ H3. In this particular 
case the condition (3), becomes 


x:iyi [ciiUi - C2iVif < 4riir2iUiUi 


Xi -f Ui 


Ki J 

1 Mi J 


(6.59) 


This c! 
maxirr 


idition is automatically satisfied if for a particular choice of rit, r 2 i, Cu and C2t, the 
.j p of the left hand side is < minimum of right hand side. i.e. 


xxlyl 


[cux*2 - C2iy*if < (^ - 1 ) ( 


Mi 


-1 


(6.60) 
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It is clear from (6.60) that for suitable choice of r,, and the above condition can 

be satisfied. In particular if we take 
rii = 0-34, ri2 = 0.35, r2i = 0.574, r22 = 0.49, 

Ki = 340. A's = 175, = 287, M 2 = 490,’ 

cii = 0.0003, C 12 = 0.0003, C 21 = 0.0001 C 22 = 0.0001, 

we get xj = 250.15, = 100, y* = 299.5, y| = 500. 

In this numerical example (6.60) is automatically satisfied. 


6.4.4 The Model under Reservoir Boundary Conditions: When 

x *2 > ."tJ and y{ > yl 

For X 2 > the rnonotonicity of steady state solution is similar as in the above case and for 
Mi > M 2 , the rnonotonicity just reverse of the above case. Hence the steady state solution 
is monotonic under both set of boundary conditions. 

In this case the linear asymptotic stability conditions for the steady state system 
under reservoir boundary conditions are similar as the above case. 


6.4.5 The Model under No-Flux Boundary Conditions: When 

xl > Xi and y^ > yl 


In this section, we study the steady state model with no-flux bormdary conditions, 

dui{0) _ Q _ du 2 {L 2 ) _ dvijO) _ p _ dv 2 {L 2 ) 
ds ~ ds ' ds ds 

As before, we assume X 2 > arj and y 2 >y\- 


i.e., 

(6.61) 


Let Pi(s,q 0 and gi(s, A), 0 < s < A be the unique solutions of (6.35) and (6.36) respec- 
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tively, for i = 1, such that 


-0 , Pi(0,ai)-Q:i 

-^-^ = 0 , gi(0,A) = A 


(6.62) 


(6.63) 


Let P 2 (s,Q^ 2 ) and g 2 (s,/?), Li < s < L 2 he the unique solutions of (6.35) and (6.36) 
respectively, for i = 2, such that 


dp2(L2,a'2) ^ ,, , 

^ = U , P2(L2,a2) = 02 

dQ2(L2,p2) . , 

-0 , 9i(L2,/%) = /% 


(6.64) 


(6.65) 


The exist, ence of the monotonic solutions follows, if we can show that there exists oi, 02 , 
and p 2 f^ach that 


Pi(Li,q:i) =P2(Li,a2) , gi(Li,^i) = g2(Li,/?2) 

(6.66) 

dpi(Li,ai) ^ dp2(Li,a2) 

r. = ^12 0 

as os 

(6.67) 

^ dqi(Li,/3i) ^ dq2(Li,l32) 

ds as 

(6.68) 


From (6.35), multiplying both side by 2dui/ds and integrating w.r.to s from 0 if i=l and 
from L 2 if i=2 , we get 

= [ [tti&(iii) - t;iPi(tii)]dui(s). 

JJli Jai 


dui 


ds 


Similarly from (6.36), we have 


dvi 


(6.69) 


ds 


2 /-"iW 


D2i 


-[ [viii(vi) - Uiqi(vi)]diii(s). 

i dPi 


(6.70) 


As in the case of reservoir boundary conditions, here also similar type of lemmas hold, as 


stated below : 
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Lemma 6.4.7 If a} > x], then 


dpi{s,ai) 

ds 


> 0 , 0 < s < Li 


Proof: From (6.69), we get, 


dpi{s,oci) 


t2 


ds 


rpiis,oii) 

bigi(Pi)-?iPi(Pi)]dpi. (6.71) 


2 /-PiCs-ai) 
Dll •'cti 


Since dpi{Q,ai)/ds = 0, pi(0,Q:i) = ai there exists si > 0 such that Pi(s,q:i) > oi on 
0 < s < .sj. If not, let sq, 0 < sq .f'li be the first positive value , if it exists, such that 
Pi(soi oi) = «i- Then by the mean value theorem there exists s such that 0 < s < Sq and 
dpi{s,ai)/i)s — 0; that is. 


0 = _ / [pigi(Pi) - 9iPi(Pi)] dpi (6.72) 

Uli Joi 

But the right-hand side of (6.72) is non-zero,since pi > cti on 0 < s < s. Also oi > x*, 
therefore, pi > x\. This implies that, Pigi(pi) - QiPiipi) < 0, giving a contradiction. 
Hence pi(s, Qi) > Oi, for all value 0 < s < Li. This implies, p(s,o;i) > ij, on 0 < s < Li. 
Hence by (6.33) and (6.71), the result follows. ■ 


Similarly by using (6.70) and (6.34) the following lemma follows: 


Lemma 6.4.8 If 0i > y*, then 

M^>0,0<s<L, 

ds 


Lemma 6.4.9 J/O < P2(5,a2) < thena2 < x*, implies p2{s,a2) < 0^2, forU < s < L2. 



6.4 Analysis of the M odel with Diffusion 


ip_a Two-Patch Habitat 


P.oof: Sin,.,. „ < ^ _ 

s.„ce U.™ ^ ^O.^N„w^ 


Again by using (6.70) and (6.34) the following 


lemma follows: 


Lemma 6.4.10 Ffi) < q2{s,P2) < y^, thenp, 


'2 < implies q^{s, /%) < for <s<Lo 


Lemma 6.4.11 Define Guioa) by Gu{ai) = Pi{Li,ai). [ 
that 

: [a:t,Q;i] -)• [xl,x*2], 
Gi 2 : [Q;2,a:2j -> 


Then there exists a,i > 0, such 


Lemma 6.4.12 Define G^iifii) by G2i(A) = QiiLi, Pi). Then there exists Pi > 0, such that 


^21 -• [yi*,/3i] H- [yl,yl] 

G22 ■• [4,^2] [yhvl] 


Proof: Same as Lemma 2.3.6 in chapter 2. 


Theorem 6.4.8 (i) There exists a continuous, monotonic solution of system ( 6 . 35 ) with 
continuous flux at Li. 



(n) There exists a continuous, 


continuous flux at Lj . 


monotonic solution of system (6.36) with 


Proof: Analogous to theorem 2.3.5 in chapter 2 


Finally, in a similar manner as in the case of t>io rooo • v 
, „ . , . . n tne case of the reservoir boundary conditions, the 

following linear and nonlinear stability theorems can be proved. 


Theorem 6.4.9 The .Urndystate, continuous, monotonic solutions of system (6 S5) and 

16.S6) u,m conUnuous flus, at the .nter/ace s = L, is asymptotically stale provided the 
following amditwns arc satisfied: 

(?) [?i.g,(ui)] < 0 , 

^ < 0 , and 


dvt 

(Hi) [p,(w,) + q,(t;.-)f <4 


dui 


«igi(wi)] 


/«'■ < «. < .T-, y; <Vi< ^ 2 * and i = 1, 2. 


dvi 




Theorem 6.4.10 The steady-state , continuous , monotonic solutions of non-linear system 
(6.35), (5.36), (6.51) and (6.61), with continuous flux at the interface s = Li is asymptot- 
ically stable in the subregion o/R: 

{a^i < «t < ^2> Vi ^ Vi ^ ^2) fo^ ^ = 1, 2} 
provided the following conditions are satisfied: 

- UigijUi) ^ Q 


( 1 ) 
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( 2 ) < 0 


?A 


( 3 ) x,]h 


’ Vi 


l2 


< 4iUiVi 


^igijXi) - UigiiUj) 
Xi - Ui 


yMui) - Vifiivi) 


Vi -Vi 


Remark: The above theorems imply that the system will settle down to a steady state 
distribution for each species, in the two patches under certain conditions, the magnitude of 
the steady state distributions of both the species being lower than it’s initial density. The 
above analysis also suggest that patchiness destabilizes the system which we can note by 
comparing the results with the uniform steady state case, [see section 6.4.1]. 


6.4.6 Both The Species have Uniform Steady State In The Sec- 
ond Patch: 

In this case, t he steady-state solutions of both the species are non-uniform in the first patch 
and constant in second patch, i.e. U 2 = x^ and V 2 = M 2 , Li < s < L 2 , t>0. As shown in 
the general case we note that here also the steady state solution is positive, continuous and 
monotonic in the first patch. For stability analysis we use the positive definite function, 



+ f"\yi-vifds+ [ {y2-V2fds (6.73) 

Jo •'Ll 

and analyze* in t he similar manner as in Theorem 6.4.3 and Theorem 6.4.7, and can prove 
the following theorem. 


Theorem 6.4.11 Let U 2 = x^ and V 2 = Then the steady-state, continuous, monotonic 
solutions of non-linear system (6.1) to (6.1) and with either (6.8) or (6.9)and with continu- 
ous flux at the interface s = Li is asymptotically stable in a subregion ofH: {Ki <xi < K 2 , 
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X- < «1 1 xr,. A/j 1 //5 . T’l < Kh }, provided the following conditions are satisfied: ' 


, x,gj(.r!) - nig](ui) ^ 


Ji - Ui 



1 ! 

<0 , 


n,ipal:r,r ^ r--i 

Xlgl(Xl)-Uigi(ui)‘ 

2/ifi(2/i) -‘yifi(ni)' 

X] Vi 

Xi-Ui 

yi - 


Similar as the prey-predator system, the conditions of the above theorem 6.4.11 are 
k«vs in xmniber then liie conditions of non-linear stability in the general nonunifonn case in 
nvo-pat.rh habitat and more in number then the uniform case. 


6.5 Numt^rical Examples 


Now we show tlu? existence of positive, monotonic, continuous steady state solutions of the 
system with reservoir boundary condition and flux matching condition at the interface by 
considering the following example corresponding to systems (6.35) and (6.36), . 


= (Oii - biiUii - CiiU2i) Uii 

(6.74) 

£?lt ^2 ^ *“ b2iU2i C2iUii) U 21 

(6.75) 

The reservoir boundary conditions are: 


till (0) = Xj, Ui2{L2) = X 2 , 

(6.76) 

«2i(0) = 2 / 1 ) ^ 22 (^ 2 ) — 2/2 

(6.77) 

The cont inuous flux matching conditions at the interface are 

duu(Li) p dui 2 {Li) 

Dn ^ dx ' 

(6.78) 



■igurc' 

,* »«* *>' 

G.3: The st< 

a n r. n 

mdy state 

solutions, when - 

= 162.162 > = 55.882, yl = 141.892 > 

— o. 

an = 

0.115 

an 

= 0.825 

021 = 0.135 

022 = 

= 0.3 

hi ~ 

0.002 

h2 

= 0.005 

hi = 0.004 

hz = 

= 0.002 

Cij = 

0.0001 

Cyj 

= 0.0001 

C21 = 0.0001 

022 = 

= 0.0001 


= 0.75 

Dr. 

! = 0.8 

D 21 = 0.75 

Dn ■ 

= 0.8 





Figtirc* 6.4: stc*ady state solutionSj when a;? = 100.0 < xX = 149.99, yo = 74.999 < 

y* = 149.99 
an = 0.275 

a 12 = 0.825 

021 = 0.63 

022 = 0.3 

bn = 0.002 

hu = 0.005 

621 = 0.004 

622 = 0.002 

cn = 0.0005 

^32 = 0.001 

C21 = 0.0003 

C22 — 0.001 

Dn = 1.2 

Aa = 1.2 

Ai = 1.2 

D 22 = 1.2 





6.5 Niitiicrirai {•'.x.-unph-s 




I igur(‘ hilrris of diffusion coefficients in steady state solutions, when xt = 55 882 < 
- 127.72.S. ^ 32.353 < = 113.614 

Oil -- 0 01 15 0j2 = 0.0625 O21 = 0.0135 022 = 0.024 

i>n = 0.0002 5 i2 = 0.0005 621 = 0.0004 622 = 0.0002 

Cn = 1 .DOE - 05 Ci^ = 1.00£? - 05 C21 = 1.00£: - 05 022 = 1.00.£^ - 05 



fi.G Smiitnary 



o„ , p fr22ai) 

ax 


dx 


(6.79) 


'‘n(il)-U,2(il), >121(1.1) =U22(i,), (gggj 

W. l,av,. s„Iv,..l sy.t,.,„ „f .quations (6.74) to (6.80) aumerically osing finite difference 
metlKKl in vnrunii. ..lesee for different set of parameters statisffng the stability conditions 

and tin* c'<vrrc;^p(tn(iinR graphs are shown in Fig.6.3, Fig.6.4 and Fig.6.5. 


Ft™,, , 1... almve Bra,.h.s it can be noted that in aU the cases the steady state solution is 
positive, continuous and monotonic. In particular from Fig.6.3 and 6.4, we see that for all 
possible roinlnnatmn of .r^s and y^’s, the steady state solutions are monotonic. In Fig.6.5, 
the effects of diffusion rocfficienfs on the non-uniform steady state distributions are shown 
and it is noted tlmt as diffusion coefficients becomes vary large in both the patches, the 
steady .sf.ate solutions tends towards linear distributions. 


6.6 Siinimary 


In this chapter, we have considered a general competition model of two species with diffusion 
in a two pat <'h environment. The model has been given by a system of two autonomous 
partial differ«*nlial etjuations. We have obtained the existence of the positive, monotonic, 
eontimnni.s .stea<iy-stat<* solution for each species, with continuous flux at the interface. 
The rondif ions for asymptotic stability both linear and non-linear model cases under both 
re.s(‘rvoir or ncfflux boundary conditions have been derived. It has been shown that in both 
cases t he noii-uniform steady-state is asymptotically stable under appropriate conditions. 
It is also shown t haf. the uniform steady-state of the system is globally asymptotically stable 
under less stringent conditions. It is further noted that the patchiness destabilizes the non- 
uniform steady state which can be seen by comparing the results ( stability conditions 
) with corresponding to uniform steady state, as in the case of prey-predator system in 
chapter-4. 
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Chapter 7 


A Competing Species Model with 
Diffusion and a Common 
Supplementary Resource for Both 
the Species in a Two-Patch Habitat 

7.1 Introduction 

As discussed in chapters 3 and 5, the study of resource-based competition model is an 
important area of research in population biology. Some experimental investigations on 
micro-organisms using the chemostat [9, 17] have been conducted in fifties and perhaps 
the best laboratory idealization of nature for population studies has been described in [22]. 
Several mathematical models of such systems, involving competition and other types of 
non-interacting populations which depend upon growth limiting nutrient in a chemostat 
with constant input and variable washout rates have been studied [1, 11, 13]. Also a few 
other mathematical investigations related to two competing populations which are wholly 
dependent on a self-renewable resource in a habitat without diffusion have been conducted 
[7, 10, 16]. 
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It may be noted here, that though the effect of diffusion on the growth and co-existence 
of interacting species have been discussed by several investigators [4, 15, 18, 19, 20, 21], very 
little attention has been paid to study such problems with alternative or supplementary 
resource and diffusion. Recently Freedman and Shukla [4] have studied the effect of an 
alternative resource for predator on a prey-predator system with diffusion in a homogeneous 
habitat. 

It may be noted that in the above studies, the effect of patchiness on population 
coexistence with diffusion in each patch has not been considered. Though, Freedman et al, 
[2, 3, 6] studied a single species diffusion model by assuming that the habitat consists of 
two patches and shown that the steady-state solutions are monotonic and the asymptotic 
stability conditions under both reservoir and no-flux boundary conditions but they did 
study the resource-based population model. They have neither considered the effect of 
patchiness on two interacting populations, nor the role of supplementary resource in patchy 
habitat. 

Keeping in view the above and noting the discussion in chapter-6, in this chapter, we 
study the effect of a common supplementary resource on the coexistence of two competing 
species with diffusion, in a two patch habitat. This chapter is organized as follows, first we 
write the complete model of two competing species with a common supplementary resource 
in a two-patch habitat in section 7.2. In the next section, we study our main model in a 
two-patch habitat for both non-uniform and uniform steady state cases under both reservoir 
and no-flux boundary conditions. 


7.2 Mathematical Model 


We consider a general model of two competing species diffusing between two homogeneous 
patches in a given habitat as discussed in the chapter-6. We consider further that there 
is a common non-diffusing self-renewable resource for both the species in the two patches. 



7.2 Mathematical Model 



Let Ri (s, t), Xi{s, t) and y.(5, t) [i = 1, 2), be the resource biomass density and populations 
densities of two species competing for this resource, such that each species logisticaUy grows 
in the absence of the other species, and the rate of growth of each species decreases due to 
the presence of the other species in the i-th patch. We assume that the density of each of 
the competing populations increases with the increase in the density of the supplementary 
resource biomass and corresponding resource biomass density decreases. The system, is 
then, governed by the following autonomous partial differential equations: 

dRi Ri, 

- aiRi{l - - aiRiXi - 0iRiyi ( 7 . 1 ) 

dxi , , , , 

— = Xigi{xi) - yiPi(xi) + eaiRiXi + Du-— (7.2) 

^ = yMvi) - mtivi) + (PPiRiVi + (7.3) 

0 < 5 < 1 / 2 , and i = 1, 2. 

where the f-th patch is assumed to lie along the spatial length Li_i < s < Li (Lq = 
0), Ci, i = 1,2 is the carrying capacity of the supplementary resource in the i-th patch 
and the constants Oj and Pi, i = 1,2, are positive interaction rate coefficients of the 
alternative resource with the first (xi) and second (yj) species respectively in the i-th patch. 
The functions gi(a:j) and fi(yi) are the respective specific growth rates, Pt(a:i), qt(yi) are 
interaction rates of Xi{s, t) and yi{s, t), and Du, D^i are the diffusion coefficients of Xi and 
yi in the i-th patch respectively. The constants 6 and (p are the conversion rate coefficients 
for the first (aii) and second (y,) species respectively. 


We assume the following assumption for gi(ri), fi(yi), Pt(a:i) and qt(yt) 
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When the first (xi) and second (y,) species have carrying capacity Ki, Mi respectively 
in the i-th patch, then gi(iiQ) = 0, {{{Mi) = 0, i = 1, 2. Further we assume that, 


AH 2 ’. 3 R*, x^*, y^* > 0, such that 
Rl = Ci[ai - aix’l* - A2/,**]/ai 

a:rgi(xr) - yrPt(^**) + OaiR^x** = 0 

yruvr) - ^r^iivD + 4>myr = o 

The model is studied under two set of boundary conditions i.e. reservoir and no-flux. 
In the case of reservoir boundary conditions, we take 

?/i(0,t) = yr>?/2(T2,i) = 2/2** 

and in the case of no-flux boundary conditions, we have 

9a;i(0,t) _ dx2{L2,t) 

ds ds ' 

dyii0,t) _ dy 2 {L 2 ,t) 

ds ds 

We also assume the continuity and flux matching conditions at the interface s = Li. 
The continuity conditions at the interface s = Li for this system are, 

Xi{Li,t) = X 2 {Li,t), yi{Li,t) =y 2 {Li,t) and Ri{Li,t) — R 2 {Li,t) (7.8) 


(7.4) 

(7.5) 

(7.6) 

(7.7) 


The continuous flux matching conditions at the interface s = L\ for Xt(s, t) and yi{s, t) 
are written as. 


^ dxx{Lx,t) ^ dx2{Li,t) 

g— 

^ dvi{Lut) ^ dv2{Li,t) 

= g_ 


(7.9) 


( 7 . 10 ) 
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Finally, the model is completed by assuming some positive initial distribution of each 
species, for i = 1,2, that is, 


0) = Xi(s) > 0, Li_i <s< Li 

yi{s, 0) = 5i{s) > 0, Li_i <s< Li 
Ri{s, 0) = Roi{s) > 0, Li_i <s< Li. 


(7.11) 

(7.12) 

(7.13) 


Our main aim is to study the long time behavior of the uniform and nonuniform 
steady state solution of the above model under both tj^e of boundary conditions (namely, 
reservoir and no flux) and the continuous flux matching conditions at the interface of the 
two patches. 


7.3 Analysis of the Model in a Two Patch Habitat 

7.3.1 The Non-uniform Steady State: Under Both Sets of Bound- 
ary Conditions 


Let Wi{s), Ui{s) and Vi{s) be the steady state solutions of the supplementary resource 
[Ri), first species (xj) and second species (t/i). Then the steady state system of the above 
proposed model becomes, 


c 

Wi = — [a,- - aiUi - 

Oi 


d?Ui 


+ Uigi(ui) - ViPi(ui) + BoiWiUi = 0 

rftn)’ 

+ ^^ifi(^^t) - ■“iqi(^t) + <i>0iWiVi = 0 


(7.14) 

(7.15) 

(7.16) 


Now substituting the value of Wi from (7.14) into (7.15) and (7.16), we get, 

ds^ 


(7.17) 



where, 


(7.18) 


ds^ ' ~ = 0 


Giiui) 

Vi{ui) 

HVi) 

Qi{vi) 


C- 

Ei{ui) + Ooii—{ai — ajitj), 
di 

Pi(wt) + 9aiPi—Ui, 

Ot 

ft(wi) 4- ({>Pi~{ai — PiVi), 

di 

C- 

+ (j>oiiPi—Vi. 
ai 


Again the functions Qiiui) and are logistic type, since 
^t(O) = gt(0) + GaiCi > 0, g((ui) = g'-(ui) - 9a]Ci/ai < 0 , Vttj ^ 0 
^.(0) = f.(0) 4 <l>P,Ci > 0, Tlivi) = fl{vi) - (l>PfCi/ai < 0, V^.- ^ 0. 

Also, since 

VM = 0, Vliur) < 0 \fui ^0 and 
Qi(0) = 0, Q'(r;i)<0Voi7^0, 

therefore the functions Vi{ui) and Qi(vi) are interaction type functions as considered in 
chapter-6. Further, the assumption H3 of the chapter 6, is now modified as, 

3 Xi*,y** > 0 , such that, 

^r^i(<*)-yrw*) = o, (7.19) 

yrHyr)-^VQiiyn = o, ( 7 . 20 ) 

We also assume as in chapter- 4 , 5 , 6 , the following conditions, 

{Xi - x;*)[xigi{xi) - yiPi(xi)] < 0 , \fxi ^ X** and min{yl*,y;*} < yi < max{y** ,yl*} 

(7.21) 

— ^ ^ ^*** ^ < max{xj^ ,X2 } (7.22) 


is valid throughout the chapter. 
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Then the system (7.17) and (7.18) with boundary conditions, 

ui (0, t) = x**, U2{L2, t) = xl* and Ui(0, t) = yl\ V2{L2, t) = y^* 

or 

-Q_ ^^2(-^2,t) dvi{ 0 ,t) dV2iL2,t) 

ds ’ ds as ’ 

and the continuity and flux matching conditions at the interface s = Li, 


(7.23) 

(7.24) . 


ui{Li,t) — u2{Li,t), vi{Li,t) = V2{Li,t) 


_ du-^ 
ds 

become exactly of the same type of steady state system, as described in chapter-6. 


= Du‘-^(L,). and 


(7.25) 

(7.26) 


Remark: Since we are only interested in the positive steady state of the system, therefore, 

Ui < aj/a, and u, < a^/ Pi and hence 
Qiiui) > gi{ui), V Ui and Pi(ui) > Pi(ui), V m, 

■^x(w.) > f,(n,.), V Ui and Qi(ui) > qi(ui), V Uj 

Thus, in presence of the non-diffusing, self-renewable supplementary resource, the 
steady state problem becomes the same type of steady state system of two competing 
species as in chapter-6, with the modified growth rates Qiiui) and ^i(vi) and the modified 
interaction rates 'Pi(uj) and Qt(^^t)- Therefore, the behavior of the steady state solutions 
will be same as in chapter-6. Now, exactly in similar manner, we can state and prove the 
following theorem, 


Theorem 7.3.1 fi) There exists a positive, continuous, monotonic steady state solution 

Ui for the Xi species xvith continuous flux at Li. 

(a) There exists a positive, continuous, monotonic steady state solution 

Vi for the yi species with continuous flux at L\. 

Now we go for the linear and nonlinear stability of the system, which wiU be different 
from the case studied in chapter-6. 
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Habitat 


Theorem 7.3.2 Th. steady sMe. continuous, monotonic 


( 7. 1)- ( 7. 2) with reservoir boundary and 


solutions of the linearized system 


continuous flux matching conditions at the interface 


s = L. is- asymptotically stable prouided the follomny condUions are satisfied: 

(0 < 0, 3^,- < 0, < 0 

(m) uf < 4 Xiyi 

(m) XiyiZi + TUiViWi < XiVj + yiWf + ZiUf. 

where 


(7.27) 

(7.28) 

(7.29) 


6t('?^t) + tttgj(ui) — Vip'flui) + 9aiWi, 
cr _ „ A 2 u;A 


^ ^ j OLiUi 

Vj = ^[(j)Vi - lOi], and 


yi - fi(t'i) + Vif'flvi) - UiCfflVi) + 
^i = -2lPi(ui) + qi{vi)], 


W,==^[9u,-wi] 


(7.30) 


Proof: Linearizing (7.1), (7.2) and (7.3) by using, 


Ri{s,t) = Wi{s) + rflsfl) 

Xi{s, t) = Ui{s) + n,(s, t) 
yi{s,t) = Vi{s)+mi{s,t) 


we have. 


dri r f 2wi\ 1 

-^ = Ti j^Oi (^1 - — j - aiUi - fliVij - niaiWi - mifliWi 


dfi' 

= ’^i[gi(^‘t) + Ui^flui) - Vip'flui) + 9 ociWi] - mig>i{ui) + rfloiUi + Du-— 

dtri" (9^772* 

- mi[fi(t>i) + Vifl{vi) - Uiq(-(ui) + (}>fliWi] - niqi(i;i) + r^fliVi + 


(7.31) 

(7.32) 

(7.33) 


(7.34) 


(7.35) 


(7.36) 
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using (7.31). ,7.32) and (7,33), the corresponding boundary and hatching conditions 
can be obtained as follows, * tuitions 

ni(0,t) = 0 = n2(L2,i), mi(0,i) = 0 = m2(L2,t) ( Reservoir case ) 


-^(0.i) = ^(L„t) (No-flux case) 

ni(Lj,t) = 0 = n2(Li,f), mi(Xi,i) = o = m2(Ri,t) 

, . , . .. dm 

ds 




ds - (js ' Qs 

Now, we consider the following positive definite function, 

r^i 1 


'"W = Er k? + n? + mads 

1 ‘'-t.i-.i z 


From w’hich we get, 

^ r-t.- 

I 

By using (7.35), (7.36) and (7.34), we get. 




dri dui dmi 


ds 


'■■(*) = E r r} 

1 

^ ^z^[gi(^i) + W:g,-(ui) - i;,-p((u,) + 9 aiWi]ds 

+ Xr /, ”^,'[fi(^^i)+u,f;-(ui)--Uiq((ui)-f 

1 '''^*“1 

2 rL( 2 L. 

+ 2^ nimi[-pi{ui)-qi{vi)]ds + J2 ninaiieui - Wi]ds 

I j JLi-i 


1 1 210, A 

fli ( 1 - j - QiUi - piVi 


ds 


(7.37) 


ds 


^ rLi 

I [£irf -f J^iuf + yimj + 2W,Tin,- -f TUiUimi + 2V,mir,-]ds 


Therefore, 

1/(0 = 
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where the functions y„ Z,, U,. V, and W. are given in (7.30). Hence V is negative 
definite if the conditions (7.27) (7.29) hold. 

We now show that there exist a small subregion of 

: {rnm{xlx;} < Xi,Ui < < y^^Vi < max{y*,y*}}, 

where the system is nonlinearly stable. For this let, 


II 

- Uigi(Ui) ^ p,(Xi) - Pi(ui) , , ^ 

r- — ii- .r » d-paiitj. 

(7.38) 


. yMVi) - vMvi) ^ qi{yi) - qi{vi) , ^ ^ ^ 

■ Vi-V, y,-v, 

(7.39) 


- Oi (l - ^ ) Oitti - PiVi, 

(7.40) 


■Mit = -^[p.-(ui) + qi(t;i)], 

(7.41) 


Ki = y[M- - i?i], 

(7.42) 


Ki^j[0ui-Ri]. 

(7.43) 


then the following nonlinear stability theorem is proved. 

Theorem 7.3.3 The positive continuous steady state solutions of the system is nonlinearly 
asymptotically stable, for i = 1, 2, if, 

(i) < 0, Myi < 0, < 0 (7.44) 

(fiz) ffxiJ^yi-hfzi d" ^ d“ ("f-4:6) 


Proof: By using (7.32), (7.33) and (7.31), we get from (7.1), (7.2) and (7.3), 


dn 

dt 


Ti 




maiRi - mi^iRi 


(7.47) 
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- iniPi(ui) 4- TiOaiUi + Di^-^ 

ds^ 
(7.48) 

Vi ~’^i ‘ _ y. ^^iRi — niqi(vi) + ri4)^iVi + Pj ,- — ~ 

J ds^ 

Her. also, we consider the same positive definite Action, as in the case of linear staW^ 


5n, 

dt 

drrii 

dt 


TTlj 


- Uigiiuj) _ p.(a;.) _ 

yMVi) - _ qiivi) - qi(vi) 


no = E£J[r? 


From which, we get, 




T 

By using (7.47), (7.48) and (7.49), we get 
V{t) ^ 


ctt dt ' dt 


ds 


+ 


E i; r. 

1 JL,^1 


My^l ^ifi(^i) ^ qiiVi) - qi{vi) , , ^ ^ 

Xj— f- (ppitq 

Vi - Vi 


Vi -Vi 


2 


ds 

ds 


rLi 




1 

2 ri, 


I ^ Ri + Wi' 

I I — J - OLiUi - PiVi \ ds 


rL>i 2 

+ XJ J rat7ni[-pi(ui) - qi(T;i)]ds + E / * 

1 ^*-1 I ‘'I'i-l 

■*■ ? C, ’■‘"‘“‘1^“' - + 1 Ou £_ + E D. 


d^TUi 

~d^ 


••"I J 

'^i'~Tr^ds 


Therefore;, 
V{fy = 


X r^i 


J2 L i^xirii + Myiml + ACif? + iN’^imnii + 2Kimiri + 2AC£r,ni ] ds 

i J Li^i 




( dirii ' 


ds 


-1 \ ds J — I ds 

where the functions A4t, A^i, A4i, Muh J^vi and Mv,i are given by (7.38) -> (7.43). Hence 
V is negative definite if the conditions (7.44) -> (7.46) hold true, for both i = 1,2. 


Now, we consider the uniform steady state case, 
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7.3.2 


The Uniform Steady State: Under Both 

Conditions 


SGts of Boundary 


Our mam purpose of this section is to find the conditions for local and global stabUity of the 
uniform steady state, i.e. JJ.-(s,f) h C. x,(s.t) h if- and y,(s,t) s M-. tor 0 < s < L,. 
t > 0, of the system, under both sets of boundary conditions. 


Theorem 7.3.4 The equilibrium {C*, K\ M*) is locally asymptotically stable, for i = l,2, 
if the following conditions are satisfied, 


(^) 

H* + eaiC* < 0 



(7.50) 


f: + < 0 



(7.51) 

{in) 

[Pi(A"*) + qi(Af*)]2 < 4[F; + eaiC*]lF* + i = 1, 2 

(7.52) 

(iv) 

[H: + 0aiC*][F* + cl>piC*] 

aid)M*' 

a 

<[F* + <j>0iC*][9K*-d>M*? + 




1 

(7.53) 


where 


H- = gj (K-) + - M’pr{K‘) 


and F; = fj(Af) + - K%(M‘) 


(7.54) 


Proof: Lineariazing the system (7.1)-(7.3), by using 


dri 

'm 


Ri{s,t) = C* + ri{s,t) 

(7.55) 

Xi{s,t) = K* + ni{s,t) 

(7.56) 

yi{s,t) = M* + mi{s,t) 

(7.57) 

Gi . 

(7.58) 


we get, 
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dt 

drut 

Ih 




)+neaiK'+Du^ (7.59) 


Consider the following positive definite function 


(7.60) 




+ {Vi - M*f + di{Ri - C'*)2] 


(7:61) 


where rf,, ?' - 1, 2 are positive constants. 
Differentiating (7.61) and using (7.58)-(7.60), we get 


V 


2 

E /,;_ n?K + fo.c-j& + j: „e|g. ^ 

T [''' A^2S£L] , ", 

v4., 4 Ci ^ + qi{M* 


^ ? £. + E £_ - d,C*Jds 




(7.62) 


Integrating by parts and using the matching conditions at the interface and for both set of 
boundary conditions, we get 

?£. = -P'‘C i^h 




d^nii 2 

-as = 


f drui ' 


ds. 


(7.64) 


Choosing for i — 1,2, such that, the coefficients of rriiri become zero, i.e. di = = 

dM‘/C\ Therefore, from (7.62), V is negative definite, if the conditions (7.50), (7.51), 
(7.52) and (7.53) are satisfied. I 

Moreover, 


Theorem 7.3.5 Let + 9aiC* > 0 and/or + (f)piC* > 0. 'Then the equilibrium 
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(C-,K\M-) is locally asymptotically stable, if the conditions (7.52) and (7.5S) along vMh 


H: + BaiC* < Du 


^ (17^ 




(7.65) 


for i — 1,2, hold. 


Proof: By using Poincare’s Inequality, we get 


D, 


I * ~ ds<Du-, — 


ds 


and 


rU 

D2i / 

JLi-i 


drui 


ds 


ds ^ D2i 


TT^ 


{Li-L, 


— r 

,-l)2 JLi-i 


rUids 


Therefore from (7.62), using (7.63) and (7.64) and choosing di = (j)M*/C\ for i = 1,2, 
get, 


we 




'Zt n^,[«t + 9aiC--Du- 

1 ( 


TT^ 


]d3 


'{Li-Li-ir 

+ zC' tnna’+d0,cr-Dsi ,^ J d 

2 rLi 


j J Li^i L/j j ./Lj— 1 


1 r^i 


+ ^ / UinailBK* - (j)M*]ds 
1 


Hence the theorem. 


We now state the global stability of the uniform steady state. 


Theorem 7.3.6 The uniform steady-state (C*,K\M*) is globally asymptotically stable if 

^ ^ Xigi(xi) - M*Pi(x,j ^ 
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Bi 


Vi-M* 


+ 

Vi 


<4 


AjBj 

^iVi 


and 




a 


Pl(3:i) ^ qt(yt) 


n2 




Vi 


(7.67) 

(7.68) 

(7.69) 


Proof: Let us consider the following positive definite function, 

r^i / 'T- \ ^ fLi 


V{x,y,R)= %Il {vi - 

i£(^ 


+ 


C* - C*ln 


C* 


(7.70) 


Differentiating (7.70) with respect to t, and using (7.1) and (7.2), we get, 


VM = E/" 

2 rL 

■ U 

A 


Xi-K*\dxi , (yi-M*\dyi , (Ik-C*\dRi 


Xi 


dt 


+ 


Vi 


dt 


+ 


K (X, - K*f 

i=l 


Xi 


Xigi{xi) - M*Pi{xi) 


^ fLi (y- - M*f [yift(yi) - K*qi{yi) 


i=l 

rLi 


Vi 


Ci\ 


ds 


Xi-K* 

yi-M* 

2 tli 


+ 6aiC* 


B-mc* 


Ri 

ds 

ds 


dt 


ds 


- E/' W-C-) 

+ E r fe - - c*) - 1)1 "fc 

1 


Pi(2:i) , qiiVi) 


-E/A {Xi-K*){yi-Mn ^ + 

1 •'ii-I L 


yi 


ds 


+ E / ‘ W - “ ')! * 

fLi Xi-K*d'^Xi^ yi-M*d^yi 


ds 


Now since, under both set of boundary conditions 

5x} 

5s 


[xi(0,t) - K*]^{0,t) = ML2,t) - -S^1^(^2,i) - 0 


(7.71) 
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(7.72) 


(7.73) 


(7.74) 


Hence ViK\M\C*) = 0 and V{x,y) < 0, if the conditions (7.66) ^ (7.69) are 
satisfied. Therefore h (x, y) is negative definite over x>0, y>0, i7>0 with respect to 
z* = A *, y* = M*, Rl = C*, proving the theorem. I 


Remark: On comparing the analysis of this chapter with chapter-6, we may conclude that 
the role of supplementary resource is to increase the level of nonuniform steady state at 
each point of the habitat. As before, the role of patchiness is destabilizing in presence of 
supplementary resource also, which we can note by comparing the cases of non-uniform 
and uniform steady states in a patchy habitat. 


7.4 Numerical Example 


Wti discuss here a numerical example, to study the behavior of the steady state solutions 
of tlu' aliove system and compare it with the case of a competing species without supple- 
menl.ary rc'sourcc in a two-patch habitat, as described in chapter-6. For this we consider 
the following particular form of functions: 

gi{ui) = ; ^ , ff(ui) = Si ^1 - , Pi{ui) = BiUi, qi{vi) = hVi, i = l,2 

also for simillJcity take, Ci = Ca = C. Then the steady state of the above system takes the 




9aiC , s 

-P (Oi - OLiUi) 

ai 


- 

9 ai( 3 iC 

-Vi 

CiUi + Ui 

ai 


= 0 (7.75) 
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n 

+ Vi 


mJ + — ft”<) 


-Ui 

, MAC' 

OiVi + — Vi 

- 

0-i 


= 0 


with reservoir boundary conditions 

ui(0) = xl\ U 2 {L 2 ) = xl\ and Ui(0) = y**, V 2 {L 2 ) = y”, 
and the continuity and flux matching conditions at the interface s = h, 

= IP.^(L0 = D.^(L0, and 

Ui{Li) = U2{Li),vi{Li) = V2{Li). 


(7.76) 


(7.77) 


(7.78) 

(7.79) 


We solve the system (7.75)-(7.79) numerically, for the following set of parameters, 
Li = 10 , Li = 20 , Du = 0 . 8 , Du = 0 . 8 , D 21 = 0 . 8 , D 22 = 0 . 8 , ri = 0.115, rj = 0.825, 

а, = 1.0, 02 = 1.0, 6 'i = 0.135, S 2 = 0.3, ei = 5.0B-05, 62 = 5.0F-05, bi = 4.0F-05, 62 = 

б. 0£^ ~ 05, h"i = 57.5, IC 2 = 165, Mi = 33.75, M 2 = 150, oi = 5.0F - 05, 02 = 2.0F - 05, 

C = GO, /A = 5.0£^ — 04, 72 = 7.0F — 04, = 0.4, and (j) = 0.8, and the plot of the steady 

state solutions are shown in Fig. 6 . 1 . By using the above values of the parameters, we get, 
x;* = 63.20 > x\ = 55.88, x^* = 176.38 > = 162.16, yj* = 35.92 > yl = 32.35, and 

y” = 181.18 > y 2 = 141.90. As in chapter-5, it is noted from the figure that, in presence 
of a renewable supplementary resource, also the steady state solutions are continuous, 
monotonic and the level of both the species are higher at each point of the habitat. 
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Figure 7.1: The steady state solutions for both the species, with and without a common 
supplementary resource for both the species 
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7.5 Summary 

In chapter-6, we have discussed the existence and stability of the non-uniform steady state 
corresponding to two competing species system in a two patch habitat with diffusion. In 
this chapter, the same ecological problem is modelled and analyzed when there is a non- 
diffusing renewable common resource for each of the competing population in both the 
patches. It has been assumed that both the competing population follow a general type of 
logistic growth in each of the patches, as in the previous chapter. Also, the self-renewable 
supplementary resource is logistically growing but non-diffusing. 

It has been shown, as in chapter-6, that there exists a positive, monotonic, continuous 
steady state solution with continuous matching conditions at the interface for both the 
species separately. The conditions for asymptotic stability in both linear and nonlinear 
cases under both set of boundary conditions have also been obtained. It has been further 
shown that in presence of supplementary resource the level of steady state distribution is 
correspondingly greater and the effect of patchiness is destabilizing even in presence of a 
common self-renewable supplementary resource. 
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